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The  simulation  of  flows  past  ships  poses  a  particular  challenge  because  of  the  size  of  the  ship  and  the 
resulting  large  values  of  dimensionless  parameters  such  as  Reynolds  number.  As  a  result,  the  only  simulation 
method  that  can  duplicate  the  required  parameters  is  computational  fluid  dynamics  (CFD).  However,  the  drawback 
of  CFD  at  these  large  scales  is  the  lack  of  data  available  for  validation  of  the  CFD  solutions.  To  provide  some 
insight  into  this  problem,  experimental  measurements  were  made  on  a  small-scale  model  and  a  large-scale  model. 
The  small-scale  model  was  0.286  m  (0.940  ft)  long  and  was  tested  in  a  water  tunnel  facility  located  at  Mississippi 
State  University.  A  method  of  obtaining  drag  values  from  wake  profile  measurements  was  developed  and  the 
measurements  were  made  using  a  laser  Doppler  velocimeter  (LDV)  system.  The  large-scale  model  was  7.03  m 
(23.0  ft)  long  and  was  tested  in  the  U.  S.  Navy  William  B.  Morgan  Large  Cavitation  Channel  (LCC). 
Computational  simulations  of  the  flows  at  both  of  these  scales  were  obtained  to  see  how  well  the  solutions  agreed 
with  the  experimental  results  over  a  very  large  change  in  length  scale.  This  was  of  particular  interest  since  the  factor 
of  the  increase  in  length  scale  in  going  from  the  LCC  model  to  the  full-scale  prototype  is  approximately  the  same  as 
the  factor  of  the  increase  in  length  scale  between  the  small-scale  and  large-scale  models.  The  drag  coefficients  of 
both  models  were  measured.  The  thrust  coefficients  for  a  propeller  on  the  large-scale  model  were  measured  over  a 
range  of  advance  ratio.  The  drag  coefficients  for  both  models  and  the  propeller  thrust  coefficients  for  the  large-scale 
model  were  obtained  from  the  computational  simulations.  The  experimental  and  CFD  drag  coefficients  were  in 
reasonable  agreement  for  the  small-scale  model.  For  the  large-scale  model,  the  computational  simulations  showed  a 
continuous  decrease  in  the  drag  coefficient  with  Reynolds  number,  typical  behavior  for  a  flow  in  which  the  drag  is 
dominated  by  skin  friction.  However,  the  large-scale  model  experiments  showed  that  the  drag  coefficient  leveled 
off  at  larger  Reynolds  number,  indicating  the  likely  presence  of  a  separated  region  that  would  cause  the  drag  to  be 
dominated  by  pressure  forces  rather  than  by  skin  friction.  LDV  measurements  of  the  flow  at  the  stem  of  the  large- 
scale  model  indicated  a  small  separated  region  with  reversed  flow.  Such  a  region  was  not  obvious  in  the 
computational  solutions.  The  computed  and  measured  propeller  thrust  coefficients  were  in  excellent  agreement.  It 
was  concluded  that  a  higher  resolution  grid  and  solution  was  probably  necessary  to  capture  fully  the  separated 
region  behind  the  submarine  and  accurately  predict  the  drag  behavior.  This  is  a  not  uncommon  problem  in  CFD, 
where  large  scale  geometries  at  large  Reynolds  number  require  large  grids  with  high  resolution  to  simulate  the  flow 
accurately. 
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b  Wake  width  scale 

Cp  Coefficient  of  drag 

Ci  Wake  width  parameter  obtained  from  curve  fits,  used  to  compute  drag  coefficient 
c„  Wake  velocity  deficit  parameter  obtained  from  curve  fits,  used  to  compute  drag  coefficient 
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1.  INTRODUCTION 


1.1  Statement  of  the  problem 

It  is  a  well-established  law  of  experimental  fluid  dynamics  that  in  order  to  achieve  dynamic  similarity 
between  prototypes  and  experimental  models,  the  values  of  certain  dimensionless  parameters  such  as  the  Reynolds 
number  have  to  be  matched  between  the  prototype  and  the  model.  It  is  also  a  well-established  fact  that  because  of 
the  large  length  scales  involved,  experiments  involving  ship  prototypes  rarely  achieve  complete  dynamic  similarity. 
It  is  usually  very  difficult  to  match  values  of  the  Reynolds  numbers,  and  almost  always  impossible  to  match 
Reynolds  and  Froude  numbers  simultaneously.  The  ship  study  problem  is  often  quoted  in  elementary  textbooks  on 
fluid  mechanics  as  an  example  of  incomplete  similarity  (e.g.  Fox  and  McDonald  1998). 

The  advent  of  computational  fluid  dynamics  (CFD)  solutions  in  recent  years  has  meant  that  it  is  now 
possible  to  simulate  flows  at  the  required  values  of  the  dimensionless  parameters,  allowing  ship  designs  to  be 
“tested”  before  construction  begins.  However,  the  issue  of  validation  of  the  computational  results  then  arises.  If  the 
new  ship  designs  are  different  enough  from  existing  models  and  prototypes,  there  are  no  data  available  by  which  the 
computational  simulations  can  be  validated.  While  CFD  has  undergone  significant  development  over  the  past  few 
years,  in  many  ways  it  is  still  in  a  developmental  stage  and  so  information  that  can  attest  to  the  validity  of  the 
solutions  is  welcomed. 

This  leads  to  the  purpose  of  the  current  study.  The  goal  was  to  study  experimentally  two  submarine 
models,  one  that  would  fit  inside  a  relatively  small  water  tunnel  located  at  Mississippi  State  University  (MSU),  and 
a  second  model  that  is  used  in  the  U.  S.  Navy’s  William  B.  Morgan  Large  Cavitation  Channel  (LCC)  located  in 
Memphis,  TN.  The  length  of  the  MSU  model  was  to  be  approximately  1  ft,  and  the  length  of  the  LCC  model  was  to 
be  approximately  23  ft,  so  that  there  would  be  over  an  order  of  magnitude  difference  in  size.  The  flows  past  these 
models  were  to  be  simulated  computationally,  in  order  to  see  how  well  the  computational  solutions  tracked  the 
measured  flows  over  large  changes  in  length  scale  and  hence  Reynolds  number.  In  particular,  the  drag  of  the 
submarine’s  bare  hull  and  the  thrust  provided  by  a  propeller  attached  to  an  appended  hull  were  to  be  measured  and 
computed,  and  the  results  compared.  The  goal  was  to  see  if  these  results  would  provide  some  information  as  to  how 
well  the  computational  solutions  would  predict  the  actual  flow  over  the  change  in  length  scale  involved  in  going 
from  the  LCC  model  to  a  full-scale  prototype. 
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1.2  Submarine  and  propeller  geometries 

As  just  mentioned,  there  were  essentially  three  studies  conducted  in  this  investigation  -  an  experimental 
study  using  a  small  scale  model  in  the  MSU  water  tunnel,  an  experimental  study  using  a  much  larger  scale  model  in 
the  LCC,  and  a  computational  study  of  both  model  geometries  and  a  full-scale  geometry.  The  common  threads  in 
these  three  different  studies  were  the  submarine  hull  geometry  and  the  propeller  geometry.  The  submarine  hull  and 
propeller  geometries  were  dictated  by  the  models  to  be  used  in  the  LCC,  so  the  LCC  hull  and  propeller  geometries 
will  be  discussed  at  some  length  here.  The  hull  geometry  was  a  standard  axisymmetric  hull  (DTMB  Model  5495-3). 
The  hull  outline  is  shown  in  Figure  1.1.  For  the  drag  studies,  the  hull  geometry  had  no  sail  or  appendages.  The 
geometry  included  a  “dummy  hub”  and  afterbody  geometry  that  corresponded  to  a  configuration  used  in  a  previous 
test  on  this  model  (see  Bridges  2004  and  the  discussion  below).  Fins  with  NACA  0012  airfoil  cross-sections  were 
added  to  the  stem  in  a  standard  cruciform  arrangement  for  the  configuration  used  in  the  thrust  studies.  A 
photograph  showing  the  arrangement  of  these  appendages  on  the  LCC  model  is  shown  in  Figure  1.2. 

The  propeller  geometry  used  in  these  studies  corresponded  to  the  David  Taylor  Model  Basin  Propeller  No. 
4381  (hereinafter  referred  to  as  P4381).  This  propeller  was  originally  designed  and  built  as  part  of  a  study  of  the 
effects  of  skew  on  marine  propellers  (Boswell  1971).  P4381  was  the  unskewed  propeller  in  this  series. 

Interestingly,  this  study  was  funded  by  several  commercial  shipping  companies,  and  so  the  propeller  design  itself 
and  the  results  of  the  study  were  unclassified.  The  fact  that  P4381  has  an  unclassified  geometry  and  that  the  actual 
model  propeller  still  exists  has  resulted  in  a  number  of  both  experimental  and  computational  studies  being 
performed  on  it,  including  the  current  study. 

The  actual  propeller  used  in  the  LCC  experiments  is  shown  in  Figure  1.3,  in  both  a  photograph  (from 
Hecker  and  Remmers  1971)  and  an  Autocad  IGS  file  (courtesy  Martin  Donnelly,  NSWCCD).  A  plan  drawing  of 
one  of  the  blades  (from  Boswell  1971)  is  also  included  in  the  figure.  The  actual  P4381  is  a  five-bladed  propeller 
with  a  diameter  of  1  ft,  an  expanded  area  ratio  of  0.725,  an  NACA  a  -  0.8  section  meanline,  and  an  NACA  66 
section  thickness  distribution  with  NSRDC  modifications  to  the  nose  and  tail  thicknesses.  The  design  advance 
coefficient  J  was  0.889,  and  the  design  thrust  loading  coefficient  Cn  was  0.534  (Boswell  1971).  Open-water  data 
for  the  thrust  coefficient  KT and  ten  times  the  torque  coefficient  \0Kq  are  shown  in  Figure  1.4  (data  from  Hecker  and 
Remmers  1971). 

P4381  was  chosen  for  the  LCC  tests  for  both  its  existence,  which  saved  the  cost  of  manufacturing  a  new 
propeller,  and  for  its  status  as  an  unclassified  geometry.  The  drawback  of  using  this  propeller,  however,  was  that  it 
was  not  designed  to  be  used  on  the  hull  model  that  was  used.  The  overall  propeller  diameter  D  was  smaller  than  that 
of  propellers  normally  attached  to  the  hull  model  used  in  the  LCC  experiments,  and  the  propeller  hub  was  not 
tapered.  The  hub  diameter  was  smaller  than  that  of  propellers  normally  attached  to  the  model  hull.  This  obstacle 
was  overcome  by  attaching  a  tapered  hull  extension  to  the  shaft  of  the  motor  in  the  model,  and  then  attaching  the 
propeller  to  this  hull  extension.  The  hull  extension  and  dummy  hub  are  shown  in  the  drawing  in  Figure  1.5  (the 
shaded  portion  of  the  figure  is  the  hull  extension  and  dummy  hub),  and  the  hull  extension  and  P4381  are  shown 
mounted  on  the  model  hull  in  the  photograph  in  Figure  1.5.  The  hull  extension  was  designed  such  that  it  rotated 
with  the  propeller.  The  hull  extension  continued  the  hull  slope  down  to  the  diameter  of  the  hub  on  P4381.  The 
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“dummy  hub,”  shown  in  the  drawing  in  Figure  1.5,  was  used  initially  for  calibration  purposes  prior  to  the  powered 
tests.  An  elliptic  tip  was  attached  behind  the  propeller,  as  shown  in  Figure  1.5.  The  powered  tests  were  conducted 
using  the  hull  extension,  propeller,  and  elliptic  tip.  The  drag  studies  were  conducted  using  the  hull  extension,  the 
dummy  hub,  and  the  elliptic  tip.  The  configuration  consisting  of  the  tapered  hull  extension,  P4381  and  the  elliptic 
tip  was  duplicated  for  the  computational  thrust  studies.  The  configuration  consisting  of  the  tapered  hull  extension, 
the  dummy  hub  and  the  elliptic  tip  was  duplicated  for  the  MSU  model  and  the  computational  hull  model  for  the  drag 
studies. 


1.3  Overview  of  report 

The  computational  studies  were  performed  by  a  graduate  student,  Allison  Cruthird  Cash,  utilizing  the 
resources  of  the  SimCenter  located  in  the  Engineering  Research  Center  at  Mississippi  State  University.  The 
experimental  study  conducted  at  MSU  was  performed  by  a  graduate  student,  John  Freudenthal,  utilizing  a  small 
water  tunnel  at  MSU.  The  theses  that  resulted  from  these  investigations  (Cash  2001;  Freudenthal  2002)  are  included 
in  this  report  as  Appendices  A  and  B,  respectively.  These  theses  form  essentially  self-contained  reports  of  the  work 
involved.  Therefore,  in  the  main  body  of  this  report,  the  significant  results  from  these  works  will  be  summarized  in 
individual  chapters,  and  the  reader  will  be  referred  to  the  appendices  for  the  details  of  the  studies.  Chapter  2  will 
summarize  the  computational  studies,  and  Chapter  3  will  summarize  the  experimental  studies  conducted  in  the  MSU 
water  tunnel.  The  experiments  involving  the  large-scale  model  in  the  LCC  were  conducted  by  the  principal 
investigator.  These  results  will  be  discussed  in  detail  in  Chapter  4,  and  then  the  results  of  the  three  different  studies 
will  be  compared  in  Chapter  5.  Conclusions  will  be  given  in  Chapter  6. 

A  word  concerning  the  experimental  study  conducted  at  MSU  is  in  order.  Initially,  as  stated  above,  the 
drag  of  a  bare  submarine  hull  and  the  thrust  of  a  powered  model  were  to  be  measured  in  both  experimental  facilities 
(MSU  and  LCC).  A  submarine  model  approximately  1 1  in  long  and  lin  in  diameter  and  matching  the  LCC  model 
geometry  was  constructed  for  the  MSU  tunnel.  As  discussed  previously,  the  propeller  that  was  to  be  used  on  the 
LCC  model  was  P438 1 .  When  the  propeller  was  actually  acquired,  it  was  discovered  that  the  propeller  diameter  was 
1  ft,  yielding  a  propeller-to-body  diameter  ratio  of  approximately  0.5  for  the  LCC  model.  This  would  have  implied 
a  propeller  diameter  of  0.5  in  for  the  MSU  model.  It  was  decided  that  it  would  be  extremely  difficult  to  manufacture 
an  accurate  model  of  the  propeller  at  that  scale,  and  so  the  decision  was  made  to  conduct  only  unpowered  studies  in 
the  MSU  water  tunnel  and  measure  only  the  drag  on  the  MSU  model  without  a  propeller.  These  results  will  be 
compared  to  the  computational  simulations  and  the  LCC  experimental  results.  The  powered  tests  were  conducted  in 
the  LCC  and  the  resulting  thrust  coefficient  measurements  will  be  compared  with  the  computational  results. 
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2.  RESULTS  OF  COMPUTATIONAL  SIMULATIONS 


2.1  Overview  of  chapter 

Because  of  the  availability  of  personnel  and  computer  time,  the  computational  study  (Appendix  A)  was 
fairly  extensive  and  actually  exceeded  the  scope  of  the  studies  described  in  Chapter  1 .  There  were  four  phases  to  the 
computational  investigation:  (1)  simulations  of  the  flow  past  the  bare  hull  geometry  at  the  three  different  length 
scales  discussed  in  Chapter  1  and  the  corresponding  drag  calculations;  (2)  simulations  of  an  isolated  Propeller  4381 
(so-called  “open  water”  studies);  (3)  simulations  of  powered  motion  with  P4381  mounted  on  the  hull  with  stem 
appendages;  and  (4)  simulations  of  the  submarine  crashback  maneuver  (forward  motion  and  reverse  propeller 
rotation).  The  computational  methods  and  other  computational  issues  will  be  summarized  here.  The  results  of  the 
first  three  phases  will  be  discussed  here  briefly,  and  phases  (1)  and  (3)  will  be  discussed  in  more  detail  in  Chapter  5, 
when  the  results  of  all  three  studies  are  brought  together.  The  reader  will  be  referred  to  Appendix  A  for  the 
discussion  of  the  crashback  simulations. 

2.2  Computational  methods 

The  flow  solver  used  in  the  computations  was  a  Reynolds-averaged  Navier-Stokes  solver  called 
“ ifNCLE ,”  for  “Unstructured  Unsteady  Computation  of  FieLd  Equations.”  It  is  capable  of  obtaining  solutions  for 
viscous,  incompressible,  high  Reynolds  number  flows  on  unstructured  grids  through  the  use  of  a  pseudo¬ 
compressibility  method.  For  more  details  on  this  solver,  see  Hyams  (2000)  and  Hyams  el  al.  (2000)*.  Two 
different  eddy-viscosity  turbulence  models  were  used,  a  one-equation  model  by  Spalart  and  Allmaras  (see  Spalart 
and  Allmaras  1992)  and  a  two-equation  model  referred  to  as  the  “q-co”  model  (see  Coakley  and  Hsieh  1985).  A 
composite  unstructured/structured  surface  grid  was  generated  for  the  hull  model  and  a  multielement  unstructured 
mesh  was  generated  for  the  volume  outside  of  the  hull.  See  pp.  15-17  of  Appendix  A  for  more  details  concerning 
the  grid  generation. 

The  simulations  were  conducted  on  a  super  cluster  of  256  dual  gigahertz  Pentium  III  computers.  All  of  the 
bare  hull  simulations  were  conducted  using  local  time  stepping  and  were  run  until  they  converged  to  steady-state 
solutions.  The  average  computation  time  for  the  bare  hull  simulations  was  approximately  10  hours.  The  bare  hull 
simulations  typically  converged  in  500  time  steps,  although  the  simulations  were  usually  run  to  1500  time  steps  to 


Both  Appendices  A  and  B  contain  lists  of  references.  Only  those  references  in  those  lists  which  are  referred  to  in  the  main 
body  of  the  report  are  included  in  the  list  of  references  for  the  main  body  of  the  report. 
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ensure  convergence.  The  isolated  propeller  simulations  were  run  with  local  time  stepping  for  approximately  400 
time  steps,  and  then  the  solver  was  switched  to  a  minimum  time  step  for  an  additional  600  to  700  time  steps. 
Convergence  was  usually  achieved  in  the  vicinity  of  700  total  time  steps,  with  a  computational  time  of  23  hours. 
The  convergence  histories  of  these  runs  can  be  seen  in  Figures  (A)3.4,  (A)3.5,  and  (A)3.8  .  Figure  (A)3.5  is  of 
interest  in  that  in  compares  the  convergence  histories  among  the  different  types  of  grids  and  turbulence  models  used 
in  the  bare  hull  studies.  This  figure  shows  that  all  of  the  different  combinations  converge  to  the  same  level  within 
approximately  700  time  steps. 

The  propelled  hull  simulations  were  a  bit  more  problematic.  The  solution  was  first  calculated  for  400  time 
steps  using  local  time  stepping.  Then,  the  stepping  was  switched  to  a  minimum  time  step.  The  rotating  propeller 
was  handled  by  locally  regenerating  the  grid  in  the  vicinity  of  the  propeller  for  each  time  step,  thus  simulating  an 
actual  rotating  propeller  on  a  fixed  body.  The  flow  was  calculated  for  800-1000  time  steps  until  it  reached  a 
periodic  state  with  a  very  small  peak-to-peak  amplitude,  as  shown  in  Figure  (A)3.11.  The  periodic  nature  of  the 
“converged”  solution  is  appropriate  for  the  thrust  force  produced  by  the  rotating  propeller.  The  total  computational 
time  for  each  of  these  solutions  was  typically  23  to  26  hours. 

2.3  Results 

The  results  of  the  computational  studies  for  the  drag  coefficient  of  the  bare  hull  are  shown  in  Figures  2.1 
and  2.2.  Figure  2.1  shows  the  results  for  just  the  LCC  model  studies,  comparing  the  two  different  turbulence 
models  used.  This  figure  shows  that  the  “q-ro”  model  consistently  predicted  a  higher  drag  coefficient  over  this  range 
of  Reynolds  number.  Figure  2.2  shows  the  results  for  all  three  ranges  of  Reynolds  number  examined.  The  MSU 
model  was  0.940  ft  long,  the  LCC  model  was  23.0  ft  long,  and  a  full-scale  prototype  was  assumed  to  be  384  ft  long, 
yielding  an  overall  Reynolds  number  range  for  the  simulations  of  7.7*105  <  ReL  <  1.6*  109.  The  computational 
simulations  for  the  MSU  model  and  the  prototype  were  conducted  using  only  the  “q-co”  model  and  the  results  for  all 
three  length  scales  seem  to  fall  consistently  on  a  single  curve  (note  the  logarithmic  scale  for  ReL).  The  continuous 
decrease  in  drag  coefficient  indicates  a  drag  dominated  by  viscous  shear  or  skin  friction  as  opposed  to  pressure  drag. 
The  results  for  the  drag  coefficients  will  be  discussed  in  more  detail  in  Chapter  5. 

The  results  of  the  computational  studies  of  the  isolated  propeller  (the  simulations  of  the  open-water  studies) 
are  shown  in  Figure  2.3.  This  figure  shows  the  computational  results  using  both  turbulence  models  for  the  propeller 
thrust  and  torque  coefficients  as  functions  of  advance  ratio  J.  These  results  are  compared  with  experimental  data 
from  Hecker  and  Remmers  (1971)  that  are  illustrated  in  Figure  1.4.  The  computational  results  are  in  very  good 
agreement  with  the  experimental  results  for  the  range  of  advance  ratio  studied.  It  should  be  noted  that  the  design 
advance  ratio  for  P4381  was  0.889  (Boswell  1971). 

The  results  for  the  simulations  of  powered  motion  using  P4381  and  an  appended  hull  are  shown  in  Figures 
(A)4.5  through  (A)4.8.  The  experimental  data  in  these  figures  are  from  the  LCC  tests  to  be  discussed  in  Chapter  4, 

*  The  figures  and  equations  in  Appendices  A  and  B  were  not  renumbered  when  these  theses  were  incorporated  into  this  report.  In 
the  main  body  of  this  report,  figure  and  equation  numbers  in  Appendices  A  and  B  will  be  preceded  by  an  “(A)”  or  “(B)”,  e.g. 
Figure  (A)3.4  is  Figure  3.4  in  Appendix  A.  A  figure  number  not  preceded  by  a  letter  will  refer  to  a  figure  in  the  main  body  of  the 
report. 
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so  a  full  discussion  of  these  results  will  deferred  until  Chapter  5.  It  is  sufficient  at  this  time  to  note  that  the 
computational  results  were  in  very  good  agreement  with  the  experimental  data.  It  was  discovered  during  these 
studies  that  the  computations  could  not  be  conducted  for  the  powered  configuration  using  the  “q-co”  turbulence 
model.  The  turbulence  model  became  unstable  and  diverged  in  a  relatively  few  time  steps.  The  suspected  culprit 
was  the  surface  grid  on  the  propeller.  The  propeller  tips  were  very  fine  and  required  a  large  number  of  points  for 
adequate  resolution.  The  grid  generator  occasionally  created  minute  “divots”  in  the  surface  topology,  especially 
near  sharp  edges.  The  “q-co”  turbulence  model  was  especially  sensitive  to  such  disturbances,  hence  the  instability 
when  these  cases  were  attempted.  Since  there  was  very  little  difference  between  the  results  for  the  “q-co”  and 
Spalart-Allmaras  turbulence  models,  the  computational  studies  were  continued  from  this  point  using  only  the 
Spalart-Allmaras  model. 

Additional  effects  of  Reynolds  number  are  explored  in  Section  4.3  of  Appendix  A.  Figure  (A)4.20 
illustrates  the  continuous  reduction  in  the  skin  friction  coefficient  for  the  bare  hull  as  Reynolds  number  increases, 
which  is  what  one  normally  expects.  Figure  (A)4.21  compares  the  flow  in  the  exit  plane  of  the  isolated  propeller  at 
two  extremes  of  the  Reynolds  number  range.  This  figure  shows  regions  of  higher  velocity  near  the  propeller  blades 
at  the  lower  Reynolds  number.  These  regions  disappear  at  the  higher  Reynolds  number.  Figure  (A)4.22  compares 
the  overall  isolated  propeller  flow  field  at  these  two  values  of  the  Reynolds  number  and  illustrates  the  fact  that  the 
propeller  has  a  greater  influence  on  the  surrounding  flow  field  at  the  lower  Reynolds  number.  Figure  (A)4.23 
compares  the  flow  in  the  region  of  the  stem  planes  at  the  two  values  of  the  Reynolds  number  for  the  appended  hull 
without  propeller,  this  time  showing  a  greater  influence  of  the  stem  planes  on  the  surface  pressure  distribution  at  the 
higher  Reynolds  number.  The  conclusion  drawn  from  these  studies  was  that  for  the  hull  and  propeller  flow  fields, 
the  Reynolds  number  had  a  localized  effect  on  flow  field  characteristics. 

Sections  4.4  and  4.5  of  Appendix  A  deal  with  hull/propulsor  interactions  and  the  crashback  maneuver, 
respectively.  While  these  results  are  of  interest,  these  studies  were  outside  the  immediate  scope  of  the  current  study. 
The  reader  is  referred  to  Appendix  A  for  the  details  of  these  results. 
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3.  EXPERIMENTAL  STUDIES  ON  SMALL-SCALE  MODEL 


3.1  Overview  of  chapter 

The  experimental  studies  on  the  small  scale  model  (Appendix  B)  were  conducted  in  a  water  tunnel  located 
at  Mississippi  State  University  (MSU)  using  a  laser  Doppler  velocimeter  (LDV)  system  to  obtain  wake  velocity 
profiles.  A  previously-developed  method  for  extracting  drag  data  from  two-dimensional  wake  velocity  profiles  by 
Dimotakis  (1977)  was  extended  to  axisymmetric  flows.  The  implementation  of  this  method  allowed  the 
determination  of  the  drag  of  the  model  solely  from  LDV  wake  velocity  measurements.  The  method  will  be 
summarized  briefly,  the  water  tunnel  and  model  will  be  described  briefly,  and  the  results  of  the  measurements  will 
be  summarized. 

3.2  Wake  profile  drag  method 

The  method  by  Dimotakis  began  by  integrating  the  momentum  equation  around  a  control  volume 
surrounding  a  two-dimensional  body.  This  procedure  was  repeated  for  an  axisymmetric  body,  as  shown  in  Figure 
(B)2.1.  The  resulting  equation  for  the  drag  coefficient  is  given  by  Equation  (B)(2.46).  This  equation  is  written  in 
terms  of  radial  integrals  of  mean  and  fluctuating  streamwise  velocity  components  and  of  the  mean  pressure 
coefficient.  The  turbulent  wake  flow  was  then  assumed  to  be  self-preserving,  and  the  appropriate  forms  of  the  mean 
and  fluctuating  velocities  and  the  mean  pressure  were  assumed,  as  shown  in  Equations  (B)(2.48)  through  (B)(2.50). 
The  result  for  the  drag  coefficient  in  terms  of  the  associated  length  and  velocity  scales  and  integrals  of  the 
dimensionless  self-similar  distributions  of  velocity  and  pressure  is  given  by  Equation  (B)(2.60).  The  length  and 
velocity  scales  were  assumed  to  vary  according  to  standard  power  law  relations  so  that  the  length  scale  increased  as 
xvi  and  the  velocity  scale  decreased  as xm ,  in  the  specific  forms 
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In  these  formulae,  u0  and  IB  are  the  velocity  and  length  scales,  respectively,  for  the  self-preserving  wake,  TJ0  is  the 
velocity  outside  the  wake,  x  is  the  downstream  distance  measured  from  the  trailing  edge  of  the  model,  x0  is  a  virtual 
origin  determined  from  fits  to  the  experimental  data,  d  is  the  model  diameter,  and  c„  and  c,  are  constants  of 
proportionality  determined  from  fits  to  the  experimental  data.  When  these  variations  are  included  in  Equation 
(B)(2.60),  and  the  results  are  extrapolated  to  very  large  x,  the  terms  involving  the  fluctuating  velocity  component 
and  the  mean  pressure  coefficient  become  negligible,  and  the  result  is  given  by  Equation  (B)(2.64).  This  equation  is 
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further  simplified  if  the  body  diameter  d  is  very  small  compared  to  the  radius  R  of  the  control  volume.  The  result, 
given  by  Equation  (B)(2.65),  is 


8c„c/2 

In  2 


The  method  was  implemented  as  follows.  The  streamwise  velocity  component  profile  in  the  radial 
direction  was  obtained  at  several  locations  downstream  of  the  model.  Each  of  these  profiles  were  fit  to  the  function 


u(z )  =  U0-u0  exp 
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using  the  “genfif  ’  routine  in  Mathcad.  In  this  formula,  u(z)  is  the  measured  streamwise  velocity  component,  and  z  is 
the  transverse  position  coordinate.  The  variable  z0  is  the  nominal  wake  centerline,  and  l0  is  the  wake  half-width 
defined  as  the  point  such  that  when  z  -  zo  =  /«,  the  wake  velocity  deficit  is  half  its  maximum  value  on  the  wake 
centerline.  The  parameters  U0,  u0,  Zo,  and  l0  were  obtained  from  the  fits.  The  axial  variations  of  the  parameters  u0 
and  I0  were  then  fit  to  the  functions  of  x  given  above,  again  using  the  “genfit”  routine  in  Mathcad,  and  the  constants 
c„  and  c,  were  obtained  from  the  fits  and  used  to  calculate  CD.  The  resulting  value  was  adjusted  to  take  into  account 
the  pressure  gradient  in  the  empty  test  section.  The  correction  for  the  contribution  to  the  apparent  drag  force  caused 
by  the  nonzero  pressure  gradient  was 


D  =  —  V 

"  dx  body 


where  Vhnjy  is  the  volume  of  the  model.  This  is  the  so-called  “buoyancy  drag”  correction,  which  is  exact  if  the 
pressure  gradient  is  constant.  Data  from  the  thesis  by  Lindsey  (1998)  were  used  to  compute  the  test  section  pressure 
gradient  and  its  associated  uncertainty. 

The  advantage  of  the  wake  profile  method  is  that  a  significant  amount  of  information  goes  into  the 
calculation  of  the  drag  coefficient.  However,  the  drawback  is  that  it  takes  a  significant  amount  of  time  to  acquire  the 
information.  This  limited  the  number  of  cases  that  could  be  tested  in  the  time  available  for  the  experiments. 


3.3  Experimental  apparatus 

The  MSU  water  tunnel  has  an  axisymmetric  test  section  that  has  a  diameter  of  1 1  in  and  a  length  of 
approximately  4  ft.  The  top  test  section  speed  is  approximately  30  ft/s.  Details  of  the  design  and  construction  of 
this  facility  may  be  found  in  the  report  by  Wells  (1964).  This  facility  was  extensively  renovated  in  the  late  1990s, 
and  the  results  of  this  renovation  and  the  initial  tests  in  the  renovated  facility  may  be  found  in  the  thesis  by  Lindsey 
(1998;  see  also  the  Appendix  to  Bridges,  2001).  Further  work  on  the  facility  was  done  for  the  current  study;  see 
Section  3.1  of  Appendix  B  for  further  discussion  of  the  facility  and  the  recent  work  done  on  the  facility. 

The  velocity  measurement  system  used  was  a  Dantec  two-component  FiberFlow  LDV  system.  Because  the 
test  section  of  the  water  tunnel  is  axisymmetric,  the  beams  had  to  enter  the  test  section  perpendicular  to  the  test 
section  wall  and  traverse  a  diameter  of  the  test  section.  This  meant  that  only  one  pair  of  beams  could  be  used  and 


9 


hence  only  one  velocity  component,  the  streamwise  component,  could  be  measured.  For  more  details  of  the  LDV 
system  and  its  calibration,  see  Sections  3.2  and  3.3  of  Appendix  B. 

As  stated  earlier,  the  submarine  model  was  11.278  in  (286.47  mm)  long.  The  details  of  the  model  are 
discussed  in  Section  3.4  of  Appendix  B,  and  the  model  as  mounted  in  the  MSU  water  tunnel  is  shown  in  Figure 
(B)3.5. 

3.4  Experimental  results 

Results  for  the  velocity  profiles  in  the  empty  test  section  are  shown  in  Figures  (B)4.1  through  (B)4.5. 
These  results,  particularly  Figure  (B)4.5,  show  some  repeatable  variations  in  the  streamwise  velocity  distributions 
that  were  believed  to  be  due  to  irregularities  in  the  Plexiglas  wall  of  the  test  section,  possibly  caused  by  distortions 
of  the  wall  when  the  tunnel  was  pressurized.  Lindsey  (1998)  had  noted  such  distortions  in  his  work.  These  will  be 
commented  on  again  below.  Figures  (B)4.6  and  (B)4.7  show  the  distribution  of  the  RMS  component  of  the 
streamwise  velocity  across  the  test  section.  These  show  a  free  stream  turbulence  level  of  slightly  over  1%,  which  is 
higher  than  desired  but  deemed  reasonable  for  the  tests  being  conducted.  Figures  (B)4.9  through  (B)4.12  show  the 
total  uncertainty  levels  estimated  for  the  empty  test  section  measurements.  Figure  (B)4.12,  in  particular,  places  the 
uncertainty  levels  in  the  proper  perspective  and  shows  that  the  total  uncertainties  in  the  empty  test  section 
measurements  were  relatively  small. 

Section  4.3  of  Appendix  B  contains  the  wake  profile  results.  It  took  essentially  one  day  of  effort  to  acquire 
each  set  of  velocity  profiles.  The  mean  velocity  profiles  and  turbulence  intensity  profiles  are  shown  in  Figures 
(B)4.14  through  (B)4.23.  The  profiles  seem  to  indicate  a  “drift”  in  the  profile  location.  This  apparent  shift  was 
caused  by  a  slight  misalignment  between  the  test  section  centerline  and  the  LDV  axial  traverse  axis,  which  was 
determined  from  the  measurements  to  be  less  than  0.7  deg.  This  misalignment  was  the  reason  for  choosing  the 
centerline  position  z0  as  one  of  the  fit  parameters.  Figures  (B)4.24  through  (B)4.26  show  representative  velocity 
profiles  that  have  been  corrected  for  the  centerline  “shift”  and  include  the  uncertainties  in  the  individual  velocity 
measurements.  The  typical  uncertainty  in  a  velocity  measurement  was  3%  of  U0  at  a  given  profile  and  the  typical 
uncertainty  in  the  z  location  was  0.675  mm. 

Figures  (B)4.27  and  (B)4.28  show  the  velocity  profiles  plotted  in  similarity  variables.  The  positive  values 
on  the  abscissa  correspond  to  the  side  of  the  wake  closest  to  the  LDV  system.  As  the  LDV  measurement  volume 
moved  across  the  test  section,  the  data  acquisition  rate  and  the  signal  quality  dropped  to  the  point  where  it  was  just 
possible  to  measure  the  complete  wake  on  both  sides,  but  not  to  completely  span  the  test  section.  This  is  why  the 
profiles  are  asymmetric  in  z  with  respect  to  the  number  of  data  points  obtained.  There  were  also  some  apparent 
asymmetries  in  the  wake  measurements  themselves  which  are  apparent  in  the  figures.  The  reasons  for  these 
asymmetries  were  not  discovered.  This  issue  will  be  discussed  briefly  below. 

The  results  for  the  axial  variations  of  velocity  deficit  and  wake  width  are  shown  in  Figures  (B)4.30  and 
(B)4.31,  respectively.  It  should  be  noted  that  there  are  no  uncertainty  bars  on  these  figures.  The  uncertainty 
analyses  were  conducted  for  the  individual  velocity  measurements,  as  discussed  in  Section  4.2  of  Appendix  B  and 
shown  in  Figures  (B)4.24  through  (B)4.26,  as  just  mentioned.  As  a  result  of  time  constraints  on  the  completion  of 
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Appendix  B,  the  uncertainty  analysis  on  the  axial  variation  fits  and  the  final  calculation  of  CD  were  not  included  in 
Appendix  B.  That  uncertainty  analysis  will  be  discussed  here. 

The  problem  was  to  take  the  uncertainties  in  u  and  z  and  compute  the  uncertainties  in  the  parameters  ulh  b, 
cu  and  ci  arising  from  the  curve  fits*,  and  then  use  these  uncertainties  to  compute  the  uncertainty  in  the  drag 
coefficient  CD.  This  was  accomplished  using  a  Monte  Carlo  method,  similar  to  the  procedure  outlined  in  Coleman 
and  Steele  (1999).  The  typical  uncertainties  in  u  and  z  were  randomly  distributed  over  the  velocity  profile  data, 
using  the  normal-distribution  random-number  generator  in  Mathcad  with  a  mean  of  0  and  a  standard  deviation  of 
0.368  (=  1/e).  The  values  of  V0,  u0,  zn,  and  b  were  then  computed  using  the  “genfit”  function.  This  procedure  was 
repeated  10,000  times.  The  resulting  uncertainties  in  u0  and  b,  taken  as  twice  the  standard  deviations  of  u0  and  b 
from  the  10,000  results,  were  on  the  order  of  5%  and  6%  of  the  mean  values  of  u0  and  b,  respectively.  The 
distributions  of  u0/U0  and  b/D  with  x/D  are  shown  in  Figures  3.1  and  3.2,  respectively.  These  are  Figures  (B)4.35 
and  (B)4.36  repeated  with  the  uncertainty  bars  added  (only  the  experimental  values  used  in  the  fits  to  obtain  cu  and  C; 
are  shown  in  Figures  3.1  and  3.2).  As  can  be  seen  in  these  figures,  the  uncertainty  in  b/D  was  somewhat  larger  than 
the  uncertainty  in  uo/U0 ,  and  the  uncertainty  in  b/D  tended  to  grow  in  the  streamwise  direction,  reaching  a  maximum 
of  approximately  7%  of  the  mean  value  of  b  at  x/D  =  28.  The  individual  values  of  b/D  tended  to  show  more  scatter 
than  the  corresponding  values  of  un  (see  e.g.  Figure  4.31  of  Appendix  B). 

The  uncertainties  in  uJUo  and  b/D  were  then  randomly  distributed  over  the  values  of  u0/U0  and  b/D  using 
the  same  normal-distribution  random-number  generator  as  before  and  the  values  of  c„  and  c;  were  computed  using 
the  “genfit”  function.  This  procedure  was  also  repeated  10,000  times.  The  resulting  uncertainties  in  cu  and  c,,  taken 
as  twice  the  standard  deviations  of  cu  and  c,  from  the  10,000  results,  were  4%  and  15%  of  the  mean  values  of  c„  and 
Ci.  The  larger  uncertainties  in  b/D  were  reflected  in  the  relatively  large  uncertainty  in  C/.  The  relatively  large 
uncertainty  in  c/  was  a  bit  surprising,  since  the  associated  uncertainties  in  the  width  parameter  b/D  were  typically  of 
the  order  of  6%,  reaching  the  peak  value  of  7%  only  at  the  most  downstream  location.  It  had  been  noted  in 
performing  the  curvefits  that  the  values  of  c„  and  particularly  C/  were  sensitive  to  small  changes  in  the  data.  This  is 
probably  due  to  the  power-law  nature  of  the  function  that  was  being  used,  and  the  fact  that  a  fit  to  a  power  law  in  x 
was  being  attempted  over  a  relatively  small  range  of  x.  Using  the  values  of  the  uncertainties  in  cu  and  c/,  the 
corresponding  uncertainty  in  the  drag  coefficient  was  found  to  be  A CD  ~  ±0.06,  or  30%  of  the  value  CD  =  0.195 
computed  using  the  procedure  described  above.  The  relatively  large  uncertainty  in  C/  was  magnified  even  further 
because  of  the  appearance  of  the  square  of  c:  in  the  formula  for  the  drag  coefficient.  It  is  believed  that  a  refinement 
in  LDV  techniques  would  reduce  the  uncertainties  in  the  results  to  acceptable  levels.  The  contribution  to  the 
uncertainty  from  the  uncertainty  in  the  pressure  gradient  term  was  negligible,  due  to  the  small  contribution  of  the 
pressure  gradient  term  to  the  drag  coefficient  itself. 

As  discussed  in  Section  4.3.3  of  Appendix  B,  the  value  of  CD  =  0.195  was  obtained  by  averaging  the  values 
of  Uo  and  b  as  functions  of  x/D  found  over  three  days  of  testing,  during  which  the  tunnel  pressure  was  allowed  to 
bleed  down  over  time.  Over  two  more  days  of  testing,  the  tunnel  pressure  was  maintained  at  a  high  value  and  a 

*  Through  an  inadvertent  switch  in  notation,  the  wake  width  parameter  is  denoted  by  “b”  instead  of  “/0”  in  Chapter  4  of  Appendix 
B.  In  the  discussion  to  follow,  the  notation  change  will  be  maintained,  so  that  the  current  discussion  is  consistent  with  the  usage 
in  Appendix  B. 
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value  of  Cd  =  0.201  was  obtained  by  averaging  the  values  of  ut)/Uo  and  b/D  found  during  these  tests.  The 
uncertainty  analysis  just  presented  was  performed  only  for  the  first  set  of  data,  due  to  the  computational  effort 
involved.  These  results  should  be  representative  of  the  results  as  a  whole. 

As  noted  above,  there  were  some  unexplained  anomalies  in  the  data.  Figure  (B)4.5  shows  the  repeatable 
velocity  variation  across  the  empty  test  section  which  was  believed  to  be  due  to  distortion  of  the  test  section  wall. 
This  figure  also  shows  an  apparent  overall  shift  in  the  velocity  of  approximately  0. 1  m/s.  This  was  believed  to  be 
due  to  some  anomaly  in  the  LDV  electronics  which  only  occurred  sporadically,  usually  at  relatively  long  time 
intervals,  but  occasionally  during  velocity  surveys.  Finally,  (B)4.27  and  (B)4.28  display  an  apparent  asymmetry  in 
the  wake,  which  was  also  thought  to  be  due  to  test  section  wall  distortion  or  possibly  a  misalignment  of  the  LDV 
traverse  in  the  direction  perpendicular  to  the  flow.  A  series  of  experiments  was  to  be  conducted  to  try  to  identify  the 
sources  of  these  various  anomalies.  The  first  was  the  apparent  shift  in  the  velocity  values.  A  long-duration  test  of 
the  test  section  speed  ahead  of  the  model  was  started.  The  results  are  shown  in  Figure  (B)4.34,  which  shows  one 
such  apparent  shift  at  a  time  of  approximately  12  minutes.  While  this  test  was  being  conducted,  the  water  tunnel 
ruptured,  and  there  was  not  enough  time  to  repair  and  continue  the  tests  before  the  project  ended. 
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4.  EXPERIMENTAL  STUDIES  ON  LARGE-SCALE  MODEL 


4.1  Overview  of  chapter 

The  large-scale  studies  were  conducted  in  the  Large  Cavitation  Channel,  as  discussed  earlier. 
Measurements  of  drag  and  LDV  surveys  of  the  velocity  field  at  the  stem  were  made  on  a  bare  hull  model  with  the 
“dummy  hub”  geometry  discussed  in  Chapter  1 .  Measurements  of  propeller  thrust  force  were  made  over  a  range  of 
advance  ratios.  Measurements  of  the  boundary  layer  at  four  axial  locations  were  also  made.  The  experimental 
apparatus  and  procedure  will  be  described  and  then  the  results  discussed  in  detail. 

4.2  Experimental  apparatus  and  procedure 

4.2.1  Facility 

The  experiments  to  be  described  in  this  chapter  were  conducted  in  the  William  B.  Morgan  Large  Cavitation 
Channel  (LCC),  located  in  Memphis,  Tennessee.  This  facility  is  part  of  the  Carderock  Division  of  the  Naval  Surface 
Warfare  Center  (NSWCCD)  and  was  first  made  operational  in  1991 .  The  LCC  has  a  test  section  that  is  3.05  m  (10 
ft)  high,  3.05  m  (10  ft)  wide,  and  12.2  m  (40  ft)  long.  The  maximum  test  section  speed  is  approximately  18  m/s  (35 
kts),  and  the  pressurization  range  for  the  test  section  extends  from  3.5  kPa  (0.05  psia)  to  415  kPa  (60  psia).  The 
freestream  turbulence  level  is  less  than  0.5%.  Details  of  the  design,  construction,  and  operation  of  this  facility  may 
be  found  in  Etter  and  Wilson  (1992;  1993).  More  recently  Park,  Cutbirth  and  Brewer  (2002)  have  done  an  extensive 
study  of  the  performance  characteristics  of  the  LCC  and  have  documented  in  detail  the  flow  temporal  and  spatial 
uniformity  and  the  turbulence  levels  in  the  facility,  as  well  as  the  data  acquisition  equipment  and  procedures  for 
velocity  measurements  in  the  LCC  test  section.  According  to  Park  et  al.  (2002),  the  temporal  stability  of  the  LCC 
test  section  flow  speed  is  ±0.15%  for  test  section  velocities  between  0.5  and  18  m/s.  The  test  section  velocity  is 
spatially  uniform  to  within  ±0.60%  for  velocities  between  3  and  16  m/s.  The  turbulence  level  in  the  facility  is 
between  0.2%  and  0.5%  for  test  section  velocities  from  0.5  to  15  m/s,  with  no  harmonics  appearing  in  the  power 
spectra  of  the  velocity  signals.  For  all  of  the  tests  conducted  during  this  study,  the  tunnel  pressure  was  set  at  30  psig 
at  the  test  section  top. 

4.2.2  Model 

As  discussed  in  Chapter  1 ,  the  model  used  was  a  standard  axisymmetric  hull  form  (DTMB  5495-3).  The 
hull  contour  is  shown  in  Figure  1.1.  For  the  drag  studies,  the  configuration  included  the  “dummy  hub”  that  is 
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represented  in  the  top  portion  of  Figure  1.5  and  did  not  include  stem  planes.  For  the  powered  studies,  P4381  and  the 
stern  planes  were  added  as  shown  in  the  lower  portion  of  Figure  1.5. 

4.2.3  Instrumentation  and  data  reduction 

The  forces  on  the  hull  model  were  measured  using  a  AMTI  6-component  dynamometer.  This  device  is  an 
internal  force  balance  that  was  attached  between  the  strut  and  the  model.  It  was  set  up  to  measure  forces  in  all  three 
directions  and  moments  about  all  three  axes.  The  dynamometer  was  calibrated  by  a  member  of  the  NSWCCD  staff 
who  oversaw  the  setup  and  the  initial  use  of  the  balance  in  the  model  at  the  LCC.  The  dynamometer  voltages  were 
converted  to  forces  and  moments  using  supplied  calibration  data.  The  systematic  errors  for  this  dynamometer  are 
shown  in  Table  4.1  (the  notation  “2o”  indicates  the  value  of  twice  the  standard  deviation). 


Table  4.1  Estimates  of  systematic  error  for  AMTI  dynamometer 


Fx 

Fy 

F 

Mx 

My 

Mz 

%  of  full  scale,  max 

0.18 

0.17 

0.6 

0.59 

0.19 

0.65 

%  of  full  scale,  min 

-0.11 

-0.28 

-0.75 

-0.49 

-0.21 

-0.67 

2o  (%  of  full  scale) 

0.08 

0.12 

0.33 

0.34 

0.12 

0.43 

full  scale,  eng.  units 

5000  lbs 

5000  lbs 

5000  lbs 

2500  ft-lbs 

25000  ft-lbs 

5000  ft-lbs 

2o,  eng.  units 

4  lbs 

6  lbs 

16.5  lbs 

8.5  ft-lbs 

30  ft-lbs 

21.5  ft-lbs 

The  forces  on  the  propeller  itself  were  measured  with  a  MicroCraft  6-componcnt  propeller  dynamometer. 
This  dynamometer  was  a  multi-component  cylindrical  assembly  that  attached  to  the  propeller  on  its  upstream  face. 
The  dynamometer  is  shown  in  the  upper  portion  of  Figure  1.5,  just  ahead  of  the  hub  extension.  Strain  gauges  were 
glued  to  12  flexures  in  order  to  measure  three  orthogonal  forces  and  three  moments.  The  gauges  were  arranged  in 
Wheatstone  bridge  circuits  with  two  to  four  arm  circuits  connected  in  parallel  for  each  channel.  The  placement  of 
multiple  bridges  was  intended  to  cancel  interactions  from  other  forces  and  moments.  This  particular  design  had 
small  interactions  (<2%)  except  for  bending  moment  interactions  on  the  orthogonal  side  forces  of  about  41%.  Care 
had  to  be  taken  to  properly  zero  the  balance  so  as  to  not  artificially  remove  possible  interactions. 

A  6  xl2  calibration  matrix  was  used  to  calculate  the  forces  and  moments  from  the  six  input  voltages.  This 
matrix  provided  first  and  second  order  interactions  of  the  various  forces  and  moments  and  yielded  a  dynamometer 
error  of  less  than  1%  of  each  channel’s  full  scale  value.  The  multi-component  construction  of  this  dynamometer 
caused  hysteresis  effects  that  were  the  chief  source  of  error.  The  estimates  of  the  systematic  error  for  this 
dynamometer  are  shown  in  Table  4.2. 

The  strain  gage  signals  passed  through  slip  rings  to  Vishay  Model  2310  Instrumentation  Amplifiers  outside 
the  tunnel.  These  amplifiers  provided  excitation  voltage,  bridge  zeroing,  gain,  and  filtering  of  the  signals.  For  these 
tests,  all  filters  were  turned  off,  resulting  in  a  wideband  setting  during  the  tests.  The  conditioned  signals  were  then 
passed  to  the  computer  A/D  inputs.  Due  to  a  limited  number  of  slip  ring  connections  it  was  necessary  to  group  some 
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of  the  bridge  excitation  lines  to  the  dynamometer.  Due  to  multiple  bridges  connected  in  parallel  the  current  draw 
along  these  excitation  lines  could  result  in  a  voltage  drop.  As  the  dynamometer  was  calibrated  with  individual 
excitation  for  each  force  or  moment  channel,  additional  calibration  checks  were  conducted  to  adjust  individually  the 
gain  of  each  channel  to  account  for  the  lower  excitation  voltage.  The  details  of  these  corrections  are  not  discussed 
here  except  to  state  that  the  corrections  adequately  reduced  any  error  introduced  to  levels  quoted  from  calibration 
data. 


Table  4.2  Estimates  of  systematic  error  for  MicroCraft  propeller  dynamometer 


Component 

Load  limit 

Error  as  %  of 
full  scale  (2a) 

fx  (thrust) 

2000  lbs 

0.40 

fy  (side  force) 

500  lbs 

0.55 

f  (normal  force) 

500  lbs 

0.68 

mx  (torque) 

500  ft-lbs 

0.17 

my  (pitch) 

500  ft-lbs 

0.17 

mz  (yaw) 

500  ft-lbs 

0.02 

The  “2a”  values  in  Tables  4.1  and  4.2  were  used  as  the  systematic  uncertainties  in  computing  the 
uncertainties  in  the  forces  measured  by  the  body  and  propeller  dynamometers.  The  propeller  and  body 
dynamometer  channels  were  sampled  at  a  rate  of  100  Hz  for  a  period  of  10  sec,  yielding  1000  samples  per  channel. 
Twice  the  standard  deviation  of  the  samples  acquired  during  the  tests  was  used  as  the  measure  of  the  random 
uncertainty. 

Measurements  of  the  velocity  field  at  the  model’s  stem  were  obtained  using  the  LCC’s  Dantec  laser 
Doppler  velocimeter  (LDV)  system.  This  system  is  described  in  detail  in  the  report  by  Park,  Cutbirth  and  Brewer 
(2002).  The  system  consists  of  4  Dantec  BSA  57N1 1  signal  processors,  3  fiber  optic  probes,  two  Spectra  Physics  6 
W  argon-ion  lasers,  a  Dantec  3-D  traverse,  and  Dantec  Flow  software.  The  LDV  system  was  calibrated  using  a 
rotating  disc.  The  resulting  uncertainty  in  the  calibration  is  reported  by  Park  et  al.  (2002)  to  be  less  than  0.018  m/s 
at  a  disc  speed  corresponding  to  a  tunnel  speed  of  15  m/s.  Based  on  the  results  presented  in  Park  et  al.,  the  overall 
uncertainty  in  individual  LDV  measurements  of  the  velocity  was  assumed  to  be  0.3%.  This  value  was  used  in  the 
calculation  of  the  force  coefficients,  Reynolds  numbers,  and  advance  ratios. 

The  density  and  viscosity  of  the  water  were  obtained  by  measuring  the  temperature  of  the  water  and  using 
curve  fits  to  compute  the  density  and  viscosity.  Calibrations  of  the  temperature  transducers  have  shown  an 
uncertainty  of  ±  2.0  deg  C,  caused  by  a  drift  in  the  calibrations  of  the  transducers.  Since  the  sensitivity  coefficients 
at  a  normal  operating  temperature  of  27  deg  C  are  0.028  %/deg  C  and  2.22  %/deg  C  for  the  density  and  viscosity, 
respectively,  the  uncertainty  in  the  temperature  leads  to  uncertainties  in  the  density  and  viscosity  of  0.056%  and 
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4.44%,  respectively*.  These  values  were  used  as  the  systematic  uncertainties  for  the  density  and  viscosity.  Because 
the  systematic  uncertainty  in  the  temperature  was  much  larger  than  random  uncertainties  recorded  in  the 
measurements  (typically  0.002  deg  C),  the  total  uncertainties  in  the  density  and  viscosity  were  essentially  equivalent 
to  the  systematic  uncertainties. 

The  uncertainty  in  the  propeller  speed  was  1  RPM  or  0.0167  s  '.  The  uncertainty  in  the  propeller  diameter 
was  assumed  to  be  0.0625  in  or  0.5%  of  the  propeller  diameter.  The  uncertainty  in  the  hull  model  length  was 
assumed  to  be  1.5  in  or  0.54%  of  the  overall  model  length.  The  uncertainty  in  the  hull  model  diameter  was  assumed 
to  be  0.125  in  or  0.51%  of  the  model  diameter. 

The  Reynolds  number  based  on  length,  ReL,  was  computed  from 


ReL 


VL 


V 


where  V  is  the  free  stream  or  tunnel  speed,  L  is  the  body  length,  and  v  is  the  kinematic  viscosity.  The  advance  ratio 
J  was  computed  from 

.7  =  22 

nD 

where  V  is  the  free  stream  or  tunnel  speed,  n  is  the  propeller  speed  in  revolutions  per  second,  and  D  is  the  propeller 
diameter.  The  drag  coefficient  CD  was  computed  from 


C°~  1 


n 


—  pV  —D. 

2  4 

where  Fx  is  the  axial  force  from  the  body  dynamometer,  p  is  the  fluid  density,  Dh  is  the  hull  diameter,  and  kD2I4  is 
the  projected  frontal  area  of  the  submarine  hull.  The  thrust  coefficient  KT  was  computed  from 

/, 


Kt  = 


2  r>4 


pn  D 


where  fx  is  the  axial  force  measured  by  the  propeller  dynamometer,  p  is  the  fluid  density,  and  D  is  the  propeller 
diameter.  Typically  the  drag  coefficient  was  plotted  as  a  function  of  Reynolds  number  and  the  thrust  coefficient  was 
plotted  as  a  function  of  advance  ratio. 


4.3  Results 

4.3.1  Drag 

The  results  for  the  drag  coefficient  as  a  function  of  Reynolds  number  are  shown  in  Figure  4.1.  Over  most 
of  the  Reynolds  number  range  tested,  the  drag  coefficient  demonstrated  a  decrease  in  value  consistent  with  a  drag 
dominated  by  skin  friction.  However,  as  the  Reynolds  number  approached  the  upper  end  of  the  range  tested,  the 
drag  coefficient  appeared  to  level  off.  This  would  indicate  that  separation  was  taking  place  and  that  the  drag 
coefficient  had  become  “locked”  at  that  value,  behavior  similar  to  that  of  the  drag  coefficient  of  a  circular  cylinder 


*  Information  from  unpublished  report  provided  to  author  by  LCC  personnel. 
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or  sphere  as  the  Reynolds  number  increases.  With  the  results  given  above  for  the  uncertainties  in  the  various 
parameters,  the  resulting  uncertainty  in  the  Reynolds  number  was  typically  on  the  order  of  5%  (caused  mostly  by  the 
systematic  uncertainty  in  the  temperature  measurement).  The  average  uncertainty  in  the  drag  coefficient  was  3.6%. 

The  evidence  for  separation  on  the  rear  of  the  model  is  shown  in  Figures  4.2  through  4.8.  These  figures 
show  the  results  of  the  LDV  surveys  in  the  vicinity  of  the  model  stem.  The  component  denoted  by  Ux  is  the 
streamwise  velocity  component.  The  component  denoted  by  Uz  is  the  transverse  velocity  component,  which  is  the 
component  perpendicular  to  the  LCC  test  section  floor  and  positive  up  ( i.e .  this  velocity  component  is  measured  in  a 
reference  frame  attached  to  the  LCC  test  section).  Both  mean  and  fluctuating  value  contours  are  shown  for  both 
velocity  components.  The  contour  that  forms  the  upper  boundary  of  each  region  is  the  outline  of  the  submarine  hull 
and  dummy  hub,  plus  the  centerline  of  the  water  tunnel.  The  measurements  were  made  on  the  lower  surface  of  the 
model,  to  avoid  any  effects  due  possibly  to  the  model  support  stmt,  hence  the  negative  values  of  z/R  reported  (R  is 
the  model  hull  radius).  Figures  4.2,  4.3,  and  4.4  are  summaries  of  the  measurements  at  tunnel  set  speeds  of  2  kts,  5 
kts,  and  10  kts,  respectively.  Figures  4.5  through  4.8  are  expanded  views  of  each  component  at  a  tunnel  set  speed  of 
20  kts.  These  figures  do  show  a  very  small  separated  region  just  behind  the  elliptic  tip  at  the  trailing  edge  of  the 
model.  There  are  negative  values  of  Ux,  indicating  reversed  flow  behind  the  elliptic  tip.  There  are  also  large  values 
of  the  fluctuating  components  along  the  edge  of  the  presumably  reversed  flow  region,  indicating  the  presence  of 
significant  mean  velocity  gradients  which  typically  bound  separated  and/or  reversed  flow  regions.  It  should  be 
noted  that  these  are  apparently  the  only  separated  flow  regions  to  be  found  for  this  model.  These  figures  show  the 
presence  of  the  separated  flow  over  the  entire  range  of  tunnel  speeds  tested.  It  is  conjectured  that  the  drag  of  the 
submarine  hull  behaves  similarly  to  other  geometries.  At  lower  Reynolds  number,  the  drag  is  dominated  by  the  skin 
friction,  which  decreases  as  the  Reynolds  number  increases.  Eventually  the  contribution  to  the  drag  from  the  skin 
friction  drops  below  that  from  pressure  drag  caused  by  the  separated  flow  at  the  stem  of  the  model.  Since  this 
pressure  drag  is  essentially  independent  of  Reynolds  number,  the  overall  drag  coefficient  levels  off. 

4.3.2  Propeller  thrust 

The  propeller  thrust  studies  were  conducted  at  LCC  tunnel  set  speeds  of  5  kts,  10  kts,  15  kts,  and  20  kts, 
and  for  propeller  speeds  ranging  from  200  RPM  to  1200  RPM.  The  results  for  the  propeller  thrust  as  a  function  of 
advance  ratio  are  shown  in  Figure  4.9.  The  rather  larger  uncertainty  bars  in  Figure  4.9a  are  the  result  of  the 
systematic  uncertainty  in  the  propeller  dynamometer.  In  Table  4.2,  the  systematic  uncertainty  is  8  lbf,  which  is  close 
to  the  mean  value  of  the  propeller  force  at  some  propeller  speeds  for  a  tunnel  set  speed  of  5  kts.  As  the  free  stream 
speed  increased,  the  propeller  force  increased  also,  and  the  relative  uncertainty  decreased.  Figure  4.9a  shows  all  of 
the  results,  and  Figure  4.9b  shows  the  results  for  10  kts,  15  kts,  and  20  kts  over  a  more  limited  range  of  advance 
ratio.  A  typical  value  of  the  uncertainty  in  the  thrust  coefficient,  neglecting  some  of  the  more  extreme  cases,  was 
approximately  8%.  The  uncertainty  in  the  advance  ratio  was  typically  1 .8%.  The  data  show  that  the  results  were  not 
dependent  on  Reynolds  number  over  the  range  tunnel  and  propeller  speeds  tested.  The  data  tended  to  fall  on  a  single 
curve  when  plotted  versus  advance  ratio.  This  demonstrated  the  repeatability  of  the  data,  since  similar  values  of 
advance  ratio  could  be  obtained  for  widely  different  values  of  tunnel  and  propeller  speed. 
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4.3.3  Boundary  layer  measurements 

In  addition  to  the  drag  and  thrust  measurements  that  were  the  main  reason  for  this  study,  the  boundary  layer 
velocity  at  several  axial  locations  toward  the  rear  of  the  model  was  measured  in  some  detail.  These  results  are 
discussed  here.  Although  these  measurements  were  not  directly  related  to  the  main  purpose  of  this  study,  it  was 
believed  that  they  might  be  of  value  to  researchers  attempting  to  validate  methods  related  to  boundary  layers  at  very 
large  values  of  the  Reynolds  number. 

The  locations  of  the  boundary  layer  surveys  are  shown  in  Figure  4.10.  Because  the  model  strut  was  located 
on  the  upper  surface  of  the  model,  the  boundary  layer  surveys  were  conducted  on  the  lower  side,  away  from  the 
strut.  The  surveys  were  conducted  in  the  model  center  plane  (the  plane  containing  the  model  axis  and  perpendicular 
to  the  test  section  floor).  In  Figure  4.10,  the  flow  is  from  right  to  left,  as  indicated  by  the  arrow,  and  the  abscissa  is 
distance  in  mm  measured  from  the  stem  of  the  model.  Because  of  the  difficulties  involved  in  attempting  to  measure 
the  boundary  layer  in  a  direction  perpendicular  to  the  model  wall,  the  surveys  were  conducted  in  a  direction 
perpendicular  to  the  submarine  model  axis,  as  shown  in  Figure  4.10.  The  x/L  location  of  each  traverse  is  noted  on 
the  figure.  Two  velocity  components  were  measured,  the  u  component  in  the  direction  parallel  to  the  model  axis, 
and  the  v  component  perpendicular  to  the  LCC  test  section  floor  (and  hence  perpendicular  to  the  model  axis).  In 
Figure  4.10,  u  is  positive  to  the  left  and  v  is  positive  up  (in  the  positive  z  direction).  Mean  and  fluctuating 
components  were  measured,  as  was  the  correlation  (u  V)l/2. 

The  free  stream  pressure  distribution  is  shown  in  Figure  4.11.  These  pressures  were  obtained  from  taps 
mounted  on  the  wall  of  the  LCC  test  section.  The  x/L  values  in  this  figure  are  identical  to  the  x/L  values  for  the 
model.  The  first  survey,  at  x/L  =  0.615,  was  conducted  in  a  region  of  essentially  constant  pressure.  The  other  three 
surveys  were  conducted  in  a  region  of  adverse  pressure  gradient. 

The  profiles  of  mean  and  fluctuating  velocity  components  and  of  the  Reynolds  stresses  are  shown  in 
Figures  4.12  through  4.15.  The  boundary  layer  thickness  6  was  defined  as  the  y  coordinate  where  the  u  velocity 
component  reached  a  value  of  0.9917,  where  U  was  the  reference  velocity  at  a  particular  location.  The  y  coordinate 
was  obtained  by  subtracting  the  local  model  radius  R,oca)  from  the  measured  z  coordinate  and  then  changing  the  sign 
(the  z  values  in  Figure  4.10  are  measured  from  the  model  centerline).  Each  caption  in  Figures  4.12  through  4.15 
gives  the  local  reference  velocity  U,  the  boundary  layer  thickness,  and  the  ratio  of  the  boundary  layer  thickness  to 
the  local  model  radius  for  each  boundary  layer  survey. 

Figure  4.12  shows  the  boundary  layer  at  x/L  =  0.615.  At  this  location,  the  wall  is  parallel  to  the  free  stream 
direction  and  so  the  v  velocity  component  is  essentially  zero,  as  one  might  expect.  In  this  case  the  boundary  layer 
thickness  as  just  defined  corresponds  to  the  standard  definition  of  the  boundary  layer  thickness.  The  boundary  layer 
appears  to  behave  as  a  typical  turbulent  boundary  layer.  As  x/L  increases,  however,  the  local  slope  of  the  model 
wall  becomes  significant  and  the  vertical  velocity  component  v  becomes  more  pronounced,  as  can  be  seen  most 
clearly  in  Figure  4.15.  Note  that  the  vertical  velocity  component  in  this  figure  is  positive.  As  the  flow  curves 
around  the  stem  of  the  model,  it  is  directed  upward,  and  thus  yields  a  positive  velocity  relative  to  the  fixed  reference 
frame  of  the  LCC  test  section.  Note  also  at  this  location  that  the  ratio  8 /Rhca,  has  reached  a  value  of  0.89,  using  the 
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definition  of  5  given  above.  In  this  region,  it  is  not  likely  that  the  boundary  layer  would  behave  as  a  two- 
dimensional  boundary  layer. 

Because  of  the  lack  of  measurements  normal  to  the  model  wall  and  the  additional  complications  of  having 
to  deal  with  an  axisymmetric  turbulent  boundary  layer,  no  further  analysis  of  the  boundary  layer  data  was  attempted. 
The  data  are  included  here  primarily  for  validation  of  computational  methods,  particularly  in  the  instance  in  which  a 
researcher  might  attempt  to  compute  the  flow  past  this  geometry  with  the  goal  of  capturing  the  boundary  layer 
accurately.  The  data  could  then  be  used  to  judge  the  success  of  the  method.  The  “raw  data”  for  the  velocity  profiles 
and  the  pressure  coefficient  distributions  are  included  in  tables  immediately  following  this  section. 


Table  4.3  Boundary  layer  data:  x/L  =  0.615,  [/=  5.42  m/s,  8  =  53.0  mm 
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5.  COMPARISON  OF  COMPUTATIONAL  AND  EXPERIMENTAL  RESULTS 


5.1  Drag 

The  results  of  the  computational  studies,  the  small-scale  experimental  studies,  and  the  LCC  large-scale 
experimental  studies  for  the  drag  coefficient  of  the  bare  submarine  hull  geometry  are  shown  in  Figure  5.1.  Figure 
5.1a  compares  all  of  the  results,  and  Figure  5.1b  concentrates  on  the  LCC  model  scale  results.  The  reasons  for  the 
large  uncertainty  bars  on  the  MSU  test  results  have  been  discussed  previously.  Figure  5.1b  is  essentially  a 
reproduction  of  Figure  (B)4.2.  However,  it  should  be  pointed  out  that  when  Figure  (B)4.2  was  prepared,  the 
pressure-gradient  correction  had  not  been  applied  to  the  LCC  drag  coefficient  data.  When  the  correction  is  added, 
the  computations  and  experiments  are  in  somewhat  better  agreement. 

The  particular  feature  to  note  in  Figure  5.1  is  the  difference  between  the  continual  decrease  in  the  drag 
coefficient  predicted  by  the  computations  and  the  leveling  off  of  the  drag  coefficient  at  a  Reynolds  number  around 
108.  If  the  conjecture  presented  in  Chapter  4  is  correct,  this  is  the  point  where  the  pressure  drag  caused  by  the 
separation  on  the  stem  of  the  model  becomes  comparable  to  the  skin  friction  drag  caused  by  the  viscous  stresses. 
There  was  no  indication  in  the  computational  results  that  separation  occurred  anywhere.  Figure  5.2  (which  is  Figure 
(A)4. 1 1  repeated)  shows  a  very  small  region  of  low- velocity  flow  at  the  stem  of  the  model  but  no  apparent  reversed 
flow.  Other  figures  in  Appendix  A  (including  those  listed  in  Table  (A)4.2  on  p.  46  of  Appendix  A)  also  indicate  no 
reversed  flow.  It  appears  from  these  results  that  perhaps  the  small  but  critical  separated  regions  were  not  captured  or 
not  computed  accurately  in  the  numerical  simulations,  and  so  the  contribution  to  the  total  drag  from  the  pressure 
drag  was  less  than  it  should  have  been.  Thus,  in  the  computational  simulations,  the  drag  behavior  appears  to  have 
been  dominated  over  the  entire  Reynolds  number  range  by  the  skin  friction  drag  arising  from  the  viscous  stresses.  It 
is  possible  that  either  the  grid  resolution  for  the  computational  simulations  was  not  sufficient  to  capture  the 
separated  region  or  that  the  absence  in  the  computational  solutions  of  a  separated  flow  region  with  reversed  flow  is 
due  to  a  problem  with  the  turbulence  model. 

An  additional  feature  to  note  in  Figure  5.1  is  that  despite  the  large  uncertainty  bars,  the  mean  values  of  the 
drag  coefficient  obtained  from  the  experiments  on  the  small-scale  model  are  reasonably  consistent  with  the 
computational  results.  This  agreement  helps  to  confirm  the  validity  of  the  wake  profile  method  obtained  in  this 
study  and  suggests  that  the  uncertainty  in  the  value  of  the  drag  coefficient  may  have  been  overestimated. 
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5.2  Propeller  thrust 

The  results  for  the  computational  and  experimental  studies  of  propeller  thrust  are  shown  in  Figures  5.3  and 
5.4.  Figure  5.3  shows  all  of  the  computational  and  experimental  results,  and  Figure  5.4  separates  the  results  for 
different  values  of  the  Reynolds  number.  As  was  pointed  out  earlier,  the  computational  results  are  in  very  good 
agreement  with  the  experimental  measurements.  This  is  a  tribute  to  the  capability  of  the  computational  scheme  to 
simulate  a  very  complex,  unsteady  flow  field.  The  computational  studies  were  not  extended  to  the  Reynolds  number 
corresponding  to  a  full-scale  prototype.  However,  the  results  for  thrust  coefficient  plotted  against  advance  ratio 
seem  to  be  relatively  independent  of  Reynolds  number,  suggesting  that  the  computations  would  be  equally  as 
accurate  at  larger  Reynolds  number.  The  one  caveat  to  keep  in  mind  in  considering  the  results  of  the  computational 
simulations  of  the  propelled  model  is  the  difficulty  experienced  by  one  of  the  turbulence  models  discussed  in 
Section  4.2  of  Appendix  A.  Both  the  “q-co”  and  Spalart-Allmaras  models  performed  well  for  the  isolated  propeller 
studies.  However,  when  the  propelled  model  study  was  attempted  using  the  “q-ro”  model,  the  model  became 
unstable.  This  instability  was  attributed  to  the  need  to  have  a  more  refined  grid  in  the  vicinity  of  the  sharp  propeller 
blade  edges. 

In  the  calculations  of  both  bare  hull  drag  and  propeller  thrust,  there  were  indications  that  further 
refinements  in  the  computational  methods  were  needed.  Either  a  higher  grid  resolution  or  a  better  (but  probably 
more  complex)  turbulence  model  was  required.  It  is  possible  that  both  were  required.  Grid  refinement  and 
turbulence  model  studies  would  probably  answer  this  question,  but  in  a  sense  the  lack  of  an  answer  is  actually  what 
the  current  study  was  looking  for.  The  idea  was  to  take  an  existing  computational  method,  use  it  “as  is”  over  a  large 
range  of  Reynolds  number,  and  then  examine  the  results  in  the  context  of  experimental  measurements.  The  fact  that 
has  to  be  kept  in  mind  when  considering  these  results  is  that  at  the  largest  Reynolds  number,  there  would  be  no 
experimental  results  with  which  to  compare.  What  we  can  leam  about  the  computational  method  at  the  lower  values 
of  Reynolds  number  has  to  be  carried  over  into  the  higher  Reynolds  number  range. 

As  just  explained,  it  is  possible  that  this  flow  required  a  computational  grid  with  a  higher  resolution  than 
the  grid  actually  used.  The  greater  grid  refinement  seemed  to  be  needed,  in  the  case  of  the  drag  studies,  to  capture 
the  separated  flow  to  produce  the  correct  drag  behavior  and,  in  the  case  of  the  propeller  thrust  studies,  to  capture 
accurately  the  propeller  edge  geometry  so  that  one  of  the  two  turbulence  models  would  work  properly.  The  need  for 
a  highly-refined  grid  on  a  (relatively)  large  geometry  is  a  common  problem  in  computational  fluid  dynamics, 
because  it  implies  the  need  for  a  very  large  number  of  grid  points,  increasing  the  computational  time  and  effort 
required  to  obtain  a  solution.  This  issue  would  have  to  be  monitored  in  attempting  to  simulate  flows  at  very  large 
Reynolds  number. 


25 


6.  CONCLUSIONS 


In  order  for  computational  solutions  to  simulate  accurately  flows  at  very  large  Reynolds  number,  the 
computational  grid  must  have  sufficient  resolution  to  capture  very  small  phenomena  in  the  flow  field,  because  these 
phenomena  can  still  have  a  significant  impact  on  the  overall  flow.  Care  must  also  be  taken  in  the  choice  and 
implementation  of  a  turbulence  model.  The  computational  methods  used  in  this  study  produced  results  for  the  drag 
that  did  not  agree  with  the  experimental  results,  and  this  was  attributed  to  the  inability  of  the  computational  method 
as  implemented  to  resolve  accurately  the  separated  flow  behind  the  submarine,  as  a  result  of  insufficient  grid 
resolution,  or  difficulties  with  the  turbulence  model,  or  both.  The  computational  method  did  an  excellent  job  of 
predicting  the  thrust  produced  by  the  propeller,  as  long  as  the  grid  resolution  did  not  interfere  with  the  turbulence 
model.  Because  the  large  Reynolds  number  is  usually  the  result  of  a  “large”  geometry,  the  computational  grids  will 
be  correspondingly  “large.”  The  large  grid  with  high  resolution  will  have  a  very  large  number  of  points,  requiring  a 
considerable  amount  of  computational  time  and  memory.  In  other  words,  little  things  count. 
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PROPELLER  4381 


Figure  1.3  Propeller  4381  geometry 


Figure  1.4  Open-' 
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Figure  2.1  Computational  results  for  bare  hull  (LCC  model  scale);  comparison  between  turbulence  models 


Figure  2.2  Computational  results  for  bare  hull  geometry  and  all  scales 
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Figure  3.1  Streamwise  variation  of  velocity  deficit  with  uncertainties  included 


Figure  3.2  Streamwise  variation  of  wake  width  parameter  with  uncertainties  included 


Mississippi  State  University 
Submarine  Near  Wake  &  Aft  Boundary  Layer 
Tunnel  Velocity  1 .05  m/s 
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Ulocal  =  1.05  m/s 
MSU02.tec 
07/12/2001 


Figure  4.2  Results  of  LD  V  surveys  of  velocity  flow  field  at  stern  of  model  at  a  tunnel  set  speed  of  2  kts  (upper 

bounding  contour  is  hull  outline  and  tunnel  centerline) 
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Mississippi  State  University 
Submarine  Near  Wake  &  Aft  Boundary  Layer 
Tunnel  Velocity  2.62  m/s 


x/L  =  1 

Umodel  =  2.62  m/s* 
U local  =  2.66  m/s 
MSU05.tec 
7/11/01 


Figure  4.3  Results  of  LDV  surveys  of  velocity  flow  field  at  stem  of  model  at  a  tunnel  set  speed  of  5  kts  (upper 

bounding  contour  is  hull  outline  and  tunnel  centerline) 


Mississippi  State  University 
Submarine  Near  Wake  &  Aft  Boundary  Layer 
Tunnel  Velocity  10.32  m/s 


x/L  =  1 

Umodel  =  10.32  m/s* 
Ulocal  =  10.51  m/s 
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Figure  4.4  Results  of  LDV  surveys  of  velocity  flow  field  at  stem  of  model  at  a  tunnel  set  speed  of  20  kts 
(upper  bounding  contour  is  hull  outline  and  tunnel  centerline) 
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Mississippi  State  University 
Submarine  Near  Wake  &  Aft  Boundary  Layer 
Tunnel  Velocity  10.32  m/s 


x/L  =  1 

Umodel  =  10.32  m/s* 
Ulocal  =  10.51  m/s 
MSU20.tec 
07/12/2001 


Figure  4.5  Detail  of  mean  streamwise  velocity  component  at  stem  of  model  at  a  tunnel  set  speed  of  20  kts 
(upper  bounding  contour  is  hull  outline  and  tunnel  centerline) 
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Mississippi  State  University 
Submarine  Near  Wake  &  Aft  Boundary  Layer 
Tunnel  Velocity  10.32  m/s 
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x/L  =  1 

Umodel  =  10.32  m/s* 
Ulocal  =  10.51  m/s 
MSU20.tec 
07/12/2001 


Figure  4.6  Detail  of  fluctuating  streamwise  velocity  component  at  stern  of  model  at  a  tunnel  set  speed  of  20 
kts  (upper  bounding  contour  is  hull  outline  and  tunnel  centerline) 
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Mississippi  State  University 
Submarine  Near  Wake  &  Aft  Boundary  Layer 
Tunnel  Velocity  10.32  m/s 
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Figure  4.7  Detail  of  mean  transverse  velocity  component  at  stem  of  model  at  a  tunnel  set  speed  of  20  kts 
(upper  bounding  contour  is  hull  outline  and  tunnel  centerline) 
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Mississippi  State  University 
Submarine  Near  Wake  &  Aft  Boundary  Layer 
Tunnel  Velocity  10.32  m/s 
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Figure  4.8  Detail  of  fluctuating  transverse  velocity  component  at  stem  of  model  at  a  tunnel  set  speed  of  20  kts 
(upper  bounding  contour  is  hull  outline  and  tunnel  centerline) 
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(a)  Entire  data  set 


(b)  Results  with  data  acquired  at  5  kts  removed  and  range  in  J  reduced 
Figure  4.9  Propeller  thrust  coefficient  variations  with  advance  ratio 


flow  direction 


L  -  x  (mm) 

Figure  4.10  Locations  of  boundary  layer  surveys  ( L  =  model  length;  numbers  indicate  x/L  positions  of 
boundary  layer  traverses) 


Figure  4.11  LCC  free  stream  pressure  distribution  obtained  from  test  section  wall  taps 
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(a)  Mean  velocity  components 


(b)  Fluctuating  components 


(c )  Reynolds  stresses 

Figure  4.12  Boundary  layer  profiles  for  x/L  =  0.615  (t/=  5.42  m/s,  5=  53.0  mm,  8/Ri„cai=  0.17) 
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(b)  Comparison  of  LCC  model  scale  results 
Figure  5.1  Comparison  of  computational  and  experimental  results  for  drag 
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fl  U/TJref  1 

Figure  5.2  Computed  xial  velocity  contours  for  ReL  =  32  million 


(b)  Expanded  scale  in  advance  ratio 

Figure  5.3  Comparison  of  experimental  and  computation  results  for  propeller  thrust  coefficient  (all  values  of 

Reynolds  number) 


(a)  &?£.«■/,=  18. 1(106),  ReLcomp  =  16(106) 
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(b)  ReL,  exp  =  35.6(106),  ReL,  cnmp  =  32(106) 

Figure  5.4  Comparison  of  experimental  and  computational  results  for  propeller  thrust  at  different  values  of  the 

Reynolds  number  (continued) 
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CHAPTER  I 


INTRODUCTION 

Today’s  submarines  travel  at  speeds  in  excess  of  25  knots  and  practice  extreme 
maneuvers  that  push  the  very  limit  of  their  design  and  structural  integrity.  It  is 
important  to  understand  the  forces  encountered  when  such  speeds  and  maneuvers  are 
executed.  There  are  two  ways  to  obtain  this  information:  experimentation  and 
computation. 

One  of  the  difficulties  associated  with  experimentation  is  the  inability  to 
effectively  test  a  full-scale  submarine.  Experiments  must  be  carried  out  on  a  smaller 
scale  and  Reynolds  scaling  used  to  determine  full-scale  forces.  The  largest  scale 
model  is  a  quarter-scale  model,  but  the  data  from  those  tests  are  classified.  The  largest 
model  scale  data  available  for  this  study  is  from  a  6%  scale  model  tested  at  the  Large 
Cavitation  Channel  (LCC)  in  Memphis,  Tennessee.  Reynolds  scaling  does  not 
account  for  all  the  differences  in  the  flow,  and  with  such  a  large  jump  between  the 
scales,  small  uncertainties  can  create  large  inaccuracies  on  the  full-scale  level. 

Another  method  to  determine  the  forces  on  a  submarine  at  extreme  conditions 
is  to  computationally  simulate  the  maneuvers.  Computationally  simulating  such 
conditions  is  a  relatively  new  avenue  of  pursuit  and  pushes  current  numerical 
algorithms  to  their  limits. 
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1.1  Previous  Studies 


Studies  have  been  performed  on  both  the  experimental  side  and  the 
computational  side  for  fully  submerged  bodies  or  submarines.  A  recent  experimental 
venture  used  to  validate  many  CFD  methods  is  the  Defense  Advanced  Research 
Projects  Agency  (DARPA)  funded  study  of  the  SUBOFF  model  [1],  Experiments  on 
the  SUBOFF  model  were  performed  in  the  David  Taylor  Model  Basin  in  1990.  The 
SUBOFF  hull  is  14.3  feet  long  and  was  tested  at  a  Reynolds  number  of  1.4  x  107. 
Force  and  moment  data,  as  well  as  pressure  and  skin  friction  coefficient  distributions, 
were  recorded.  The  model  was  tested  in  several  configurations,  but  all  were 
unpropelled. 

The  computational  investigation  of  complex  hydrodynamic  configurations  has 
evolved  from  CFD  studies  involving  aerodynamic  bodies  that  began  in  the  1980s. 
These  studies  have  focused  on  fully  submerged  prolate  spheroids  and  SUBOFF 
configurations  because  of  the  availability  of  the  experimental  data  for  validation. 

Such  studies  have  included  force  and  moment  validation  [2],  and  full-scale  studies  [3] 
Little  has  been  done  with  propelled  simulations,  and  even  less  with  the  validation  of 
those  simulations.  Only  within  the  past  five  years  have  propelled  studies  with  the 
SUBOFF  model  been  investigated.  Some  work  has  even  involved  six-degree-of- 
freedom  simulations  [4],  but  these  are  still  in  the  development  stages.  This  study  is 
one  of  the  first  to  validate  fully  submerged  and  propelled  body  simulations  with 


experimental  data. 


1.2  Present  Study 
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This  study  will  focus  on  several  aspects  of  the  problem  described  above.  The 
propulsor  and  body  interaction  will  be  investigated,  as  well  as  the  differences  in 
computed  flow  fields  as  Reynolds  number  varies.  Turbulence  model  effects  and  the 
integrity  of  the  flow  solver  to  run  long  duration  unsteady  problems  will  also  be  noted. 

There  are  four  major  parts  to  this  study:  1)  simulation  of  the  flow  past  an 
isolated  marine  propeller  with  0°  skew,  2)  simulation  of  the  flow  past  a  unappended 
submarine  at  three  different  length  scales  with  large  corresponding  differences  in 
Reynolds  number,  3)  simulation  of  the  flow  past  a  forward  propelled  submarine  with 
asymmetric  stem  appendages  at  various  advance  coefficients,  and  4)  simulation  of  the 
flow  past  a  submarine  with  asymmetric  stem  appendages  in  crashback.  Reynolds 
numbers  for  this  study  are  based  on  the  length  of  the  hull.  All  simulations  are 
validated  by  experimental  data  where  available. 

The  bare  hull  and  isolated  propeller  studies  provide  the  groundwork  to  analyze 
the  propulsor  and  body  interaction.  The  forward  propelled  and  appended  submarine 
cases  will  be  compared  against  the  component  cases  in  order  to  analyze  what  effect 
the  impingement  of  hull  boundary  layer  and  appendage  wakes  on  the  propeller  flow 
field  have  on  the  overall  flow  field.  This  relationship  is  important  because  it  can 
provide  information  to  better  design  a  quieter  submarine. 


1.3  Chapter  Overview 

The  next  chapter  will  discuss  the  experiments  and  the  flow  solver  that  are 
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pertinent  to  the  sufficiency  of  this  study.  Chapter  III  discusses  the  methodology 
behind  the  four  major  components  of  this  study.  Chapter  IV  presents  the  results  of  the 
study  and  discusses  the  relevancy  of  each.  The  last  chapter  discusses  the  study  as  a 
whole.  Concluding  remarks  and  recommendations  for  future  work  are  made. 


CHAPTER  II 


BACKGROUND 

2.1  Description  of  the  Experiments 

2.1.1  Isolated  Propeller 

The  isolated  propeller  tests  were  carried  out  in  the  Naval  Ship  Research  and 
Development  Center’s  (NSRDC)  deep-water  basin  in  1971  [5].  Open  water 
propulsion  tests  were  conducted  on  four  marine  propellers  with  varying  degrees  of 
skew.  The  propeller  of  interest  to  the  current  study  is  the  P438 1 .  The  P438 1 
propeller  has  five  blades  with  0°  skew. 

Forward  and  backing  propulsion  tests  were  performed.  The  forward  tests  are 
of  primary  interest  to  this  study.  A  propeller  boat  was  instrumented  with  a  gravity 
dynamometer  for  the  forward  tests.  All  forward  tests  were  run  at  a  rotational  speed  of 
7.8  rotations  per  second  (rps)  with  the  speed  of  advance  varying  from  2.75  fps  to 
10.14  fps.  Thrust  and  torque  data  were  recorded. 

2.1.2  Submarine  Models 

The  6%  scale  model  submarine  experiments  were  conducted  at  the  Large 
Cavitation  Channel  (LCC)  in  Memphis,  Tennessee.  Three  sets  of  experiments  were 
performed:  bare  hull,  forward  propelled,  and  crashback.  All  of  these  experiments 
used  an  existing  axisymmetric  hull  (referred  to  as  the  LCC  hull),  which  totaled  22.964 
feet  in  length. 
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Figure  2.1  shows  the  layout  of  the  model  in  the  test  chamber.  The  view  of 
Figure  2.1  is  from  the  side.  The  submarine  model  was  suspended  from  the  ceiling  of 
the  test  chamber  by  a  strut. 


For  the  propelled  cases,  the  LCC  hull  was  outfitted  with  a  one-foot  diameter 
P4381  propeller  and  four  asymmetric  stem  appendages  set  in  a  conventional  cruciform 
pattern.  The  bare  hull  tests  were  conducted  on  the  LCC  hull  without  the  propeller  and 
stem  appendages. 

Although  the  P4381  is  a  marine  propeller  designed  for  surface  ships,  it  was 
chosen  for  these  experiments  because  it  was  readily  available,  unclassified  and 
previous  data  were  available.  The  P4381  still  maintains  many  of  the  characteristics 
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associated  with  a  submarine  propeller.  Importantly,  it  is  designed  to  operate  within 
the  boundary  layer  and  wake  of  a  much  larger  body. 

For  this  particular  hull,  the  one-foot  P4381  is  slightly  undersized  and  mounted 
further  downstream  from  the  stem  appendages  than  is  typical.  This  does  not 
significantly  affect  the  experiments,  but  must  be  considered  in  analysis. 

For  the  bare  hull,  61  runs  were  made  at  various  Reynolds  numbers  (based  on 
body  length)  ranging  from  4.5  x  106  to  1 .23  x  108.  A  Laser  Doppler  Velocimeter 
(LDV)  was  used  to  measure  the  tunnel  velocities  with  the  reference  velocity  taken  at 
x/L  of  0.866,  where  x  is  the  axial  position  and  L  is  the  length  of  the  body.  Force  and 
moment  data  were  also  recorded. 

The  forward  propelled  cases  were  run  at  various  values  of  the  advance 
coefficient  (J,  defined  as  J  =  V/nD  where  V  is  the  forward  speed,  n  is  the  propeller 
rotational  speed  and  D  is  the  propeller  diameter)  for  four  different  tunnel  speeds:  5, 
10,  15,  and  20  knots.  These  speeds  yield  Reynolds  numbers  (based  on  body  length) 
of  1.6  x  107,  3.2  x  107,  4.8  x  107,  and  6.4  x  107.  Thrust  and  drag  measurements  were 
taken  for  each  case.  A  total  of  35  runs  were  made  for  the  forward  propelled 
submarine.  Fourteen  of  those  cases  are  investigated  in  this  study. 

Thirty-one  runs  were  made  for  the  submarine  in  crashback.  The  crashback 
simulation  involves  forward  tunnel  speed,  but  with  the  submarine  propeller  rotating  in 
reverse.  These  cases  were  ran  at  various  advance  coefficients  for  four  different  tunnel 
speeds:  2.5,  5,  7.5,  10  knots.  These  speeds  yield  Reynolds  numbers  of  8  x  106,  1.6  x 
107,  2.4  x  107,  and  3.2  x  107.  Force  and  moment  data  were  recorded  for  the  crashback 
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experiments.  Two  cases  at  a  Reynolds  number  of  1 .6  x  107  and  three  cases  at  a 
Reynolds  number  of  3.2  x  107  are  investigated  in  this  study. 

Additional  bare  hull  experiments  are  currently  being  carried  out  on  a  LCC  hull 
model  that  is  1 1.278  inches  in  length  and  is  identical  in  hull  contour  to  the  model  used 
in  the  LCC  tests.  These  tests  are  performed  in  the  water  tunnel  located  in  Patterson 
Engineering  Laboratories  at  Mississippi  State  University.  The  Reynolds  numbers 
from  these  tests  range  from  7  x  105  to  2.5  x  106.  The  purpose  of  these  experiments  is 
to  provide  an  additional  Reynolds  number  range  for  comparison  purposes. 

The  previous  descriptions  are  intended  to  provide  context  for  the 
computational  simulations.  A  complete  description  of  the  experiments  can  be  found  in 
[6], 

2.2  Description  of  the  Flow  Solver 

The  flow  solver  used  in  this  study  is  U2NCLE  (Unstructured  Unsteady 
Computation  of  Fie/d  Equations).  U2NCLE  is  a  Reynolds-averaged  Navier-Stokes 
(RANS)  incompressible  flow  solver  capable  of  performing  viscous,  high  Reynolds 
number  flow  simulations  using  unstructured  grids.  In  [7]  the  U2NCLE  solver  is 
demonstrated  for  large  scale  meshes,  and  is  shown  to  be  an  effective  solver  for 
complex  hydrodynamic  applications. 

In  this  section,  the  fundamentals  behind  the  flow  solver  are  explained.  If 
greater  detail  is  desired,  [8]  explains  the  finer  points. 


2.2.1  Governing  Equations 
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The  pseudo-compressibility  method  is  used  to  put  the  governing  equations  into 
a  time-dependent  form.  This  transforms  the  equations  into  a  hyperbolic  system  and 
allows  for  the  use  of  compressible  flow  algorithms  [9].  The  integral  form  of  the 
Navier-Stokes  equations  is  the  basis  for  finite  volume  flow  solvers  [10].  Written 
below  is  the  conservative  integral  form  of  the  RANS  equations 


\fm=  Jc 


V^-  +  \FM=  |  G-hdl 
dt  -  an 


(2.1) 


an 


where  «  is  the  outward  unit  normal  to  the  control  volume  V. 

Given  below  are  the  vector  of  dependent  variables  and  the  components  of  the 
inviscid  and  viscous  flux  vectors  with  the  pseudo-compressibility  factor  (P). 


(2.2) 


/B 

u 0  +  nxP 
V©  +  hyP 
w©  +  hzP 


G  ■  ft  = 


hxt„  +  nyrxy  +  nzrxz 

^ x Tvx  ^ xy  ”z*yz 

^ x  ^ xx  Hy^zy  ^z^zz 


(2.3) 


(2.4) 


The  variables  in  the  above  equations  are  normalized  by  the  characteristic 
length  and  free  stream  values  of  velocity,  density,  and  viscosity.  The  dependent 
variables  u,  v,  and  w  are  the  velocity  components  in  the  x,  y,  and  z  directions 
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respectively.  P  is  the  pressure.  The  unit  normal,  n,  is  broken  into  its  x,  y,  and  z 
components  for  the  inviscid  and  viscous  flux  vectors.  The  term  /?  in  the  inviscid  flux 
vector  is  the  pseudo-compressibility  factor  and  is  equal  to  15  for  this  study.  ©  is  the 
velocity  normal  to  the  control  volume  face.  The  viscous  stresses  from  Equation  2.4 
are  given  as 


%  =(.v  + 1‘, ) 


1  (dui  Suj 


Re  I  Sty  3x, 


(2.5) 


where  fi  is  the  molecular  viscosity  and  /u,  is  the  eddy  viscosity.  The  Reynolds  number 
is  denoted  Re  and  is  defined  as  LV/v,  where  L  is  the  length  of  the  body,  V  is  the 


forward  velocity  and  v  is  the  kinematic  viscosity  of  water. 

2.2.2  Numerical  Method 


The  basic  flow  solver  is  a  node-centered,  finite  volume,  implicit  scheme 
applied  to  general  unstructured  grids  with  non-simple  elements.  The  solver  is  based 
on  domain  decomposition  for  concurrent  solution  within  subdomains  assigned  to 
multiple  processors.  This  parallelization  allows  for  quick  turnaround  of  a  solution. 

The  basic  numerical  method  involves  reconstruction  of  the  solution  states, 
evaluation  of  the  residuals,  and  time  evolution  of  the  solution  in  each  control  volume 
[11].  Reconstruction  is  done  by  extrapolating  the  solution  at  the  vertices  to  the  faces 
of  the  surrounding  control  volume  using  a  higher-order  spatial  method.  The  gradients 
at  the  vertices  are  computed  using  the  unweighted  least  squares  method  and  then  the 
variables  at  the  interface  are  computed  using  a  first-order  Taylor  series.  The 
governing  equations  are  discretized  using  a  finite  volume  technique.  The  surface 
integrals  from  Equation  2.1  are  approximated  by  a  quadrature  over  the  surface  of  the 
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control  volume  of  interest.  The  evaluation  of  the  discrete  residual  is  performed 
separately  for  the  inviscid  and  viscous  terms.  After  the  spatial  terms  have  been 
discretized,  the  time  derivative  is  approximated.  For  the  time  evolution,  a  Newton 
iterative  scheme  is  used.  It  requires  the  solution  of  a  sparse  linear  system  at  each 
nonlinear  subiteration.  A  bi-directional  Gauss-Seidel  algorithm  is  used  to  solve  the 
system. 

2.2.3  Boundary  Conditions 

All  the  boundary  conditions  are  dealt  with  in  an  implicit  manner.  The  farfield 
conditions  are  managed  by  a  characteristic  variable  reconstruction.  Viscous 
conditions  are  enforced  by  modifying  the  linear  system  such  that  no  change  is  allowed 
in  the  velocity,  and  the  pressure  is  driven  according  to  the  imbalance  in  the  continuity 
equation  in  the  boundary  control  volume. 

A  symmetry  plane  boundary  condition  is  dealt  with  by  creating  a  layer  of 
phantom  cells  that  is  a  mirror  image  of  the  cells  inside  and  connected  to  the  symmetry 
plane.  The  control  volumes  on  the  symmetry  plane  are  closed  and  behave  just  as 
interior  control  volumes. 

2.2.4  Turbulence  Models 

There  are  two  turbulence  models  available:  1)  the  one-equation  Spalart- 
Allmaras  model  [12],  and  2)  the  two-equation  q-co  model  [13].  The  Spalart-Allmaras 
model  formulates  a  transport  equation  for  the  turbulent  Reynolds  number,  which  is 
then  related  to  the  turbulent  viscosity.  The  q-co  model  uses  a  transport  equation  each 
for  the  velocity  scale  and  length  scale  to  specify  the  distribution  of  the  eddy  viscosity. 


Although  the  two  equation  q-oo  model  is  very  robust  and  more  accurate,  it  is  also 
more  expensive  to  run  due  to  the  number  of  equations  to  be  solved.  The  q-ro  model 
also  very  sensitive  to  the  surface  grid,  while  the  Spalart-Allmaras  model  is  more 
forgiving.  For  this  study,  both  turbulence  models  are  employed. 


CHAPTER  III 


DESCRIPTION  OF  COMPUTATIONAL  CASE  PREPARATION 

The  flow  around  a  powered  submarine  is  very  complex  due  to  the  interactions 
among  the  hull  boundary  layer,  appendage  boundary  layers  and  wakes,  and  the 
propulsor.  This  study  breaks  the  problem  down  into  bare  hull,  isolated  propeller  and 
appended  and  powered  simulations.  The  flow  fields  about  the  bare  hull  and  isolated 
propeller  will  be  analyzed  to  determine  which  flow  characteristics  are  unique  to  the 
components.  The  flow  field  about  the  powered  and  appended  submarine  will  be 
compared  against  the  component  flow  fields  to  determine  what  effect  the  impingement 
of  the  hull/appendage  flow  on  the  propeller  flow  has  on  the  overall  flow  field. 

Another  aspect  of  this  study  is  the  comparison  of  the  hull  (bare  as  well  as 
appended)  flow  field  at  three  different  length  scales  to  examine  differences  in  flow 
characteristics.  It  has  been  observed  in  previous  studies  that  significant  differences  in 
the  submarine  flow  fields  exist  between  model  and  full  scale  flows  [14]  [15].  It  is  also 
noted  that  the  differences  are  widespread  and  not  localized  to  a  specific  location.  This 
leads  to  the  need  to  fully  simulate  the  full  scale  Reynolds  number  flow  in  order  to 
investigate  the  flow  in  its  entirety. 

Several  simulations  are  computed  to  provide  the  groundwork  for  the  above 
analyses.  The  LCC  model  and  isolated  propeller  simulations  are  corroborated  with 
available  experimental  data,  thus  serving  to  further  validate  the  flow  solver. 
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The  simulation  of  a  submarine  in  crashback  serves  as  the  ultimate  test  case 
for  this  study.  It  demonstrates  the  integrity  of  the  code  to  run  long  duration  unsteady 
simulations  yielding  results  in  good  agreement  with  experimental  data.  It  also 
underscores  the  level  of  complexity  the  interaction  of  the  component  flows  can  attain. 
It  is  hoped  that  the  information  gleaned  from  the  crashback  simulations  will  yield  new 
insight  into  the  nature  of  the  flow.  It  should  be  noted  that  for  this  study,  all 
simulations  are  computed  at  zero  angle  of  attack. 

The  table  in  Appendix  A  breaks  down  all  of  the  cases  and  gives  the  conditions 
for  each.  In  the  remainder  of  this  report,  the  cases  will  be  referred  to  by  the 
designation  given  to  them  in  Appendix  A.  Multiple  cases  are  referred  to  by  their  root 
case  name,  followed  by  an  x  that  represents  the  alphanumeric  that  designates  a 
specific  case,  i.e.  case. lx  stands  for  case. la,  case. lb,  case.lc,  etc. 

3.1  Bare  Hull 

The  bare  hull  simulations  provide  a  baseline  flow  field  for  the  propulsor  and 
hull  interaction  study.  These  cases  are  also  computed  for  three  different  length  scales 
for  the  scaling  analysis.  The  three  length  scales  correspond  to  the  1 1 .278-inch  small 
model  tested  in  Patterson  labs,  the  22.9614-foot  large  model  tested  in  the  LCC,  and  a 
383.5-foot  full-scale  submarine  [16].  Figure  3.1  shows  the  bare  hull  geometry 


configuration. 


Figure  3.1:  Bare  LCC  hull. 


The  first  step  in  performing  these  simulations  is  to  build  an  accurate 
computational  model  for  the  hull.  Geometry  preparation  and  surface  grid  generation 
are  performed  using  SolidMesh  [17]  with  AFLR  (Advancing  Front/Local 
Reconnection)  surface  grid  generation  [18].  The  computational  grids  are  multielement 
unstructured  meshes  generated  with  an  advancing  normal  methodology  for  the 
boundary  layer  elements,  and  an  AFLR  methodology  for  the  isotropic  volume 
elements  as  given  in  [18].  For  clarity  of  discussion,  the  axis  orientation  for  the 
submarine  is  as  follows:  x-axis  is  from  nose  to  tail,  y-axis  is  through  the  dorsal  (top) 
of  the  hull  and  the  z-axis  is  through  the  port  side  (left  looking  forward). 

A  composite  unstructured/structured  surface  grid  is  generated  for  the  LCC  hull 
and  a  multielement  unstructured  mesh  is  generated  for  the  volume.  A  composite  grid 
is  often  used  in  cases  where  there  is  a  large  axisymmetric  section  that  does  not  require 
high  resolution.  Usually,  with  a  reasonable  aspect  ratio  (less  than  5)  for  the  structured 
elements,  the  point  spacing  will  be  broader  than  with  all  unstructured  elements. 
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Figure  3.2  shows  the  surface  grid  with  the  structured  grid  elements  making 
up  most  of  the  hull.  The  sections  in  red  represent  the  unstructured  portions,  and  the 
section  in  blue  represents  the  structured  portion  of  the  surface  grid. 


Figure  3.2:  Composite  unstructured/structured  surface  grid. 

The  symmetry  of  the  body  is  exploited  to  save  points,  and  consequently, 
computational  time.  A  symmetry  plane  is  used  to  solve  for  only  half  of  the  hull.  The 
force  and  moment  data  are  multiplied  by  a  factor  of  two  to  take  this  into  consideration. 

One  grid  is  used  for  all  Reynolds  numbers  for  each  of  the  three  length  scales. 
The  initial  absolute  off- wall  spacing  of  the  volume  grid  from  the  body  is  5  x  10"7, 
which  yields  a  y+  distribution  of  approximately  1  over  the  entire  hull  for  the  high 
Reynolds  number  cases.  This  indicates  good  viscous  sublayer  resolution.  The  volume 
grid  around  the  stem  of  the  sub  body  is  shown  in  Figure  3.3. 


Figure  3.3:  Volume  element  distribution  about  the  LCC  stern. 


For  the  Reynolds  scaling  analysis,  the  axial  force  coefficient  or  drag 
coefficient  is  of  major  concern.  For  the  bare.2x  cases,  the  drag  coefficient  is  validated 
with  experimental  data  [Bridges],  The  sensitivity  of  axial  force  (drag)  to  the 
streamwise  grid  density  is  explored  by  building  two  additional  grids:  a  fully 
structured  grid  and  a  fully  unstructured  grid.  These  grids,  along  with  the 
structured/unstructured  composite  grid,  provide  three  different  axial  densities  and 
topologies  to  examine  the  sensitivity  of  axial  drag.  Table  3.1  gives  a  breakdown  of  the 
grid  statistics.  It  is  noted  that  the  fully  unstructured  and  the  unstructured/structured 
composite  grid  have  about  the  same  number  of  points. 


Table  3.1:  Grid  statistics. 
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Grid  Type 

Nodes 

Tetrahedra 

Pyramids 

Prisms 

Hexahedra 

Unstructured 

397302 

410857 

419 

617760 

0 

U  nstructured/Structured 

410715 

363008 

434 

662984 

0 

Structured 

1181037 

0 

120 

8160 

1143780 

All  bare  hull  simulations  were  run  using  local  time  stepping  and  converged  to 
obtain  steady-state  solutions.  Runs  are  made  with  both  turbulence  models.  Specific 
run  details  are  outlined  by  case  in  Appendix  A. 

The  bare  hull  cases  were  typically  run  for  1500  time  steps.  As  seen  in  the 
figures  below,  convergence  was  apparent  by  500  time  steps,  but  the  solution  was  run 
further  for  insurance.  All  simulations  were  carried  out  using  Lakota,  a  super  cluster  of 
256  dual  gigahertz  processor  Pentium  III  computers.  The  average  computation  time 
for  the  bare  hull  computations  was  approximately  10  hours. 
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Figure  3.4:  Convergence  history  by  scale  for  q-co  turbulence  model. 


Figure  3.4  shows  the  q- vector  convergence  histories  that  were  obtained  for  the 
various  Reynolds  numbers  using  the  q-co  turbulence  model.  As  Reynolds  number 
increases,  more  time  is  required  for  convergence. 
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Figure  3.5:  Convergence  history  by  grid  type  and  turbulence  model. 


It  is  noted  from  Figure  3.5  that  all  of  the  grids  converge  at  about  the  same  rate. 
The  noise  in  each  of  the  residuals  is  due  to  the  use  of  a  mean  flow  limiter  to  help  with 
solver  stability.  It  is  also  noted  that  there  is  no  significant  difference  in  the 
convergence  between  turbulence  models.  Although  convergence  appears  to  be 
relatively  better  with  the  structured  grid,  it  requires  more  computational  resources. 

The  typical  load  per  processor  is  approximately  100,000  grid  points.  The  structured 
grid  requires  1 2  processors  as  opposed  to  the  four  required  by  the 
unstructured/structured  and  unstructured  grids. 
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An  additional  case,  denoted  apphull.l,  was  run  to  isolate  the  effects  of  the 
stem  appendages  on  the  flow  field.  Figure  3.6  shows  the  volume  grid  in  the  vicinity 
of  the  stem  appendages.  From  this  picture,  the  asymmetry  of  the  mdders  can  be  seen. 
Boundary  layer  clustering  is  also  clearly  visible. 


Figure  3.6:  Volume  element  distribution  in  the  vicinity  of  stem  appendages. 

3.2  Isolated  Propeller 

Isolated  propeller  studies  are  performed  using  the  P4381  marine  propulsor. 

The  five  bladed  propeller  is  mounted  on  a  generic  hub.  Figure  3.7  shows  the  geometry 
used  for  the  isolated  simulations.  An  unstructured,  multielement  mesh  is  generated  for 
the  propeller  geometry.  This  mesh  consists  of  2.4  million  nodes,  16,000  pyramids,  2.1 
million  prisms,  and  7.9  million  tetrahedra. 


Figure  3.7:  Geometry  and  surface  grid  of  the  isolated  P4381. 


The  isoprop  1.x  and  isoprop2.x  cases  are  used  to  validate  the  simulation  and 
compare  the  turbulence  models.  The  isoprop3.x  cases  are  run  at  a  Reynolds  number 
of  32  million  for  the  advance  coefficients  that  correspond  to  the  10-knot  forward 
propelled  cases  from  the  LCC  experiments.  The  thrust  coefficients  and  velocity 
profiles  from  these  cases  will  be  compared  to  the  fwdsub.xx  cases. 

The  reference  velocity  used  for  non-dimensionalization  of  variables  for  flow 
solver  input  is  typically  the  free  stream  velocity,  U„.  For  the  isolated  propeller  cases, 
the  typical  non-dimensionalization  resulted  in  very  large  tip  speeds  with  respect  to  the 
forward  speed.  This  created  stability  problems  with  the  solver.  To  make  the  input 
values  more  manageable  for  the  solver,  the  rotational  tip  speed  was  used  as  the 
reference  velocity. 
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The  isolated  propeller  cases  were  run  unsteady  with  local  time  stepping  for 
approximately  300-400  time  steps.  The  solver  was  then  switched  to  minimum  time 
stepping  that  corresponded  to  1.5°  of  rotation  per  time  step  and  run  for  an  additional 
600-700  time  steps.  As  seen  in  the  given  force  history  (Figure  3.8),  around  time  step 
700  the  solution  has  reached  steady  state.  The  average  computational  time  required 
was  approximately  23  hours. 


o  100  200  300  400  500  600  700  800  900  1000 

Time  Step 


[  Ahead  J=0.95 - Ahead  J=0.97,  Re  =  32  million] 

Figure  3.8:  Force  history  for  isolated  propeller  thrust  coefficient.  J  =  0.95  correlates 
to  experimental  data;  J  =  0.97  correlates  to  hull/propulsor  study. 
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3.3  Forward  Propelled  Appended  Submarine 

The  forward  propelled  simulations  model  the  composite  flow  field  generated 
by  an  appended  and  propelled  submarine  in  straight  and  level  forward  motion.  The 
model  for  these  simulations  is  the  LC'C  hull  with  four  asymmetric  stern  appendages 
and  a  one-foot  P4381  propeller.  Figure  3.9  depicts  the  geometry  just  described. 


Figure  3.9:  Appended  LCC  Hull  with  P4381  Propeller 


An  unstructured  grid  is  generated  to  model  the  geometry  used  in  the  LCC 
experiments  minus  the  support  stmt.  The  grid  is  approximately  3.1  million  points  and 
yields  a  viscous  sublayer  resolution  for  all  four  Reynolds  numbers  of  y‘  of 
approximately  I .  The  unstructured  mesh  consists  of  19500  pyramids,  4.3  million 
prisms,  and  5  million  tetrahedra.  The  volume  grid  in  the  vicinity  of  the  stern 


appendages  is  shown  in  Figure  3.10.  This  figure  shows  the  point  clustering  in  the 
boundary  layer  region. 


Figure  3.10:  Volume  grid  in  the  vicinity  of  the  stem  appendages. 

Typical  nondimensionalization  by  free  stream  quantities  is  used  for  flow  solver 
inputs.  The  propelled  submarine  cases  are  run  with  local  time  stepping  for 
approximately  the  first  400  time  steps.  Time  stepping  is  then  switched  to  minimum 
with  step  size  set  to  1 .5°  rotation  per  time  step.  The  rotating  propeller  is  handled  by 
locally  regenerating  the  grid  in  the  vicinity  of  the  propeller  for  each  time  step.  This 
simulates  an  actual  rotating  propeller  on  a  fixed  body.  Three  Newton  iterations  are 
used  to  ensure  time  accurate  solutions.  These  cases  are  run  for  800-1000  time  steps 
where  a  periodic  state  has  been  well  attained.  The  force  history  given  in  Figure  3.10 
shows  convergence  to  a  periodic  state  for  cases  fwdsub.lc,  fwdsub.2d,  fwdsub.3c  and 
fwdsub.4c.  The  cases  are  all  for  an  advance  coefficient  of  J®1 .  Note  that  at  a  similar 
advance  coefficient,  the  solutions  converge  to  a  similar  thrust  coefficient. 
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S  Knots,  J*0.98 - to  Knots, J  =  0.98 - 15  Knots,  J  -  0  9S - 20  Knots.  J  -  1.17  | 


Figure  3.11:  Force  history  for  forward  propelled  submarine,  J~1 . 

The  inset  depicts  what  is  meant  by  periodic  state.  Due  to  the  unsteady  nature 
of  a  rotating  propeller,  the  thrust  coefficient  will  not  converge  to  a  constant  value  but 
will  converge  to  a  periodic  pattern  about  a  mean  value.  To  allow  for  this,  force 
coefficients  are  averaged  over  the  last  propeller  revolution. 

It  should  be  noted  that  due  to  the  orientation  of  the  computational  model  (x- 
axis  from  nose  to  tail),  thrust  is  shown  as  negative  and  drag  is  shown  as  positive. 
Typical  sign  designations  of  thrust  as  positive  and  drag  as  negative  are  applied  for  the 
results  section. 

Utilizing  3 1  processors,  the  computation  time  required  per  time  step  averaged 
1.25  minutes.  This  totaled  to  23-26  hours  of  total  run  time  for  a  converged  solution. 


3.4  Appended  Submarine  in  Crashback 
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The  crashback  maneuver  in  this  study  is  a  forward  moving  submarine  with  the 
propeller  rotating  in  reverse.  This  results  in  a  very  unsteady  flow  with  the 
establishment  of  an  asymmetric  ring  vortex  around  the  propeller  and  a  nonuniform 
pressure  distribution  on  the  propeller  blades.  The  complexity  and  the  unsteadiness  of 
this  flow  make  it  a  particular  challenge  to  model  numerically. 

Of  utmost  importance  is  the  accuracy  of  modeling  this  flow.  Crashback 
maneuvers  put  the  submarine  under  large  forces  and  high  stress,  and  complicate  the 
issue  of  controllability.  The  understanding  of  these  forces  and  stresses,  and  the 
physics  of  the  ring  vortex  can  greatly  aid  in  the  effective  design  and  safety  of  a 
submarine. 

Crashback  is  difficult  to  analyze  experimentally.  The  flow’s  unsteadiness 
makes  Fourier  analysis  difficult,  and  complicates  long  duration  testing.  It  is  also  hard 
to  simulate  the  deceleration  that  occurs  during  crashback.  If  computational 
simulation  can  yield  accurate  results,  it  simplifies  the  problem  of  data  acquisition  and 
generates  more  data. 

The  grid  used  for  the  forward  propelled  cases  is  used  for  the  crashback 
simulations.  The  major  difference  between  the  crashback  and  forward  propelled 
simulations  is  the  direction  of  propeller  rotation.  The  crashback  cases  also  require  a 
considerable  amount  of  time  to  converge  to  a  quasi-periodic  state,  if  they  ever 
converge.  These  cases  are  run  long  duration  for  1 8-20  propeller  revolutions.  Previous 
studies  have  only  run  8-10  propeller  revolutions,  so  this  pushes  the  limits  of  the  flow 
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solver.  Computation  time  per  time  step  averages  1 .25  minutes.  The  total  run  time 
for  a  crashback  simulation  is  approximately  1 20  hours  for  the  super  cluster  described 
previously.  Figure  3.12  shows  the  force  history  for  a  representative  crashback 
solution.  Notice  that  the  switch  from  local  to  minimum  time  stepping  occurs  around 
time  step  400  and  convergence  is  apparent  after  3000  time  steps. 


Figure  3.12:  Force  history  for  crashback  showing  total  force  in  the  x-direction,  case 

cbsub.2b. 


CHAPTER  IV 


DISCUSSION  OF  RESULTS 

The  first  step  in  evaluating  the  computational  results  is  to  validate  them  with 
available  experimental  data.  The  next  section  presents  the  validation  of  the  bare  hull, 
isolated  propeller,  and  forward  propelled  simulations  with  the  available  experimental 
data  described  in  Section  2.1.  A  comparison  of  the  two  turbulence  models  is  then 
given.  The  real  essence  of  this  chapter  lies  in  Sections  4.3  and  4.4  where  Reynolds 
scaling  and  hull/propulsor  flow  interaction  is  discussed.  Section  4.5  discusses  the 
results  of  the  crashback  simulation  and  presents  a  few  of  the  more  interesting  flow 
characteristics.  The  last  section  of  this  chapter  discusses  possible  sources  of  error  for 
the  computational  results. 

4.1  Validation  of  Computational  Results 
Drag  coefficient  is  of  paramount  interest  in  the  bare  hull  simulations. 
Theoretically,  as  the  Reynolds  number  increases,  the  drag  coefficient  should  decrease. 
This  occurs  because  the  boundary  layer  becomes  thinner  with  increase  in  Reynolds 
number,  thus  diminishing  the  effective  contour  the  flow  sees.  As  shown  in  Figure  4.1, 
this  trend  holds  true  for  the  computational  as  well  as  the  experimental  results. 
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Figure  4.1 :  Drag  coefficient  versus  Reynolds  number  with  15%  error  bars.  LCC 
model  scale  =  22.964  feet.  Normalized  by  length. 


It  is  also  noted  from  Figure  4.1  that  the  agreement  between  the  computational 
results  and  the  experimental  results  degrades  with  the  increase  in  Reynolds  number.  A 
potential  explanation  for  this  is  discussed  at  the  end  of  this  chapter  with  possible 
sources  of  error. 
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Figure  4.2:  Drag  coefficient  versus  Reynolds  number.  LCC  model  scale.  Normalized 

by  projected  frontal  area. 


Figure  4.2  is  similar  to  Figure  4.1,  except  that  the  drag  coefficient  is 
normalized  by  projected  frontal  area  rather  than  hull  length.  The  same  trend  is 
apparent  in  both  Figures.  Drag  coefficient  decreases  with  Reynolds  number,  and  the 
agreement  between  the  experimental  and  computational  results  degrades  with  increase 
in  Reynolds  number. 

The  drag  coefficients  from  the  bare.2f,  bare.4f,  bare.5a,  bare.5b,  bare.6a,  and 
bare.6b  cases  are  compared  against  the  experimental  drag  coefficient  in  Table  4.1. 

The  drag  coefficients  are  normalized  by  hull  length.  The  percent  contribution  from 
viscous  forces  to  the  drag  coefficient  is  also  given.  The  large  viscous  contribution 
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illustrates  the  need  to  include  the  viscous  terms  in  the  computation  of  these 
simulations  despite  the  increased  computational  cost. 

Table  4. 1 :  Drag  coefficient  comparison  between  grids  and  turbulence  models. 


Case  Name 

Drag  Coefficient 

%  Viscous  Forces 

Experimental 

0.0008274 

bare.  2  f 

0.0006546 

84 

bare.4f 

0.0006254 

83.5 

bare.  5  a 

0.0005838 

89.5 

bare.  5  b 

0.0005816 

88 

bare.  6a 

0.0006252 

82 

bare.  6b 

0.0006458 

83.5 

The  tightest  axial  spacing  and  boundary  layer  resolution  is  with  the  structured 
grid,  although  it  is  noted  that  the  structured  grid  had  the  poorest  agreement  with  the 
experimental  data.  The  structured  grid  was  also  run  with  a  structured  flow  solver 
using  the  k-s  turbulence  model  [ramesh].  The  drag  coefficient  computed  was  0.00056. 
The  computational  drag  coefficients  are  consistently  low.  Three  different  grid 
topologies  and  densities,  three  different  turbulence  models,  and  two  different  flow 
solvers  were  used  to  investigate  this  phenomenon.  Further  investigation  is  warranted, 
but  is  beyond  the  scope  of  this  study.  The  unstructured/structured  results  had  the  best 
agreement  and  will  be  used  in  further  analyses. 

The  next  two  plots  show  the  agreement  between  the  computational  and 
experimental  results  for  thrust  and  torque  coefficients  for  the  isolated  propeller. 


| - Experimental  ■  Spalart-Allrraras  a  g-orrega 


Figure  4.3:  Validation  of  thrust  coefficient  for  isolated  propeller. 


Figure  4.4:  Validation  of  torque  coefficient  for  isolated  propeller. 
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Although,  a  discussion  on  the  turbulence  models  follows  later,  it  is  noted  from 
these  plots  that  the  Spalart-Allmaras  one  equation  turbulence  model  maintains  better 
agreement  with  the  experimental  data  for  the  thrust  coefficient.  For  the  torque 
coefficient,  the  q-co  two-equation  turbulence  model  has  better  agreement  with  the 
experimental  data.  However,  there  is  not  a  significant  difference  between  either 
turbulence  model. 

With  good  agreement  for  the  bare  hull  and  isolated  propeller  computations,  the 
forward  propelled  cases  are  next.  Excellent  agreement  with  the  experimental  data  is 
exhibited  by  all  forward  cases.  Some  difference  is  noted  for  the  higher  advance 
coefficients. 

The  next  four  plots  show  the  propeller  thrust  coefficient.  The  computational 
data  are  plotted  against  the  experimental  data.  Four  tunnel  speeds  are  given  with 
varying  advance  coefficients  for  each.  The  design  advance  coefficient  for  the  P4381 
is  J  =  0.889.  It  is  clear  that  performance  degrades  with  increase  in  advance 
coefficient,  and  for  advance  coefficients  less  than  the  design  advance  coefficient,  no 
real  advantage  is  gained  because  the  thrust  curve  begins  to  flatten  out.  It  is  noted  that 
the  propeller  actually  produces  drag  at  the  higher  advance  coefficients.  Axial  velocity 
contours  can  be  found  in  Appendix  B  for  each  speed  and  advance  coefficient. 


-Experimental  ■  Computational  j 


Figure  4.5:  Thrust  coefficient  for  installed  propeller  versus  advance  coefficient. 
Tunnel  speed  is  5  knots.  Reynolds  number  =  16  million. 
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Figure  4.6:  Thrust  coefficient  for  installed  propeller  versus  advance  coefficient 
Tunnel  speed  is  10  knots.  Reynolds  number  =  32  million. 
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4.2  Turbulence  Model  Comparisons 

The  choice  of  turbulence  model  has  a  direct  effect  on  the  accuracy  of  the 
solution.  Chapter  2  gave  a  description  of  the  two  turbulence  models  available  for  use. 
In  general,  the  Spalart-Allmaras  turbulence  model  is  more  forgiving  in  relation  to  the 
computational  geometry.  The  q-co  model  is  more  accurate,  but  due  to  number  of 
equations  to  be  solved  is  more  expensive.  It  also  tends  to  be  very  sensitive  to  the 
quality  of  the  surface  grid. 

It  was  found  that  the  propelled  simulations  would  not  run  using  the  q-co  model. 
It  is  thought  that  the  surface  grid  on  the  propeller  is  the  culprit.  The  propeller  blade 
tips  are  very  fine  and  require  many  points  for  resolution.  Surface  triangulation  by  the 
grid  generator  sometimes  creates  minute  ‘divots’  in  the  surface  topology,  particularly 
for  sharp  edges.  Due  to  the  sensitivity  of  the  q-co  model,  these  ‘divots’  appear  to  be 
part  of  the  geometry  itself.  In  trying  to  account  for  these  geometry  characteristics,  the 
turbulence  model  becomes  unstable  and  blows  up  in  relatively  short  order.  With 
special  attention  paid  to  the  propeller  tips,  the  geometry  could  be  refined  so  the  q-co 
model  could  be  used.  This  was  deemed  unnecessary  since  the  Spalart-Allmaras  model 
yielded  accurate  results.  Also,  the  amount  of  time  required,  and  the  additional 
computational  resources  that  would  be  needed  were  considered  too  great  for  this 
study.  As  shown  in  the  previous  section,  there  is  not  a  significant  difference  between 
the  results  from  either  model.  The  agreement  between  computational  and 
experimental  results  for  the  forward  propelled  solutions  is  extraordinary. 
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Both  models  were  used  for  the  isolated  propeller  and  bare  hull  simulations.  As 
stated  earlier,  both  models  show  excellent  agreement  with  the  experimental  data  for 
the  isolated  propeller.  For  the  bare  hull  simulations,  the  q-co  has  slightly  better 
agreement  with  the  experimental  results,  particularly  at  the  lower  Reynolds  numbers. 
This  agreement  degrades  with  the  increase  in  Reynolds  number.  For  the  higher 
Reynolds  numbers,  the  results  from  both  turbulence  models  are  off  by  15-20%. 

It  is  important  to  understand  the  differences  in  the  flow  field  due  to  the 
turbulence  models.  The  Reynolds  scaling  study  uses  the  q-co  model  since  it  had  the 
best  agreement  with  the  experimental  data.  The  hull/propulsor  interaction  study  uses 
the  Spalart-Allmaras  model  because  of  its  use  in  the  propelled  simulations.  Given 
below  are  the  pressure  and  skin  coefficient  distributions  across  the  hull  for  each 
turbulence  model.  The  first  plot,  Figure  4.9,  shows  that  both  turbulence  models  yield 
approximately  the  same  pressure  distribution.  Although  in  Figure  4. 10,  it  is  clear  that 
the  Spalart-Allmaras  model  has  a  lower  skin  friction  distribution  than  does  the  q-co 
model. 
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Figures  4. 1 1  and  4,12  show  the  xy-plane  axial  velocity  contours  at  a  Reynolds 
number  of  32  million  for  each  turbulence  model. 


O  TT/TImr  X 

Figure  4.1 1:  Xy-plane  axial  velocity  contour,  case  bare.4d;  Spalart-Allmaras. 
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Figure  4. 12:  Xy-plane  axial  velocity  contour,  case  bare.2d;  q-co. 
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The  xy-plane  axial  velocity  contours  for  both  turbulence  models  are  very  similar. 
Upon  close  inspection,  it  is  found  that  the  q-ca  solution  has  a  thicker  boundary  layer. 
The  next  two  figures  are  velocity  contours  in  the  yz-plane  at  an  x/L  location  of  0.9987. 
The  view  is  from  the  rear  of  the  submarine  looking  forward. 


Figure  4.13:  Velocity  contour  at  x/L  =  0.9987  in  yz-plane,  case  bare.4d;  Spalart- 

Allmaras. 
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Figure  4.14:  Velocity  contour  at  x/L  =  0.9987  in  yz-plane,  case  bare.2d;  q-ro. 
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The  difference  in  die  boundary  layer  thickness  can  more  easily  be  seen  in  the 
yz-plane  axial  velocity  contours  than  in  the  xy-plane  axial  velocity  contours.  This 
difference  is  small,  but  consistent  with  the  drag  coefficient  plots.  The  thicker 
boundary  layer  would  produce  more  drag. 

The  isolated  propeller  cases  were  also  run  with  the  two  turbulence  models. 

The  next  few  plots  show  the  xy-plane  and  yz-plane  axial  velocity  contours  for  the 
isolated  P4381  at  an  advance  coefficient  of  J  =  0.95.  Flow  field  differences  are  minute 
at  best.  This  is  consistent  with  Figures  4.3  and  4.4;  there  is  not  much  difference 
between  the  two  turbulence  models. 
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Figure  4. 1 5:  Xy-plane  axial  velocity  contour,  J  =  0.95;  Spalart-Allmaras. 
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Figure  4.16:  Xy-plane  axial  velocity  contour,  J  -  0.95;  q-co. 
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Figure  4. 1 7:  Yz-plane  axial  velocity  contour,  J  =  0.95;  Spalart-Allmaras. 
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Figure  4. 1 X:  Yz-plane  axial  velocity  contour,  .1  =  0.95;  q-w. 

As  can  be  seen  in  Figures  4. 1  1  4.1  X,  the  differences  in  specific  flow 

characteristics  between  the  turbulence  models  are  small.  The  most  noticeable 
difference  overall  is  with  the  skin  friction  distribution.  It  is  decided  that  consistency  in 
the  use  of  a  turbulence  model  is  more  important  than  which  particular  model  is  used. 

4.3  Reynolds  Scaling  Analysis 

This  section  aims  to  explore  some  of  the  particular  differences  in  flow 
characteristics  between  Reynolds  scales.  The  three  scales  under  investigation  are  the 
small  model  scale  (1 1.27X  inches),  the  LCC  model  scale  (22.964  feet),  and  full  scale 
(3X3.5  feet).  These  three  scales  generate  Reynolds  numbers  from  770,000  to  1.6 
billion.  It  has  previously  been  discussed  that  simply  applying  Reynolds  scaling 
principles  to  the  force  components  is  not  sufficient  because  it  does  not  take  into 
account  localized  differences  [paper].  The  first  plot  (Figure  4.19)  given  is  drag 
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coefficient  versus  Reynolds  number.  This  is  similar  to  Figure  4.1,  but  includes  only 
the  computational  results  and  covers  a  broader  range  of  Reynolds  numbers.  As 
expected,  the  same  trend  of  drag  coefficient  decreasing  with  Reynolds  number 
increase  is  readily  seen.  It  is  noted  that  the  most  significant  drop  happens  across  the 
lower  Reynolds  numbers. 


Small  Model  (0.94  ft)  — B—  LCC  Model  (22.964  ft)  -  -  a-  - Full  Scale  (383.5  ft)  | 

Figure  4.19:  Computational  drag  coefficient  versus  Reynolds  number,  three  length 

scales.  Normalized  by  length. 

Table  4.2  is  a  table  of  bare  hull  xy-plane  axial  velocity  contours.  Presenting 
these  contours  side  by  side  makes  it  easier  to  see  the  differences  between  the  Reynolds 
scales.  As  the  Reynolds  number  increases,  the  boundary  layer  becomes  thinner.  The 
Reynolds  number  and  case  number  are  given  below  the  contour.  It  is  noted  that  the 
bare.3x  contours  have  some  peculiar  characteristics.  An  explanation  is  offered  in 


Section  4.5. 


Table  4.2:  Xy-plane  axial  velocity  contours. 
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Table  4.3  compares  the  yz-plane  axial  velocity  contours  for  all  the  Reynolds 
numbers.  Again,  the  progressive  diminishing  of  the  boundary  layer  can  be  seen  across 
the  Reynolds  numbers.  The  view  is  aft  looking  forward. 


Table  4.3:  Yz-plane  axial  velocity  contours. 
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The  boundary  layer  diminishing  process  seems  to  become  more  rapid  and  more 
pronounced  at  the  higher  Reynolds  numbers.  This  seems  counter  to  the  drag 
coefficient  plot,  Figure  4.19,  where  the  decrease  in  drag  coefficient  is  more  rapid  and 
pronounced  at  the  lower  Reynolds  number.  As  stated  previously,  the  diminishing 
boundary  layer  and  decreasing  drag  coefficient  go  hand  in  hand. 


Figure  4.20:  Skin  friction  coefficient  distribution  along  bare  hull  for  various  Reynolds 

numbers. 

As  seen  in  Figure  4.20,  the  skin  friction  coefficient  also  decreases  with 
increase  in  Reynolds  number.  The  decrease  in  skin  friction  coefficient  is  indicative  of 
the  decrease  in  boundary  layer  thickness.  This  is  related  to  the  corresponding  decrease 
in  drag  coefficient.  Note  that  there  is  a  lot  of  noise  present  at  the  lowest  Reynolds 
number.  This  may  be  due  to  the  off-wall  grid  spacing  being  too  tight  for  that  low  of  a 


Reynolds  number.  Since  one  grid  was  used  for  all  the  Reynolds  numbers,  the  lowest 
Reynolds  numbers  had  the  tightest  off-wall  spacing. 

Next  is  to  compare  the  propeller  flow  field  at  different  Reynolds  numbers,  but 
for  the  same  advance  coefficient.  It  has  been  shown  in  Figure  3.10  that  similar 
advance  coefficients  have  similar  thrust  coefficients.  The  two  xy-plane  axial  velocity 
contours  and  the  two  yz-plane  axial  velocity  contours  given  below  compare  the 
propeller  flow  fields  of  an  isolated  P4381  propeller  for  an  advance  coefficient  of 
J-0.95  at  Reynolds  numbers  of  2.4  million  and  32  million.  For  the  yz-plane  contour 
pictures,  the  view  is  aft  looking  forward.  The  cut  is  at  an  x/L  of  0.9987,  which  is  at 
the  rear  of  the  propeller  blades. 


Figure  4.21(a):  J  =  0.95,  Figure  4.21(b):  J  -  0.98, 

Re  =  2.4  million  Re  =  32  million 
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Figure  4.22(a):  J  =  0.95, 
Re  =  2.4  million. 


Figure  4.22(b):  J  =  0.98, 
Re  =  32  million. 


For  the  Reynolds  number  of  2.4  million,  the  propeller  has  a  greater  influence 
on  the  surrounding  flow  field.  There  is  a  ‘bulge’  of  high-speed  flow  in  between  blade 
passages.  The  region  of  high-speed  flow  is  more  confined  in  the  higher  Reynolds 
number  case.  This  confinement  limits  the  propeller’s  region  of  influence.  This  effect 
is  very  similar  to  the  reduction  in  boundary  layer  thickness  of  the  bare  hull  cases. 
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Figure  4.23(a):  Pressure  contours  for  Figure  4.23(b):  Pressure  contours  for 
stem  region  at  Re  =  16  million.  stem  region  at  Re=  64  million. 

Figures  4.23(a)  and  4.23(b)  show  the  differences  in  the  region  of  influence  of 
the  stem  appendages  for  two  Reynolds  numbers.  The  view  is  from  the  top  looking 
down.  The  stagnation  point  on  the  rudder  at  a  Reynolds  number  of  64  million  exhibits 
a  greater  upstream  influence  than  at  a  Reynolds  number  of  16  million.  The 
downstream  influence  is  also  more  pronounced.  The  effect  of  the  stem  appendages  is 
more  widely  felt  at  the  higher  Reynolds  number  because  of  the  decreased  boundary 
layer  thickness. 

The  results  given  in  this  section  illustrate  the  concept  of  boundary  layer 
thickness  reduction,  and  thusly  reduction  in  drag  coefficient,  with  increase  in 
Reynolds  number.  It  has  been  shown  that  for  the  hull  and  the  propeller,  the  Reynolds 
number  has  a  localized  effect  on  flow  field  characteristics. 


4.4  Hull/Propulsor  Interaction 
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The  complex  flow  surrounding  a  powered  submarine  is  analyzed  by  breaking 
the  flow  up  into  its  component  pieces  and  examining  the  characteristics  unique  to  each 
component.  A  discussion  of  the  flow  field  interactions  due  to  the  hull,  appendages 
and  propeller  ensues. 

The  three  figures  below  show  graphically  the  differences  in  the  xy-plane  axial 
flow  field  due  to  the  stem  appendages  and  the  propeller  at  a  Reynolds  number  of  32  x 
106.  The  addition  of  the  stern  appendages  thickens  the  boundary  layer  upstream,  but 
lessens  the  influence  of  the  hull  boundary  layer  downstream  of  the  submarine.  The 
stem  appendages  increase  the  drag  on  the  hull  by  125%.  The  drag  coefficient  of  the 
bare  hull  is  0.00074;  with  the  addition  of  the  stern  appendages,  the  drag  coefficient 
increases  to  0.00 1 7.  In  the  third  picture,  the  boundary  layer  can  be  seen  “necking 
down”  just  before  the  propeller  hub.  The  propeller  drastically  affects  the  submarine’s 
downstream  influence,  but  doesn’t  appear  to  have  any  influence  upstream  of  the  stem 
appendages.  It  is  apparent  that  the  propeller  is  immersed  in  the  wake  of  the 


appendages. 
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Figure  4.24:  Influence  of  stern  appendages  and  propeller  on  the  hull  velocity  field. 

The  next  figure  shows  the  upstream  influence  of  the  stern  appendages  and  the 
propeller.  The  view  is  a  cut  in  the  yz-plane  just  before  the  appendages.  The  axial 
location  is  x/L  =  0.887.  The  view  is  forward  looking  aft. 


Hull  with  Stem  Appendages 


Hull  with  Stem  Appendages  and  Propeller 


Figure  4.25:  Upstream  influence  of  the  stem  appendages  and  propeller. 
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It  is  clearly  seen  in  Figure  4.25  that  the  stem  appendages  have  an  upstream 
influence.  The  lack  of  noticeable  differences  in  the  flow  field  of  the  hull  with 
appendages  and  the  hull  with  appendages  and  propeller  shows  that  the  appendages  are 
the  predominant  upstream  influence.  The  propeller  has  little,  if  any,  effect  on  the  flow 
upstream  from  the  appendages.  This  may  be  due  to  the  propeller  being  slightly 
undersized  and  further  downstream  than  is  typical. 

Figure  4.26  shows  that  the  stem  appendages  maintain  a  significant  influence 
downstream  of  their  position.  This  view  is  at  approximately  the  middle  of  the 
propeller.  The  axial  location  is  x/L  =  0.986  and  the  view  is  aft  looking  forward. 


Hull  with  Stem  Appendages  and  Propeller 
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Figure  4.26:  Downstream  influence  of  the  stem  appendages  and  propeller. 
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The  wake  from  the  stem  appendages  squares  off  the  interior  velocity  field 
downstream.  Their  location  can  be  identified  by  the  cross  pattern  in  the  flow  field. 

The  swirl  pattern  seen  in  the  third  picture  is  due  to  the  propeller  blade  passages.  The 
white  area  is  the  accelerated  flow  that  provides  the  thrust. 

The  next  series  of  graphics  are  a  side-by-side  comparison  of  the  velocity  fields 
for  the  isolated  propeller  and  the  installed  propeller.  Figure  4.26  shows  yz-plane  axial 
contours  at  the  rear  of  the  propeller  blades  with  the  view  aft  looking  forward.  Both 
solution  sets  were  ran  at  the  same  conditions  for  a  Reynolds  number  of  32  x  106.  The 
contour  scales  have  been  adjusted  to  show  the  most  flow  characteristics.  The  advance 
coefficient  and  case  name  are  given  below  each  picture. 


Table  4.4:  Comparison  of  isolated  propeller  and  installed  propeller  yz-plane  axial 
velocity  fields  showing  the  influence  of  the  appendage  wakes  and  hull  boundary  layer. 


Table  4.4  (continued):  Comparison  of  isolated  propeller  and  installed  propeller  yz- 
plane  axial  velocity  fields  showing  the  influence  of  the  appendage  wakes  and  hull 

boundary  layer. 


-0.17  U/Uzvf  n.ft-2 


memm- 

-0.3  U/Urcf  1^15 


J  =  1,3,  isoprop.3c 


J  =1.3,  fwdsub.2c 
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Table  4.4  (concluded):  Comparison  of  isolated  propeller  and  installed  propeller  yz- 
plane  axial  velocity  fields  showing  the  influence  of  the  appendage  wakes  and  hull 

boundary  layer. 


The  cross  pattern  due  to  the  appendage  wakes  is  apparent  in  several  of  the 
above  contours.  For  J  =  1 .3,  the  appendage  wakes  and  hull  boundary  layer  create 
almost  a  low  velocity  ring  around  the  propeller.  This  phenomenon  is  not  present  in 
the  corresponding  isolated  propeller  contour.  The  combination  of  the  wakes  and  the 
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hull  boundary  layer  confine  the  flow  field  generated  by  the  propeller.  All  cases  show 
similar  flow  characteristics  between  the  two  configurations.  The  tip  vortices  are  very 
noticeable  in  both  cases  for  J  =  3.83.  The  big  difference  between  the  two  J  =  3.83 
contours  is  the  immersion  of  the  installed  propeller  in  the  hull  and  appendage 
boundary  layers.  The  interior  propeller  flow  field  is  mostly  low  velocity.  On  the  other 
hand,  the  isolated  propeller  interior  flow  field  has  pockets  of  high  velocity  flow. 

Notice  for  J  =  1.3,  the  impingement  of  the  appendage  wakes  on  the  propeller  flow 
field  smears  out  the  low  velocity  pockets  to  a  blurred  low  velocity  ring.  The  presence 
of  the  appendage  wakes  is  less  apparent  at  the  lower  advance  coefficients,  but  it  still 
has  an  effect.  The  isolated  propeller  contour  for  J  =  0.78  is  fairly  round.  The 
appendage  wakes  and  hull  boundary  layer  constrain  the  propeller  flow  field,  which  is 
exhibited  by  the  flattened  sides  of  the  high  velocity  region. 

Table  4.5  is  a  side-by-side  comparison  of  the  xy-plane  axial  velocity  contours 
for  the  isolated  propeller  and  the  installed  propeller.  It  is  interesting  to  note  that  many 
of  the  same  flow  characteristics  appear  in  both  sets  of  cases. 
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Table  4.5:  Comparison  of  isolated  propeller  and  installed  propeller  xy-plane  axial 
velocity  fields  showing  the  influence  of  the  appendage  wakes  and  hull  boundary  layer. 


Table  4.5  (continued):  Comparison  of  isolated  propeller  and  installed  propeller  xy- 
plane  axial  velocity  fields  showing  the  influence  of  the  appendage  wakes  and  hull 

boundary  layer. 
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Table  4.5(concluded):  Comparison  of  isolated  propeller  and  installed  propeller  xy- 
plane  axial  velocity  fields  showing  the  influence  of  the  appendage  wakes  and  hull 

boundary  layer. 


The  predominant  difference  between  the  two  flow  fields  is  the  inflow 
conditions.  The  isolated  propeller  has  close  to  free  stream  inflow  conditions  while  the 
propeller  on  the  body  has  the  benefit  of  a  reduced  inflow  due  to  the  hull  and 
appendage  boundary  layers.  Another  note  is  that  the  isolated  propeller,  even  at  the 
lower  advance  coefficients  never  generates  a  concentrated  jet  of  flow.  The  reduced 
inflow  and  the  containment  by  the  body  and  appendage  boundary  layers  allow  for  the 
body-mounted  propeller  to  produce  that  jet  of  flow. 

It  is  becoming  clear  that  rather  than  hurting  the  performance  of  the  propeller, 
the  hull  and  appendage  boundary  layers  actually  provide  the  necessary  environment 
for  the  propeller  to  do  its  job.  This  makes  sense  given  that  marine  propellers  are 
designed  for  optimum  performance  with  a  body. 
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To  support  the  previous  statement  further,  Figure  4.27  compares  the  thrust 
coefficients  computed  from  the  isoprop. 3x  and  fwdsub.2x  cases.  As  shown  in 
Appendix  B,  both  sets  of  cases  are  run  at  identical  conditions. 


Figure  4.27:  Thrust  coefficient  comparison  between  P4381  mounted  on  appended 
LCC  hull  and  isolated  P4381;  Re  =  32  million. 

Note  in  Figure  4.27,  the  propeller  mounted  on  the  body  performs  better  than 
the  isolated  propeller.  At  first,  this  may  seem  backwards,  but  as  stated  previously, 
marine  propulsors  are  designed  to  operate  mounted  to  a  larger  body.  The  propeller  on 
the  appended  LCC  hull  is  buried  in  the  boundary  layer  from  the  hull  and  the 
appendages.  Due  to  this,  the  effective  inflow  that  the  mounted  propeller  sees  is  less 
than  that  of  the  isolated  propeller.  The  reduced  inflow  velocity  increases  the  lift  and 
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therefore  the  thrust  produced  by  the  propeller  blades.  This  graph  illustrates  the  point 
that  the  hull  and  appendage  boundary  layers  actually  aid  the  propeller  in  performing 
its  function. 

4.4  Submarine  Crashback 

In  this  section,  the  crashback  computations  will  be  shown  to  be  valid  and  a  few 
of  the  more  interesting  flow  characteristics  will  be  examined.  The  first  plot  is  of  the 
total  axial  force  versus  advance  coefficient  for  both  the  cback.lx  and  cback.2x  cases. 


-^—Experimental  -  5  Knots  —^Experimental  - 1  0  Knots 

±  Computational  -  5  Knots  •  Computational  -  10  Knots 

Figure  4.28:  Validation  of  crashback  computations  cback.lx  and  cback.2x. 
Experimental  has  1 5%  error  bars. 

The  results  for  the  ten  knot  cases  are  much  better  than  the  results  for  the  five 
knot  cases.  The  computational  results  were  averaged  over  two  propeller  revolutions  in 
an  attempt  to  attain  the  best  possible  force  representation.  Most  of  the  results  are 


64 

within  the  15%  error  bound  which  is  very  good  agreement  for  such  an  unsteady  case. 
The  ring  vortex  presents  a  unique  challenge  in  simulating  this  flow;  it  is  constantly 
changing.  Figure  4.29  is  an  intrinsic  swirl  isosurface  of  value  ~2  [19].  The  propeller 
blades  are  colored  by  pressure.  Notice  the  asymmetry  of  the  ring  vortex.  Some  of  the 
asymmetry  is  due  to  the  nature  of  the  ring  vortex  itself,  but  it  is  believed  that  the 
asymmetry  of  the  stern  appendages  and  resulting  wakes  have  some  influence  on  the 
ring  vortex  structure  as  well. 


Figure  4.29:  Swirl  parameter  isosurface  showing  the  instantaneous  ring  vortex 

structure. 


Figures  4.30  and  4.3 1  show  the  ring  vortex  in  relation  to  the  entire  submarine. 
Figure  4.30  is  another  isosurface  picture  showing  the  establishment  of  the  ring  vortex 
at  the  rear  of  the  submarine.  Figure  4.3 1  is  a  particle  trace.  The  swirl  in  the  trace  is 
very  descriptive  of  the  nature  of  the  ring  vortex. 


Figure  4.30:  Ring  vortex  in  relation  to  submarine. 


Figure  4.3 1 :  Particle  trace  of  ring  vortex. 
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The  next  set  of  figures  is  a  comparison  between  the  crashback  velocity 
contours  and  a  similar  set  of  forward  propelled  contours.  The  crashback  velocity 
contours  are  for  the  cbsub.2c  case.  The  related  forward  case  is  the  fwdsub.2c  case. 


The  first  two  pictures.  Figure  4.32,  are  yz-plane  axial  velocity  contours  of  the 
propeller.  Figure  4.33  is  a  set  of  xy-plane  axial  velocity  contours  showing  the 


reversed  flow  region  around  the  propeller  of  the  crashback  case. 
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cbsub.2c  J  =  - 1 .45  fwdsub.2c  J  =  1.3 

Figure  4.32:  Yz-plane  axial  velocity  contour  comparison. 


cbsub.2c  J  =  -1.45 


fwdsub.2c  J  =  1.3 


Figure  4.33:  Xy-plane  axial  velocity  contour  comparison. 
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There  are  several  things  to  note  in  the  above  velocity  contours.  The 
impingement  on  the  boundary  layer  upstream  of  the  propeller  is  clearly  shown  in  the 
crashback  contour  of  Figure  4.33.  In  the  crashback  contour  of  Figure  4.32,  the 
instability  of  the  flow  field  and  resulting  asymmetry  can  be  seen.  The  propeller  in 
crashback  also  has  a  much  larger  region  of  influence  than  does  the  propeller  rotating 
forward.  The  crashback  flow  field  is  very  complex  and  is  an  intriguing  flow  to  model. 
The  results  presented  in  this  section  show  good  agreement  with  the  experimental 
results  and  several  of  the  pictures  clearly  show  the  physical  phenomenon  associated 
with  crashback  flow. 

4.6  Possible  Sources  of  Error 

There  are  a  couple  of  possible  sources  of  error.  As  with  any  study,  the  results 
are  not  infallible,  but  considering  the  good  to  excellent  agreement  with  the 
experimental  results,  the  overall  quality  of  the  computations  is  felt  to  be  very  high. 

Considering  the  bare  hull  drag,  one  source  of  error  might  be  the  absence  of  the 
test  strut.  At  the  onset  of  the  study,  it  was  decided  that  the  strut  was  probably  a 
negligible  source  of  drag.  But  it  is  recognized  that  the  strut  would  have  a  contribution 
to  the  recorded  drag.  Also,  the  tunnel  walls  were  not  modeled.  All  of  the 
computations  were  run  as  open  water  simulations.  It  was  decided  that  the  tunnel 
blockage  was  not  sufficient  as  to  warrant  the  computational  resources  that  would  be 
required  for  such  a  large-scale  simulation. 

Another  possible  source  would  be  the  geometry  definition  for  the  hull.  This 
possibility  did  not  show  up  until  the  full-scale  simulations  were  run.  As  shown  in 
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Figure  4.34,  with  the  highest  Reynolds  numbers  (as  in  those  for  the  full  scale 
submarine),  a  numerical  artifact  is  present  in  the  velocity  flow  field.  This  artifact  is 
represented  as  the  dark  blue  area  in  the  boundary  layer  region.  Although  this  artifact 
affects  the  boundary  layer  region  of  the  velocity  flow  field,  it  doesn’t  seem  to  have  a 
major  effect  on  the  overall  flow  field. 


Figure  4.34:  Numerical  artifact  in  high  Reynolds  number  flow  boundary  layer. 

Re  =  1  billion. 

When  the  curve  definition  for  the  hull  is  analyzed  by  taking  the  second 
derivative,  it  is  noted  that  the  curve  is  not  entirely  smooth.  The  top  line  is  the  hull 


definition  and  the  bottom  line  is  the  second  derivative. 
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Figure  4.35:  Hull  curve  and  2nd  derivative.  Large  plot  is  close  up  of  stem  region. 

The  ‘wiggles’  apparent  in  the  hull  definition  may  have  caused  some  bumps  in 
the  computational  model  of  the  hull.  While  at  the  lower  Reynolds  number  this  didn’t 
seem  to  have  much  of  an  effect,  it  did  affect  the  boundary  layer  for  the  higher 
Reynolds  numbers.  At  the  higher  Reynolds  numbers,  the  boundary  layer  is  so  thin 
that  the  flow  becomes  much  more  sensitive  to  the  computational  geometry.  To  prove 
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that  the  hull  definition  is  the  culprit,  the  curve  would  need  to  be  smoothed,  a  new 
computational  model  generated,  and  the  solution  run  again. 

The  last  source  to  be  discussed  is  the  grid  resolution  behind  the  stem 
appendages.  Grid  refinement  to  the  rear  of  the  stem  appendages  would  help  resolve 
the  wakes.  The  tighter  resolution  of  the  wakes  might  yield  some  more  interesting  flow 
interactions  with  the  propeller. 


CHAPTER  V 


CONCLUSIONS  AND  RECOMMENDATIONS 

5.1  Conclusions 

Computational  simulation  is  growing  by  leaps  and  bounds.  Technological 
advances  are  making  it  easier  and  faster  to  model  complex  geometries  and  physics. 
This  study  has  shown  that  the  flow  solver,  U2NCLE,  can  produce  valid  results  for 
isolated  propeller,  bare  hull,  forward  propelled,  and  crashback  submarine  simulations. 
With  the  exception  of  the  bare  hull  configuration,  all  simulations  produced  solutions 
with  excellent  agreement  to  experimental  data.  Although  the  agreement  between  the 
experimental  drag  coefficient  and  the  computational  drag  coefficient  for  the  bare  hull 
would  not  be  considered  excellent,  it  would  still  be  considered  good  agreement. 

The  flow  solver  has  also  been  shown  to  run  long  duration  unsteady  problems, 
such  as  1 8-20  propeller  revolutions  for  a  submarine  in  crashback.  Previous  propelled 
studies  have  only  run  8-10  revolutions. 

With  respect  to  the  turbulence  models,  it  has  been  determined  that  the 
consistency  of  turbulence  model  use  is  more  important  than  which  particular  model  is 
used.  While  the  q-co  turbulence  model  generally  produced  better  agreement  with 
experimental  data,  it  was  also  more  difficult  to  run.  Its  sensitivity  to  the  surface  grid 
can  make  it  difficult  to  work  with  depending  upon  the  type  of  geometry  being 
modeled. 
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The  importance  of  the  computational  geometry  itself  has  also  been  pointed  out. 
The  Spalart-Allmaras  model  has  been  shown  to  be  both  forgiving  and  reliable.  The 
results,  particularly  for  the  forward  propelled  simulations,  produced  with  the  Spalart 
model  all  had  good  to  excellent  agreement  with  the  experimental  results. 

The  Reynolds  scaling  study  showed  that  the  drag  coefficient  for  the  bare  hull 
decreased  with  the  increase  in  Reynolds  number.  Related  to  this,  the  boundary  layer 
became  increasingly  thin  with  Reynolds  number  as  well.  The  change  in  the  velocity 
profiles  is  something  that  cannot  be  taken  into  consideration  with  Reynolds  scaling  of 
forces.  Depending  on  which  flow  characteristics  are  of  great  importance, 
computationally  simulating  the  entire  flow  field  becomes  essential,  particularly  with 
respect  to  the  propeller.  Reynolds  number  affects  the  extent  of  the  propeller’s 
influence  on  the  surrounding  flow.  The  propeller  flow  field  has  some  localized  flow 
characteristics  that  cannot  be  accounted  for  through  basic  Reynolds  scaling.  Thus,  the 
need  to  computationally  simulate  the  entire  flow  field  for  higher  Reynolds  numbers 
has  been  supported. 

It  has  been  shown  that  the  stem  appendages  have  the  greatest  influence  on  both 
the  hull  and  propeller  flow  fields.  The  propeller  is  more  efficient  mounted  on  a 
submarine  than  isolated  by  itself.  By  itself,  it  actually  creates  more  drag  because  of 
the  effective  inflow.  The  boundary  layer  interaction  between  the  hull,  appendages  and 
propeller  is  actually  a  good  thing.  It  allows  the  propeller  to  take  advantage  of  a 
reduced  inflow  thereby  increasing  its  thrust  production. 
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The  submarine  crashback  has  been  discussed  and  pictorially  described. 

Specific  flow  characteristics  have  been  pointed  out,  and  insights  gained  in  the  nature 
of  this  very  complex,  and  unsteady  flow.  The  validation  of  the  code  to  run  crashback 
simulations  is  also  very  important. 

This  study  has  made  a  significant  advance  towards  the  accurate  prediction  of 
forces  and  flow  characteristics  for  a  fully  submerged  body.  It  is  one  of  the  first  studies 
to  validate  fully  submerged  and  propelled  simulations  with  experimental  data.  This 
study  has  hopefully  shed  a  little  light  on  understanding  the  physics  associated 
submarine  flow  fields.  This  understanding  might  help  with  the  design  of  a  faster, 
quieter  submarine. 

5.2  Recommendations 

Several  recommendations  are  made: 

•  Complete  a  grid  topology  and  density  study  on  the  bare  hull  drag 
coefficient  to  determine  why  the  agreement  with  experimental  data 
degrades  with  increase  in  Reynolds  number. 

•  Smooth  hull  definition  and  re-run  full-scale  bare  hull  simulations. 

•  Run  additional  crashback  simulations  to  see  if  agreement  is  consistent. 

•  Refine  the  propeller  tips  to  run  a  full-scale  propelled  simulation. 
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TABLE  OF  CONDITIONS  FOR  EACH  COMPUTATIONAL  CASE 
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Case  Name  Scale*  Reynolds#  Velocity  1  Propeller 

(f 


_ 78_ 

Turbulence 

Model 


isoprop. 2a 


isoprop.2b 


isoprop. 2c 


isoprop.  3  a 


isoprop.3b 


isoprop. 3c 


isoprop. 3d 


isoprop. 3e 


fwdsub.la 

LCC 

fwdsub.lb 

LCC 

fwdsub.lc 

. 

LCC 

fwdsub.2a 

LCC 

fwdsub.2b 

LCC 

fwdsub.2c 

LCC 

fwdsub.2d 

LCC 

fwdsub.2e 

LCC 

fwdsub.3a 

LCC 

fwdsub.3b 

LCC 

fwdsub.3c 

LCC 

6.93 


7.41 


16.87 


16.87 


16.87 


16.87 


16.87 


Submarine  (appended  and  propelled  LCC  hull  -  P4381) 


Spalart 

Allmaras 


Spalart 

Allmaras 


Spalart 

Allmaras 


Spalart 

Allmaras 


Spalart 

Allmaras 


Spalart 

Allmaras 


Spalart 

Allmaras 


Spalart 

Allmaras 


8.42 

3.33 

Spalart 

Allmaras 

8.42 

4.32 

Spalart 

Allmaras 

8.42 

8.63 

Spalart 

Allmaras 

16.87 

4.25 

Spalart 

Allmaras 

16.87 

8.5 

Spalart 

Allmaras 

16.87 

12.75 

Spalart 

Allmaras 

16.87 

17 

Spalart 

Allmaras 

16.87 

21.25 

Spalart 

Allmaras 

25.32 

6.35 

Spalart 

Allmaras 

25.32 

12.68 

Spalart 

Allmaras 

25.32 

25.38 

Spalart 

Allmaras 
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Case  Name 

Scale* 

Reynolds  # 

Velocity 

(fps) 

Propeller 

rPs _ 

Turbulence 

Model 

fwdsub.4a 

LCC 

64  x  106 

33.75 

8.43 

Spalart 

Allmaras 

fwdsub.4b 

LCC 

64  x  106 

33.75 

16.88 

Spalart 

Allmaras 

fwdsub.4c 

LCC 

64  x  \06 

33.75 

28.33 

Spalart 

Allmaras 

cbsub.la 

LCC 

16  x  106 

8.42 

-3.35 

Spalart 

Allmaras 

cbsub.lb 

LCC 

16  x  106 

8.42 

-8.33 

Spalart 

Allmaras 

cbsub.2a 

LCC 

32  x  106 

16.87 

-3.33 

Spalart 

Allmaras 

cbsub.2b 

LCC 

32  x  106 

16.87 

-8.33 

Spalart 

Allmaras 

cbsub.2c 

LCC 

32xlOfl 

16.87 

-11.67 

Spalart 

Allmaras 

*  Small:  hull  length  is  1 1.278  inches,  LCC:  hull  length  is  22.964  feet,  Full:  hull  length  is  383.5  feet 
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VELOCITY  CONTOURS  FOR  FORWARD  PROPELLED  SIMULATIONS 
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1  45  U/Ur«I  -0  3 

fwdsub.2c;  Yz-plane  axial  velocity;  J=  1 .3 


145  O/Urcl  -0.3 

fwdsub.2c;  Xy-plane  axial  velocity;  J=1.3 


Figure  B.2:  Velocity  contours;  10  knot  tunnel  speed;  Re  =  32  million. 
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1.3  U/Urcf  -0.6 

fwdsub.2d;  Yz-plane  axial  velocity;  J=0.98 


1.3  U/Uref  -0  6 

fwdsub.2d;  Xy-plane  axial  velocity;  J=0.98 


fwdsub.2e;  Yz-plane  axial  velocity;  J=0.78  fwdsub.2e;  Xy-plane  axial  velocity;  J=0.78 


Figure  B.2  (concl’d):  Velocity  contours;  10  knot  tunnel  speed;  Re  =  32  million. 
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Pi 
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mm  ■ 


lwdsiib.3u;  Yz-plane  axial  velocity;  J-3.88 


lwdsub.3b;  Yz-plane  axial  velocity;  J=1.95 


fwdsub.3c;  Xy-plane  axial  velocity;  J=0.98 


fwdsub.3c;  Yz-plane  axial  velocity;  J=0.98 


Figure  B.3:  Velocity  profiles;  15  knot  tunnel  speed;  Re  =  48  million. 
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Figure  B.4:  Velocity  contours;  20  knot  tunnel  speed;  Re  =  64  million. 
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CHAPTER  I 


INTRODUCTION 

In  recent  years,  the  study  of  Reynolds  number  effects  has  become  an  increasingly 
important  area  of  research.  The  extrapolation  of  data  from  a  large  experimental  model  to 
a  full-scale  prototype  is  extremely  critical  in  predicting  the  performance  of  the  prototype. 
Often  the  large  model  is  still  much  smaller  than  the  prototype.  This  is  particularly  true  of 
naval  vessels,  including  submarines.  The  Office  of  Naval  Research  funded  a  Reynolds 
number  study  to  examine  the  capability  of  computational  fluid  dynamics  (CFD)  solutions 
to  track  large  changes  in  Reynolds  number.  Three  length  scales  of  submarine  hull  forms 
were  to  be  used  in  this  study.  A  small  model  was  to  be  constructed  and  tested  in  a  water 
tunnel  facility  located  at  Mississippi  State  University.  This  model  was  to  be  identical  in 
shape  to  a  model  tested  in  the  U.  S.  Navy's  Large  Cavitation  Channel  (LCC),  but  the  LCC 
model  was  to  be  over  twenty  times  larger.  The  flow  conditions  of  speed  and  pressure 
were  to  be  the  same  in  both  experiments.  CFD  solutions  were  to  be  obtained  for  the 
small  and  large  models,  as  well  as  for  the  full-scale  submarine,  which  is  over  twenty 
times  larger  than  the  LCC  model.  It  was  hoped  that  the  information  obtained  concerning 
how  well  the  CFD  solutions  tracked  the  changes  in  Reynolds  number  when  going  from 
small  model  to  large  model  would  be  useful  in  predicting  how  well  the  CFD  solutions 
would  track  the  next  jump  in  order  of  Reynolds  number  magnitude  when  going  from  the 
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large  LCC  model  to  the  full-scale  submarine.  The  objective  of  the  current  experiment 
was  to  obtain  drag  coefficient  measurements  of  the  smaller  model  submarine. 

The  Mississippi  State  University  Department  of  Aerospace  Engineering  has 
performed  hydrodynamic  research  in  conjunction  with  the  United  States  Navy  for  several 
years.  The  water  tunnel  in  the  Mississippi  State  Patterson  Engineering  Laboratories  has 
been  renovated  to  perform  research  work  with  the  LCC.  The  first  phase  of  research 
studied  the  effect  of  pressure  tap  hole  size  on  static  pressure  measurements.  After  this 
first  phase  of  testing,  the  water  tunnel  needed  many  repairs.  The  protective  coating  was 
replaced  in  many  sections  of  the  water  tunnel.  Also,  the  modification  of  the  interior  drive 
shaft  bearing  assembly  and  the  addition  of  a  reverse  thrust  bearing  solved  a  problem  that 
had  developed  with  the  bearings  being  burned  out  during  operation.  These  modifications 
were  completed  before  the  current  experiment  was  undertaken. 

The  Aerospace  Engineering  Department  acquired  a  Laser  Doppler  Velocimeter 
(LDV)  for  use  with  the  water  tunnel,  among  other  things.  A  calibration  device  consisting 
of  a  precisely  controlled  rotating  disk  was  built  for  the  LDV.  Once  the  water  tunnel  was 
operational,  the  LDV  was  used  to  measure  the  flow  characteristics  of  the  water  tunnel  test 
section.  After  these  tests  were  completed,  a  submarine  model  was  mounted  in  the  test 
section  and  wake  velocity  profiles  were  measured  at  several  stations  downstream.  The 
drag  coefficient  of  the  submarine  model  was  then  computed  using  these  wake  velocity 
profiles.  The  calculated  drag  coefficients  of  the  model  submarine  were  in  fairly  good 
agreement  with  the  CFD  solutions.  During  the  final  stages  of  the  experiment,  a  section  of 
the  water  tunnel  failed  under  pressure,  effectively  ending  the  experiment. 
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CHAPTER  II 


BACKGROUND 


2.1  Total-Head  Loss  Measurement  of  the  Drag  of  a  Body 

The  classical  approach  to  the  drag  measurement  of  a  body  has  been  to  use  a  pitot 
probe  to  measure  the  total  pressure  distribution  of  the  wake  behind  an  experimental 
model.  Two  of  the  earliest  methods  were  derived  by  Betz  and  by  Jones,  as  discussed  in 
Modem  Developments  in  Fluid  Dynamics  Betz  was  the  first  to  develop  the  technique 
of  determining  the  drag  of  an  airplane  wing  by  measuring  the  pressure  distribution  of  the 
wake  at  a  short  distance  behind  the  wing’s  trailing  edge.  From  the  momentum  equation, 
the  drag  of  a  wing  can  be  calculated  using  the  flow  properties  through  two  parallel  planes 
at  right  angles  to  the  mean  flow: 

D=  JJ| (Pl+pU?)dSt  -  JJ (p2  +  pU22) dS2 ,  (2.1) 


where  the  subscripts  1  and  2  denote  the  front  plane  and  rear  plane  respectively.  Betz 
introduced  the  total-head  equations  and  a  hypothetical  flow  to  transform  the  above 
equation  to  an  integral  requiring  measurements  on  the  rear  plane  only: 


D=  H{H]-H2)dS  +  ^P\\(ut2-U2)(u;+U2-2U0)dS  + 

T  ^  T 

+  \pMV?+Wi)-W+Wl)}dS- 


(2.2) 
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where 


(2.3) 


and 


u;=. 


2  (Hx-p2) 


(2.4) 


and  the  region  T  is  the  region  of  the  wake  only.  The  sum  of  the  first  two  integrals  of 
equation  (2.2)  yields  the  profile  drag,  and  the  third  integral  yields  the  induced  drag.  The 
second  term  in  equation  (2.2)  is  a  correction  factor  that  is  important  in  regions  close  to 
the  body.  It  is  small  when  the  measurement  plane  is  at  a  considerable  distance  from  the 
body.  If  there  is  no  lift  being  produced  by  the  body  and  the  measurement  plane  is  at  a 
certain  distance  downstream  of  the  body,  the  profile  drag  is  simply  the  integration  of  the 
head  deficit  of  the  wake: 

D=  \\{H,~H2)dS.  (2.5) 

T 

B.  M.  Jones  developed  a  similar  method  using  measurements  in  a  plane  far 
downstream  of  a  body.  For  a  body  at  zero-lift  conditions,  at  a  certain  distance 
downstream  the  pressure  is  uniformly  equal  to  the  freestream  pressure  p„  and  the 
velocity  is  uniformly  equal  to  the  freestream  velocity  U„ ,  except  in  the  wake.  The 
profile  drag  of  the  body  is  simply  an  integration  over  the  wake  region: 

D  =  pjjui{U„-U3)dS3,  (2.6) 
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where  U3  is  the  velocity  in  area  element  dS^  of  the  wake.  The  total  head  at  this 
downstream  plane  is 

H2  =  '-pU]+p..  (2.7) 

Defining  the  function  g  such  that 

r, 

8=  1 - ;  rj2  ’ 

L  2PU~  J 

equation  (2.6)  reduces  to  the  form 

D  =  \pUl\\2gm(\-g'n)dSs.  (2.9) 

At  a  certain  distance  downstream  of  the  body,  H2  — >  and  g  — >  1 ,  and  the  above 

equation  (2.9)  can  be  written  approximately 

D=X2PUl\\{\-g)dS,  =  \\{H„-H2)dS,  ,  (2.10) 

which  is  the  same  form  as  equation  (2.5). 

2.2  Momentum  Defect  Measurement  of  the  Drag  of  an  Axisymmetric  Body 

2.2.1  Derivation  of  Drag  Coefficient 

In  the  current  experiments,  it  was  not  possible  to  make  pressure  measurements  in 
the  wake  of  the  test  model.  Only  velocity  measurements  were  possible  with  the  LDV 
system,  which  makes  local  velocity  measurements  at  the  intersection  of  two  crossed  laser 
beams.  Because  the  test  section  of  the  water  tunnel  was  cylindrical,  the  LDV  was  only 
able  to  measure  points  lying  on  the  horizontal  transverse  plane  that  bisected  the  test 
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section  through  the  sides.  This  limitation  was  imposed  by  the  requirement  that  the  each 
laser  beam  must  have  the  same  incidence  angle  when  entering  the  cylindrical  test  section. 
Since  a  full  velocity  field  cannot  be  measured  with  the  cylindrical  test  section,  the  flow 
was  necessarily  assumed  to  be  axially  symmetric.  Because  of  these  limitations,  it  was 
necessary  to  derive  a  formula  for  the  drag  coefficient  of  an  axisymmetric  body  that  only 
made  use  of  velocity  measurements  to  compute  the  drag  coefficient.  The  following 
derivation  closely  follows  the  procedure  used  by  Dimotakis  2  for  a  Cartesian  coordinate 
system  and  a  nominally  two-dimensional  flow. 


Figure  2.1 :  Control  volume  representation. 

Consider  an  axisymmetric  body  located  along  the  x-axis.  The  body  is  surrounded 
by  a  cylindrical  volume  ■¥-  along  the  x-axis  whose  surfaces  are  ABCDE,  as  shown  in 
Figure  2.1.  Surface  A  is  the  transverse  plane  upstream  of  the  submerged  body,  where 
uniform  freestream  conditions  are  assumed  to  apply.  Surface  B  is  a  streamtube 
surrounding  the  submerged  body;  surface  C  is  the  transverse  plane  downstream  of  the 
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body;  surface  D  is  the  fluid  surface  adjacent  to  the  body;  and  surface  E  is  a  cylindrical 
surface  of  infinitesimally  small  diameter  extending  along  the  x-axis  from  surface  D  to 
surface  C. 

The  Navier-Stokes  equations  are 

div  u  =  0,  (2.11) 

and 

^H+V«»+-(Vg-Vr)  =  0,  (2.12) 

where  ii  is  the  velocity  vector,  p  is  the  fluid  density,/?  is  the  pressure,  and  T  is  the 

viscous  stress  tensor.  Integrating  the  momentum  equation  with  respect  to  the  interior 
volume  V-  and  applying  Gauss’  theorem,  this  becomes 

y  \udT+  \uu  dS  +  -\{p-p„)dS--\TzdS  =  Q,  (2.13) 

dt  ^  $  P  5  P  $ 

where  dS  is  a  vector  directed  along  the  outward  surface  normal,  of  magnitude  equal  to 
the  surface  element  dS;  and  the  freestream  pressure  p„  is  chosen  as  the  origin  in 
pressure.  The  various  operators  and  vectors  in  equation  (2.13)  are  defined  in  cylindrical 
coordinates: 


Gradient  operator 

„  d  ee  d  „  d 
v  =  e'&+73 ~e+e’Jr 

(2.14) 

Velocity  vector 

u  =  exu  +  egv  +  erw 

(2.15) 

Surface  element  vector 

dS  =  exdSx  +  egdSg  +  erdSr 

(2.16) 
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The  surface  element  magnitudes  are  further  defined: 

dSx  -  rdddr ; 

dSg  =  drdx ; 
dSr  =  rdOdx . 

The  viscous  stress  tensor  T  is  defined  for  incompressible  flow  as 

t  =  h\Vu  +  (Vu)T  , 


(2.17) 

(2.18) 
(2.19) 


(2.20) 


where  fl  is  the  absolute  viscosity.  When  T  is  expanded  in  cylindrical  coordinates,  it  is 
written  as: 


dw 


1  dw  dv\ „  ( dw  du 


t  =  H\2 -erer+{--  +  -)eree+^-+dr^ 


e  e  + 

r  x 


f  dv  1  dw 
ydr  rdd 


eger  + 


(2  dv  ' 

f dv  1  du ) 

QJ 

S 

OJ 

+ 

Krd0) 

|  dgSg  4-  j 

Kdx  rddj 

hK  + 

[dr  dx  J 

exer  + 


(2.21) 


+ 


1  du  dv'l  „  „  _  du  „  „ 


r  dO  dx  J  dx 


When  the  cylindrical  coordinate  operators  and  vectors  are  substituted  into 
equation  (2.13)  and  dotted  with  ex,  the  x-component  of  the  momentum  equation  is 
obtained: 


—  judV-  +  ^uwdSr  +  juvdS0  +  juudSx  +  —  j(p  —  p„)  dSx  + 
dt  v  $  $  $  P  s 


■If. 


du  dw 

i\Yr  +  dx 


V  rf  1  du  dv  V 

dSr  -v  -  —  +  —  dS 
J  Ar  d&  dx  J 


=  0 


(2.22) 
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Taking  the  time  average  of  the  equation  eliminates  the  time  derivative  term: 


j ~uwdSr  +  jm’dSg  +  j \uudSx  +—  J(j v~  P~)dSx  —  vj 
s  P  s 


^  du 

i  dr  dx j 


tis  + 


(\  du  dv^ 


r  I  au  av  r 

—v - +  —  dS0-v\ 

•  f  rtf)  ^ 


sKrde  dxj 


r2du^ 

V  dxJ 


(2.23) 


ks  =0. 


At  this  point  it  is  helpful  to  decompose  each  surface  integral  to  the  sum  of  the 
integrals  over  each  sub-surface  A,  B,  C,  D,  E;  i.e., 

f=f+f+f+f+f.  (2.24) 

S  A  B  C  D  E 


The  integrals  over  surface  E  are  zero  because  surface  E  has  been  chosen  to  be  a 
streamtube  of  infinitesimally  small  diameter  across  which  there  is  no  flux.  The  integral 
over  surface  D  represents  the  drag  force  per  unit  mass  per  unit  volume  that  the  fluid 
control  volume  ¥•  exerts  on  the  body  in  the  x-direction.  The  resulting  equation  is: 


—  =  -j ~uwRdddx-\  J+  yru2rd6dr-  —  \  J+  j| (p  -  pS) rdOdr  + 

P  b  cj  P  c J  ^2  25) 

+yl^+^}dMx+2vlj+j}sirMdr  • 

where  R  is  the  radius  of  the  cylindrical  control  volume.  If  the  maximum  body  diameter  d 
is  the  characteristic  length  of  the  submerged  body,  and  equation  (2.25)  is  divided  by 

—  Ul—d1,  then: 

2  4 
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- - -  = — —  [uwRdddx  ■ 

2  PUt  f  d2  xd2Ull 


Kd2Ul 


J+  y?u2rd6dr  + 


~Tdkr{\+i\(-p-p~)rMdr* 

8v  /  du  \6v 

nd2Ul  J  dr  dx  )  Kd2Ul 


(2.26) 


i+  l\^rd0dr- 


The  drag  coefficient  of  the  body  is  defined  as 


.pUljd2 


(2.27) 


and  the  time  averaged  pressure  coefficient  c  is  defined  as 


-  _  P~P- 
p  hpul 


(2.28) 


Applying  these  definitions,  equation  (2.26)  becomes: 


C°  TCd1 


8  Hz*'**  -  s?  {/+ 1)  mrdedr  ll  'rrMdr + 


8v  ff  du  dw  ,  \6v 

nd  Ul  il  dr  dx  j  nd  Ul 


16v  r  r \  du 


(2.29) 


h  i\^rd0dr  ■ 


Equation  (2.1 1)  is  integrated  over  the  control  volume  ¥- ,  and  Gauss’  Theorem  is 


applied  to  obtain: 


jVudf^=  judS  =  0. 


(2.30) 


Expanding  equation  (2.30)  and  multiplying  by  a  factor  of  — - - ,  the  continuity 

ndU _ 


equation  is  written  in  terms  of  integrals  over  the  sub-surfaces  A,  B,  and  C: 
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(2.31) 


8 


|+  ^iurdOdr  4 - - - jwRdddx  =  0 


7Cd2U 


U  c 


7td  U„ 


Equation  (2.31)  can  be  subtracted  from  the  Co  equation  (2.29): 


Cn  =■ 


8  ruw  -Uw 


ndl 


J 


ul 


RdOdx  - 


8 


nd 2 
81/ 


H 


M2  -t/_« 


..j  C 

A 


Ul 


rd8dr  + 


4-- 


7td 2  ( 

16r 

;r</2t/2 


[4-  \ — r  f 

J  p  ^2f/2  J 


A  C  J 


W 


du  d 
dr  dx 


IRdddx  + 


(2.32) 


Surface  A  has  been  chosen  to  be  far  enough  upstream  of  the  body  such  that  the 
velocity  and  pressure  are  assumed  to  be  uniform  across  the  surface  and  equal  to  the 


freestream  values  £/„  and  p„ ,  respectively.  Thus  the  second  and  third  terms  of  equation 
(2.32)  are  equal  to  zero  on  surface  A. 

The  velocity  components  u  and  w  can  be  separated  into  a  mean  value  and  a 
fluctuating  part,  i.e.: 


u  =  u  4-  u' ; 

(2.33) 

w  =  w  4-  w  ; 

(2.34) 

2  —  2  ,  n 

u  -u  +u  ; 

(2.35) 

uw  =  uw  +  u'w' . 

(2.36) 

The  product  uw  appears  only  in  the  integrals  over  surface  B.  Surface  B  has  been  defined 
to  be  a  streamtube  that  lies  outside  of  the  regions  where  disturbances  from  the  wall  or  the 
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body  are  assumed  negligible,  which  implies 

w  =  u'w'  =  ~~  =  0  (2.37) 

ox 

on  surface  B. 

Applying  the  conditions  on  surfaces  A  and  B,  and  assuming  that  the  surfaces  are 
far  enough  away  from  the  body  that  all  mean  velocity  gradients  are  negligible,  equation 
(2.32)  is  reduced  to: 


Now  define  the  displacement  thickness  8]  at  the  downstream  surface  C  such  that 


{nR1  -n8\^U0  =  j 'urdGdv, 

c 


(2.39) 


•sMjf 


~  \ 
1--0- 
V  UoJ 


rdOdr\ 


(2.40) 


where  U0  is  defined  as  the  velocity  outside  of  the  wake  at  surface  C.  Since  there  is  no 
mass  flux  through  surface  B,  the  flow  rate  through  surfaces  A  and  C  must  be  equal: 

ttR2U„  =  [ttR2  - k82  ) UQ .  (2.41) 


Equation  (2.41)  can  be  rearranged  to  the  form: 


(2.42) 
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The  effects  of  freestream  acceleration  due  to  displacement  thickness  Sx  on  the 
pressure  coefficient  can  be  separated  from  the  pressure  differences  across  the  wake  by 


dividing  the  pressure  coefficient  into  two  separate  components: 


7  =P-P~=Po-P»  ,  P-Po  -c  +c 

P  {pul  1 pul  {pul  p" 


(2-43) 


where  cp>  is  the  pressure  coefficient  of  the  flow  outside  of  the  wake,  and  cK  is  the  time- 

averaged  pressure  coefficient  of  the  flow  within  the  wake.  Bernoulli’s  equation  can  be 
applied  outside  of  the  wake  between  surfaces  A  and  C: 


Po+\pu2o=P~+\pU~- 


(2.44) 


Using  equation  (2.42)  with  equation  (2.44),  cpo  can  be  defined  in  terms  of  the 


displacement  thickness: 


So=- 


2- 


R 


1- 


(2.45) 


Using  equations  (2.45),  (2.40),  and  (2.42)  in  equation  (2.38),  the  equation  for  the  drag 
coefficient  becomes: 
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2.2.2  Self-Preserving  Flow 


The  generation  of  a  wake  behind  a  body  is  an  example  of  a  developing  flow.  The 
velocity  deficit  of  the  wake  is  greater  and  the  width  of  the  wake  is  smaller  near  the  body, 
but  as  the  flow  moves  downstream  the  velocity  deficit  decreases  and  the  wake  width 
increases  as  the  wake  turbulence  dies  out.  An  assumption  of  self-preserving  flow 
stipulates  that  the  transverse  velocity  profiles  retain  the  same  functional  forms  at  points 
further  downstream  of  the  body,  but  that  the  length  scales  and  mean  flow  property  values 
scale  according  to  the  distance  downstream  of  the  body.  According  to  Townsend  ,  the 
definition  of  a  self  preserving  flow  is  one  in  which  the  variation  of  a  mean  value  quantity 
Min  the  flow  direction  must  be  of  the  form: 

M  =  M,  +  m0(x)func(r / l0(x)) ,  (2-47) 

where  the  length  scale  /0  and  the  scale  of  quantity  m0  are  functions  of  x,  and  Mj  is  the 
reference  mean  value. 

Assume  u  ,  un  ,  and  p  have  the  forms: 

u(r)  =  U0  -  u0  (x)f(r  /  lg  (x))  ;  (2.48) 

u'2  =Uo(x)g(r/l0(x));  (2.49) 

p  =  pt)+~pu20(x)n(r/ljx)).  (2.50) 

Introduce  the  non-dimensional  coordinate  rj : 


(2.51) 
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Then, 


u(r)  =  U0-u0(x)f(rj), 


u  =u0(x)g(Tj),  and 


p  =  p0+^pul(x)n(rf)- 


(2.52) 


(2.53) 


(2.54) 


Applying  these  forms  to  equation  (2.40)  and  integrating  in  0  results  in  the  displacement 
thickness  in  terms  of  the  length  and  velocity  scales,  and  an  integration  of  the  function 

f{v)- 

;2 /  \  R/k 
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(2.55) 


Now  consider  the  time  averaged  pressure  coefficient  of  the  wake  c  : 
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Substituting  the  equations  (2.42)  and  (2.54)  into  equation  (2.56)  yields 


(2.56) 


(2.57) 


When  the  functional  forms  of  equations  (2.51)  through  (2.55),  and  equation  (2.57)  are 
applied  to  the  drag  coefficient  equation  (2.46),  the  resulting  equation  is: 
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where  uo  and  lo  are  implicitly  functions  ofx. 

White  4  states  that  the  velocity  defect  of  a  circular  wake  can  be  expressed  as  a 
Gaussian  distribution.  The  function  / ( TJ )  is  assumed  to  have  the  Gaussian  form: 


=  (2.59) 

The  presence  of  the  (In  2)  factor  in  the  exponent  of  equation  (2.59)  references  the  length 
scale  to  the  wake  half  width  b,  which  is  defined  as  the  width  of  the  wake  where  the 
velocity  deficit  is  half  of  the  maximum  velocity  deficit.  Substitution  for  / (rj)  in 


(2.60) 

According  to  White  4,  the  self-similarity  parameters  ug  and  lo  of  a  circular  wake 

_2/ 

with  small  velocity  defect  are  described  by  the  power  law  relations  u0  x  73  and 
/„  OC  .  A  further  fit  is  introduced  such  that 


equation  (2.58)  and  integration  results  in: 
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where  xo  is  the  origin  from  the  curve  fitting  process.  The  drag  coefficient  equation 
becomes: 
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(2.63) 


By  extrapolating  far  enough  downstream,  the  terms  containing 


the  equation  for  drag  coefficient  becomes: 
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If  it  is  assumed  that  the  body  diameter  d  is  much  less  than  the  control  volume  radius  R, 
the  terms  containing  (d/R)2  may  be  neglected,  and  the  equation  for  CD  reduces  to 
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CHAPTER  III 


EXPERIMENTAL  APPARATUS 

3.1  Water  Tunnel 

3.1.1  Water  Tunnel  History  and  Description 
The  Aerospace  Engineering  Department  Water  Tunnel  is  located  in  Patterson 
Engineering  Laboratories  on  the  Mississippi  State  University  main  campus.  A  layout 
drawing  of  the  water  tunnel  is  shown  in  Figure  3.1.  The  preliminary  design  of  the  water 
tunnel  was  begun  in  1961  by  Dr.  Graham  Wells  5.  The  tunnel  was  constructed  from 
fiberglass  by  Aerospace  Engineering  Department  shop  technicians,  with  the  largest 
section  diameter  being  40  inches,  the  limit  of  shop  capabilities.  The  test  section  is  a  4-ft 
length  of  1 1-inch  inner  diameter  Plexiglas®  tube.  This  results  in  a  contraction  ratio  of 
13.22:1.  The  length  of  the  contraction  is  5  feet,  and  the  contraction  was  designed  for 
minimal  adverse  pressure  gradient  and  a  uniform  exit  velocity  profile.  The  original  drive 
system  was  a  Buick  V-8  automotive  engine  with  belt  drive,  but  this  was  later  changed  to 
an  electric  motor  retaining  the  old  belt  drive,  and  finally  to  an  efficient  direct-drive 
electric  motor.  The  drive  section  is  located  on  the  first  floor  of  Patterson  Laboratory,  and 
the  test  section  is  located  on  the  second  floor.  The  current  motor  has  a  remote  control 
panel  so  that  the  tunnel  can  be  operated  from  the  second  floor.  The  3-bladed  propeller 
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Figure  3.1 


was  reportedly  taken  from  a  Second  World  War  amphibious  assault  vehicle.  The  water 
fill  line  is  located  on  the  second  floor.  The  valve  to  this  water  line  is  opened  to  pressurize 
the  water  tunnel  to  the  line  pressure,  which  is  usually  about  50  psig.  The  tunnel 
pressurization  and  the  hydrostatic  head  suppress  cavitation  from  occurring  in  the  tunnel 
drive  section.  The  tunnel  contains  no  flow  straighteners,  but  does  have  turning  vanes  on 
each  of  the  comers.  Another  set  of  tunnel  segments  consisting  of  contraction,  test 
section,  and  diffuser  were  built  to  study  nominally  two-dimensional  flows,  but  these 
sections  were  never  fitted  into  the  tunnel  and  are  in  storage. 

The  primary  motivation  for  the  construction  of  the  water  tunnel  at  MSU  was  to 
attract  hydrodynamic  research  from  the  submarine  and  shipbuilding  industries,  which 
were  booming  in  the  1950’s  and  1960’s  as  the  modem  nuclear  submarine  was  being 
developed.  The  facility  never  saw  any  of  this  research,  and  for  most  of  the  1 980’s  and 
the  1990’s  the  tunnel  was  not  operated  nor  maintained.  It  was  renovated  in  the  late 
1990’s  to  support  research  at  the  new  U.S.  Navy  Large  Cavitation  Channel  (LCC)  in 
Memphis,  Tennessee.  The  first  experiments  performed  after  the  renovation  studied  the 
effect  of  tap  hole  size  on  static  pressure  measurements  in  the  test  section  6.  During  these 
experiments  a  number  of  problems  developed  with  the  tunnel  that  required  further  work, 
and  the  tunnel  was  taken  apart. 

3.1.2  Repainting  of  Comers,  Turning  Vanes,  and  Drive  Section 

At  the  end  of  the  previous  testing  on  the  water  tunnel,  it  was  noticed  that  certain 
things  had  to  be  repaired.  Although  the  entire  tunnel  interior  had  received  a  new  finish  of 
four  coats  of  Duratec®  vinyl  ester,  the  finish  on  the  tunnel  comers,  turning  vanes,  and 
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drive  section  had  peeled  off.  It  was  necessary  to  recoat  these  comer  sections  to  protect 
the  fiberglass  structure  from  water  penetration.  Additional  coats  of  Interlux 
InterProtect®  epoxy  paint  were  applied  to  the  tunnel  drive  section  and  the  third  comer 
downstream  of  the  test  section.  InterProtect®  is  formulated  to  prevent  water  penetration 
into  fiberglass  boat  hulls. 

3.1.3  Drive  Bearing  Modifications 

A  more  serious  problem  that  developed  during  the  previous  testing  of  the  tunnel 
was  bearing  overheating  and  burnout.  The  water-cooled  Cutless  bearing  inside  the  drive 
section  of  the  tunnel  was  completely  burned  out.  When  the  drive  section  was 
disassembled,  it  was  discovered  that  the  step  of  the  drive  shaft  had  been  cut  into  by  the 
edge  of  the  brass  Cutless  bearing.  It  was  determined  that  the  location  of  the  shaft  had 
shifted  downstream  relative  to  the  support  bearings,  causing  the  step  to  butt  against  the 
brass  Cutless  bearing.  This  not  only  cut  off  the  water  flow  through  the  bearing,  but 
damaged  the  stainless  steel  shaft  as  well.  To  ensure  adequate  water  flow,  a  Little  Giant® 
1/40  hp  water  pump  was  placed  into  the  bearing  cooling  water  feed  line.  Previously, 
water  was  forced  through  the  feed  line  by  the  differential  pressure  between  the  third 
comer  (where  the  water  is  picked  up)  and  the  drive  section. 

It  was  also  noted  that  a  hydraulic  press  was  needed  to  press  the  Cutless  bearing 
into  and  out  of  its  stainless  steel  holder.  Measurement  of  the  holder  inner  diameter  and 
the  Cutless  bearing  outer  diameter  showed  that  the  holder  diameter  was  0.015  inch  too 
small,  and  out  of  round.  It  was  speculated  that  local  buckling  of  the  bearing  as  it  was 
pressed  into  the  holder  caused  the  inner  diameter  to  shrink  enough  to  cause  the  bearing  to 
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grip  the  drive  shaft  too  tightly.  A  new  bearing  slipped  onto  the  shaft  easily,  but  after 
being  pressed  into  the  bearing  holder  it  could  only  be  forced  onto  the  shaft.  It  was 
decided  to  bore  out  the  bearing  holder  to  a  diameter  that  was  0.002  inch  less  than  the 
bearing  outer  diameter.  The  parts  were  prepared  for  installation  by  placing  the  bearing  in 
a  freezer  and  by  heating  the  bearing  holder  in  an  oven.  The  chilled  bearing  was  quickly 
inserted  into  the  heated  holder.  Afterwards,  as  the  temperatures  equalized,  the  bearing 
expanded  and  the  holder  shrinked,  providing  a  solid  grip  on  the  bearing. 

Other  modifications  were  made  to  the  bearing  to  ensure  adequate  water  flow. 

First,  the  bearing  length  was  reduced  by  0.4  inch  to  fit  completely  within  the  holder 
sleeve.  Previously,  the  bearing  had  protruded  from  the  holder  end  and  had  severely 
restricted  the  flow  into  the  holder  from  the  feed  line  outlets  on  either  side  of  the  holder 
end.  Another  minor  modification  to  the  bearing  was  the  filing  of  ramps  into  the  brass  and 
rubber,  to  aid  the  water  flow  into  the  rubber/shaft  interface.  Furthermore,  the  drive  shaft 
steps  in  this  region  were  moved  back  0.875  inch  to  insure  no  further  contact  with  the 
bearing  end. 

Later,  during  tunnel  assembly,  it  was  finally  determined  why  the  drive  shaft  had 
drifted  downstream  and  caused  the  Cutless  bearings  to  fail.  During  the  last  renovation 
and  upgrade  of  the  tunnel,  a  new  50  hp  electric  motor  was  installed.  Power  was  directly 
transmitted  to  the  shaft  by  means  of  a  flexible  connection.  It  was  discovered  during 
tunnel  assembly  that  the  previous  power  train  lacked  a  reverse  thrust  bearing.  The  power 
system  before  that  had  a  belt  drive  system  that  incorporated  a  reverse  thrust  bearing.  A 
thrust  bearing  and  locking  collar  were  fitted  onto  the  shaft  so  that  the  new  thrust  bearing 
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worked  against  the  existing  outer  bearing  casing.  This  prevented  the  shaft  from  shifting 
longitudinally  when  reverse  thrust  was  applied.  Although  the  motor  was  not  run  in 
reverse,  reverse  thrust  was  applied  when  the  tunnel  was  slowed  or  stopped,  causing  the 
propeller  to  work  against  the  water  flow.  It  was  speculated  that  in  the  previous  power 
train,  the  shaft  was  pulled  downstream  just  enough  to  bind  the  shaft  step  against  the 
Cutless  bearing  during  these  transitional  phases  of  operation. 

3.1.4  Other  Tunnel  Modifications 

Several  other  modifications  were  made  to  the  tunnel  to  improve  tunnel  installation 
and  maintenance.  A  new  and  larger  drain  line  was  placed  into  the  bottom  of  the  third 
comer.  A  Culligan®  HF-150  water  filter  rated  at  125  psi  was  fitted  into  the  fill  line  of 
the  water  tunnel.  Culligan®  SI  A  20-micron  sediment  water  filter  cartridges  were  used  to 
filter  out  the  great  amounts  of  rust  and  sediment  in  the  fill  line,  providing  cleaner  water  in 
the  tunnel.  To  further  aid  in  keeping  the  tunnel  clean,  common  pool  algaecide  was  added 
to  the  tunnel  to  prevent  organic  pollution.  Ball  valves  were  added  onto  the  percolator 
domes,  which  are  installed  on  each  comer  of  the  water  tunnel  to  collect  air  that  is 
pressurized  out  of  solution  with  the  water  during  operation.  Algaecide  and  seeding 
particles  can  be  added  to  the  filled  water  tunnel  by  use  of  the  funnel  on  the  percolator  on 
the  first  comer. 

The  bearing  bolts  that  hold  the  bearing  holder  in  the  drive  section  of  the  tunnel 
were  modified.  Previously,  these  bolts  were  very  short  and  difficult  to  reach.  Also,  the 
bolt  had  to  pass  through  three  separately  threaded  lugs.  All  but  the  bearing  holder  lug 
holes  were  stripped  of  their  threads.  5/16x3-1/8  inch  Allen  head  bolts  were  then  used  to 
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secure  the  holder.  The  extra  length  of  these  new  bolts  extended  the  bolt  head  to  the 
outside  of  the  drive  fairing,  which  greatly  simplified  the  assembly  of  the  tunnel  drive 
section.  These  bolts  had  to  be  custom  made  from  3/8  x  3-1/8  inch  bolts.  The  ends  of 
these  bolts  were  turned  down  on  a  lathe  to  5/ 16-inch  diameter  and  threaded. 

3.1.5  Propeller  Knocking 

After  final  assembly  of  the  tunnel  and  initial  run-up,  a  loud  knocking  noise  was 
heard  in  the  propeller  area.  Because  the  clearance  between  the  propeller  and  the  wall  is 
less  than  a  paper  thickness,  the  most  likely  source  of  the  noise  was  the  propeller  striking 
the  wall.  After  running  the  tunnel  unpressurized  for  some  time,  the  noise  was  reduced  at 
the  higher  normal  operating  pressures.  It  was  speculated  that  the  propeller  had  worn 
itself  in  during  the  unpressurized  run.  The  reduced  noise  was  deemed  acceptable,  and 
was  ignored  during  further  testing. 

3.2  Laser  Doppler  Velocimeter 

The  laser  Doppler  velocimeter  (LDV),  which  is  also  called  a  laser  Doppler 
anemometer  (LDA),  employs  the  Doppler  theory  to  measure  particle  velocities  at  a  point 
in  a  fluid  flow.  The  measurement  probe  of  the  system  is  a  lens  that  focuses  two  laser 
beams  of  the  same  wavelength  onto  a  point  determined  by  the  focal  length  of  the  probe 
lens.  The  intersection  of  the  two  beams  forms  an  ellipsoid-shaped  measurement  volume 
where  a  distinct  pattern  of  light  is  formed.  The  amplitudes  of  the  two  laser  beams  add 
together  to  form  alternating  regions  of  high  amplitude  and  low  amplitude  light — this  is 
called  the  fringe  pattern.  If  a  light-reflecting  particle  were  to  pass  through  this 
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measurement  volume  in  a  direction  normal  to  the  fringe  pattern,  it  would  reflect  pulses  of 
light  as  it  passes  through  the  high  amplitude  fringes  of  the  measurement  volume.  The 
wavelength  and  the  angle  of  intersection  of  the  beams  are  known,  which  determines  the 
distance  between  each  fringe  in  the  measurement  volume.  The  LDV  system  thus 
measures  the  time  interval  between  reflected  bursts  of  light,  and  is  able  to  calculate  the 
particles  speed  through  the  measurement  volume  in  the  direction  normal  to  the  fringe 
pattern. 

The  above  description  assumes  that  the  frequencies  of  the  two  beams  are  the 
same,  and  are  in  “sync”.  This  is  achieved  by  optically  splitting  a  single  laser  beam. 
However,  it  is  undesirable  to  have  both  beams  of  the  same  frequency.  If  the  frequency  of 
one  of  the  beams  is  shifted  by  a  slight  amount,  the  fringe  pattern  “scrolls”  in  one 
direction.  This  allows  the  LDV  system  to  discern  between  negative  and  positive  velocity 
directions.  The  LDV  system  does  this  through  the  use  of  a  Bragg  cell.  The  Bragg  cell  is 
a  device  that  shifts  one  of  the  laser  beams  frequencies  by  directing  the  beam  through  a 
prism  that  is  electro-mechanically  vibrated  at  the  desired  frequency  shift. 

The  laboratory  LDV  is  a  Dantec  FiberFlow  system.  The  laser  source  is  a  Spectra 
Physics  Stabilite  2017  Argon-ion  laser.  This  unit  consists  of  a  laser,  a  laser  power 
supply,  and  a  remote  control  module.  The  laser  and  the  power  supply  run  on  220  V  and 
require  water-cooling.  The  output  beam  of  this  system  is  directed  into  the  FiberFlow 
transmitter,  which  contains  the  beam  splitter  and  Bragg  cell.  The  beam  splitter  splits  the 
primary  laser  beam  into  up  to  three  pairs  of  beams.  The  three  pairs  of  beams  may  be 
476.5  nm  wavelength  (violet),  488  nm  (blue),  or  514.5  nm  (green).  One  pair  of  beams 
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represents  one  velocity  component  that  can  be  measured.  Currently,  the  optics  are  not 
installed  on  the  violet  beams,  rendering  the  system  capable  of  only  two-component 
measurements.  The  Bragg  cell  puts  a  40-Mhz  frequency  shift  into  one  of  the  beams  of 
each  pair.  The  laser  beams  are  then  transmitted  to  a  2-D  probe  by  means  of  fiber  optic 
cables.  The  probe  has  a  80  mm  diameter  lens  with  a  focal  length  of  800  mm.  The  probe 
is  equipped  with  a  1 .95  ratio  beam  expander.  The  system  is  computer-controlled  and  the 
data  are  collected  and  processed  using  the  Dantec  Flow  Velocity  Analyzer  (FVA)  and  its 
associated  software. 


3.3  LDV  Calibration  Apparatus 

3.3.1  Calibration  Apparatus  Description 
The  LDV  was  calibrated  by  measuring  the  tangential  velocity  at  a  point  of  radius 
r  on  a  disk  with  a  rough  surface  rotating  at  a  precisely  controlled  angular  velocity.  The 
disk  was  rough-cut  from  a  sheet  of  0. 25-inch  Plexiglas  and  was  turned  down  on  a  lathe  to 
a  radius  of  12  inches.  An  aluminum  hub  was  attached  to  the  backside  of  the  disk  by  two 
screws.  The  disk  assembly  was  not  balanced,  nor  was  the  disk  perfectly  flat;  the  disk 
wobbled  ±0.381  mm  perpendicular  to  the  disk  plane  near  the  extreme  radius.  The  wobble 
should  not  have  affected  the  accuracy  of  the  calibration  measurements,  but  rather  it  just 
moved  the  disk  surface  back  and  forth  in  the  axial  direction  through  the  measurement 
volume,  which  had  a  dimension  of  4.303  mm.  A  rough  disk  surface  was  created  by 
attaching  a  12-inch  diameter,  80-grit  sandpaper  disk  with  a  spray-adhesive  such  as  3M 
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Super77®.  This  rough  surface  simulated  a  particle  flow  for  the  LDV  system  to  measure. 
A  pinhole  was  located  in  the  center  of  the  disk  for  laser  alignment  purposes. 

The  calibration  disk  was  rotated  by  a  Compumotor  SM233BE-NTQN  type  servo 
motor.  The  motor  was  mounted  onto  aluminum  T-stock,  which  was  then  clamped  onto 
the  lower  water  tunnel  test  section  tie  bars  with  four  “c”-clamps.  These  tie  bars  were 
located  just  beneath  the  test  section.  This  location  allowed  an  LDV  calibration  to  be 
performed  without  physically  moving  the  traverse  system  out  of  alignment  with  the  water 
tunnel.  It  was  necessary  to  shim  the  motor  mount  2.25  inches  above  the  lower  support 
bars  in  order  to  place  the  center  of  the  disk  within  the  vertical  limits  of  the  LDV  traverse. 
The  TQ10X  Servo  Controller  was  mounted  on  an  aluminum  beam,  which  acted  as  a  heat 
sink  and  allowed  the  unit  to  hang  from  the  lower  tie  bars. 

3.3.2  Uncertainty  of  Disk  Radius  Measurement 

The  first  task  in  the  calibration  process  was  to  locate  the  LDV  measurement 
volume  on  the  center  of  the  calibration  disk.  The  current  laboratory  arrangement  did  not 
allow  the  LDV  measurement  volume  to  measure  more  than  a  few  centimeters  below  the 
calibration  disk  center,  so  an  alignment  procedure  was  developed  using  only  the  top  half 
of  the  disk  surface.  First,  the  LDV  measurement  volume  was  visually  centered  by 
positioning  the  probe  so  that  the  measurement  volume  “disappeared”  into  the  pinhole  on 
the  calibration  disk  center.  The  centering  accuracy  of  this  visual  method  was  only  ±1 .0 
mm  A  more  accurate  centering  alignment  was  achieved  by  using  the  LDV  beam  system 
and  the  traverse  to  measure  the  disk  radius,  and  to  center  the  traverse  by  measuring  off  of 
the  top  and  sides  of  the  calibration  disk. 


-27- 


Designate  the  left  side  of  the  disk  as  x/,  and  the  right  side  as  X2-  The  radius  R  of 
the  disk  is  then: 


(3.1) 


After  the  LDV  measurement  volume  was  visually  centered  as  closely  as  possible,  the 


blue  beams  were  used  to  measure  the  radius  of  the  disk  at  xj  and  X2. 

According  to  the  methods  advocated  by  Coleman  and  Steele  7,  the  uncertainty  of 
the  disk  radius  is  given  by  the  following  equation: 


B]  +2 

*2 


dR)  dR 


dx,  Jl  dx. 


(3-2) 


where  P  denotes  the  random  uncertainty  associated  with  xk,  BXk  denotes  the  systematic 


uncertainty  associated  with  x*,  and  B  denotes  the  correlated  systematic  uncertainties  of 


the  variables  x/  and  X2-  Because  the  same  procedure  was  used  to  measure  both  x/  and  X2, 
the  same  random  and  systematic  uncertainties  were  applied  to  each  x-coordinate. 

Only  one  source  of  random  uncertainty  was  found  in  the  procedure,  and  that  was 
the  judgment  of  the  operator  as  to  the  exact  alignment  of  the  beams  onto  the  edge  of  the 
disk.  The  possible  error  in  this  alignment  was  estimated  to  be  half  of  a  beam  diameter,  or 
±0.675  mm.  Possible  random  uncertainty  of  the  traverse  was  considered  negligible,  as 
the  traverse  specification  reported  a  resolution  of  6.25  pm . 


There  were  a  number  of  systematic  uncertainty  sources.  One  systematic 
uncertainty  source  was  the  specified  accuracy  of  the  traverse,  which  was  ±0.2  mm. 
Another  systematic  uncertainty  was  in  the  radius  of  the  Plexiglas®  disk  itself.  Using  a 
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dial  indicator,  a  variation  in  disk  radius  of  ±0.3429  mm  was  measured.  The  total 
systematic  uncertainty  of  the  location  of  xj  and  X2  then  was  ±0.3967  mm. 

Because  the  measurement  of  xj  and  jq  was  by  the  same  method,  the  systematic 
uncertainties  of  these  two  variables  were  completely  correlated.  The  correlated 
systematic  error  term  is 

Bxx  =BXBX  .  (3.3) 

Using  these  uncertainty  estimates,  the  uncertainty  in  the  measurement  of  disk  radius  R 

was  ±0.4773  mm.  The  actual  measured  value  of  the  calibration  disk  radius  using  this 

procedure  was  152.9875  mm,  or  6.023  inches. 

3.3.3  Uncertainty  in  Radial  Positioning  During  Calibration 
All  random  and  systematic  uncertainties  of  the  radial  measurement  used  in  the 
calibration  procedure  were  associated  with  the  centering  of  the  measurement  volume  in 
the  y-direction  (vertical).  This  was  done  by  visually  locating  the  green  laser  beams  on 
the  top  edge  of  the  calibration  disk,  which  introduced  a  random  error  of  ±0.675  mm,  as 
discussed  previously.  The  measurement  volume  was  then  moved  down  a  distance  equal 
to  the  disk  radius,  which  introduced  a  systematic  error  of  ±0.4773  mm.  Accounting  for 
the  additional  systematic  error  introduced  by  the  traverse  accuracy  (±0.2  mm)  and  the 
finite  diameter  of  the  measurement  volume  (±0.0995  mm),  the  total  systematic 
uncertainty  in  the  radial  position  measurement  was  ±0.5270  mm.  Because  the  traverse 
was  able  to  resolve  to  within  6.25  pm  of  a  given  position,  and  the  step  sizes  during  the 
calibration  tests  were  on  the  order  of  millimeters,  the  total  random  uncertainty  in  the 
radial  position  measurement  was  ±0.675  mm. 
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3.3.4  Positioning  the  Measurement  V olume  in  z 


Before  proceeding  with  the  calibration,  the  measurement  volume  must  be  located 
such  that  its  center  is  bisected  by  the  surface  of  the  calibration  disk.  After  the 
measurement  volume  was  centered  in  x  and  y,  the  volume  was  moved  to  a  point  close  to 
the  center  of  the  disk,  but  not  in  the  pinhole.  With  the  calibration  disk  rotating,  the  z 
position  was  adjusted  to  maximize  the  LDV  data  rate.  This  was  done  by  enabling  the 
high  voltage  in  the  FVA  software  and  connecting  an  oscilloscope  to  the  appropriate 
channel  on  the  FVA.  The  z  position  was  then  adjusted  to  maximize  the  signal-to-noise 
ratio  as  was  seen  on  the  oscilloscope.  To  further  improve  the  calibration  data,  the  high 
voltage  settings  were  adjusted  to  increase  the  signal-to-noise  ratio,  again  using  the 
oscilloscope  as  a  guide.  It  was  found  that  a  high  voltage  setting  of  800-900  worked  well 
for  this  purpose. 

3.3.5  Uncertainty  in  Angular  Velocity  Measurement 
The  angular  velocity  of  the  servo  motor  was  precisely  controlled  and  reported  to 
the  computer  control  software  every  few  seconds.  A  stroboscopic  tachometer  was  used 
to  verify  the  accuracy  of  the  rotational  velocity  at  low  and  high  speeds.  There  was  no 
discrepancy  between  the  angular  velocity  reported  by  the  motor  controller  and  the 
velocity  measured  by  the  tachometer.  The  systematic  uncertainty  of  the  rotational 
velocity  was  set  at  the  resolution  limits  of  the  stroboscopic  tachometer,  which  was  ±0.05 
rad/s. 
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During  the  actual  calibration  measurements,  data  were  only  taken  when  the  motor 
controller  reported  a  stable  reading  of  the  target  speed.  However,  at  some  speeds 
(especially  lower  speeds)  a  random  fluctuation  of  disk  speed  was  noticed  in  the  digital 
readout  of  the  motor  control  program.  This  flicker  was  on  the  order  of  0.03  rev/s.  Thus 
the  random  uncertainty  for  the  disk  speed  was  estimated  to  be  ±  0.03  rev/s,  or  ±0.18 
rad/s. 


3.3.6  Random  Uncertainty  of  LDV  V elocity  Measurement 
The  random  uncertainty  of  the  mean  velocities  measured  by  the  LDV  was 
determined  to  be  ±0.004  m/s.  This  estimate  was  developed  by  taking  25  measurements 
of  a  reference  speed,  and  determining  the  standard  deviation  of  the  mean  velocities  of 
these  measurements.  The  random  uncertainty  of  an  LDV  mean  velocity  measurement  is 
twice  the  standard  deviation.  These  measurements  were  done  for  several  reference 
velocities,  as  shown  in  Figure  3.2. 
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Figure  3.2:  Random  uncertainties  of  LDV  mean  velocity  measurements. 

3.3.7  Calibration  Results 

Several  calibration  tests  were  made.  In  one  type  of  test,  the  measurement  volume 
was  placed  at  a  specific  radius  and  the  disk  angular  velocity  was  stepped  through  a  range 
of  speeds.  In  a  second  type  of  test,  the  disk  angular  velocity  was  set  to  a  specific  speed 
and  the  measurement  volume  was  moved  through  a  range  of  radii  using  the  traverse. 

Care  was  taken  during  these  procedures  to  not  overload  the  LDV  system’s  photo 
multiplier  tubes  (PMT).  As  a  guideline,  FVA  software  high  voltage  settings  of  800-900 
worked  well,  with  a  low  laser  current  such  as  1 0  Amps. 

The  results  for  the  series  of  four  calibration  tests  from  1/25/02  are  shown  below. 
The  method  of  constant  angular  velocity  was  used  for  these  tests  because  it  could  be 
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executed  quickly.  The  angular  velocity  was  set  at  1 5  revolutions  per  second  to  avoid 
velocity  fluctuations.  Tangential  velocities  were  measured  in  a  range  of  radii  from 
-5  mm  to  127  mm.  These  combinations  of  radius  and  angular  velocity  enveloped  the 
range  of  expected  velocities  in  the  water  tunnel,  and  included  small  negative  velocities. 
The  radii  positions  were  stepped  sequentially  up  and  down,  to  detect  any  possible 
hysteresis.  No  hysteresis  was  detected  in  the  measurements,  and  the  results  were  very 
consistent. 

The  comprehensive  data  set  consisted  of  140  points.  A  linear  regression  was 
applied  to  the  data  to  return  a  slope  of  m  =  0.9681  and  an  intercept  of  b  -  -0.0352  m/s. 

An  uncertainty  analysis  was  performed  analytically  on  the  regression  equation,  and  a 
second  order  curve  fit  was  applied  to  the  uncertainty  to  obtain  an  estimate  of  the 
regression  uncertainty  Uu-regress  for  new  values  of  calibrated  velocity  uc: 

uc  =0.968  \uFVA-  0.0352  (3.4) 

Uu_regress  =  0.0507  +  (-2.4825  xl0^)wc  +  (2.4639  xl0~5)wc2  (3.5) 

The  linear  regression  with  uncertainty  bands  is  plotted  in  Figure  3.3  for  the  entire 
calibration  range.  The  linear  regression  plot  is  repeated  in  Figure  3.4  for  the  range  of 
expected  water  tunnel  velocities.  This  second  plot  also  includes  the  calibration  data 
points  with  error  bars. 
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Velocity  uFVA  Reported  by  FVA  (m/s) 


Figure  3.3:  LDV  calibration  plot  with  upper  and  lower  uncertainty  bands. 


Figure  3.4:  Calibration  plot  of  expected  water  tunnel  test  speed  range,  with  uncertainty 
bands  and  error  bars  on  the  data  points. 
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3.4  Submarine  Model  and  Mount 


The  submarine  shape  that  was  tested  was  identical  to  one  that  was  used  in  a 
number  of  experiments  at  the  LCC.  The  submarine  hull  shape  is  plotted  in  Figure  3.5. 
The  model  constructed  for  the  current  experiment  was  286.47  mm  (1 1 .278  in)  long  with  a 
diameter  of  25.4  mm  (1  in),  which  resulted  in  a  fineness  ratio  of  1 1 .278.  The  hull  form 
had  no  appendages  or  fins,  but  did  incorporate  a  propeller  hub  shape  on  the  end.  The 
model  was  constructed  of  two  pieces  of  stainless  steel  bar  stock.  The  hull  shape  was 
turned  down  on  a  CNC  lathe  in  the  MSU  Mechanical  Engineering  shop.  The  two  pieces 
were  screwed  together  to  form  the  complete  single  model. 
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Figure  3.5:  Submarine  model  drawing  showing  mounting  details. 
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The  model  was  mounted  along  the  test  section  centerline  by  means  of  a  single 
stainless  steel  strut.  The  strut  was  1.25  inches  wide  and  0.125  inches  thick,  which  was 
just  wide  enough  for  4-40  size  threads.  Two  4-40  size  flat  head  screws  were  used  to 
mount  the  model  on  one  end  of  the  strut.  Two  more  4-40  size  flat  head  screws  were  used 
to  mount  the  other  end  of  the  strut  onto  a  stainless  steel  mounting  plate.  This  mounting 
plate  was  in  turn  mounted  onto  the  inside  of  the  test  section  wall  by  two  flat  head  10-32 
screws,  which  were  threaded  into  two  holes  in  the  test  section  wall.  Nuts  were  put  on  the 
outside  to  further  secure  the  mount.  Silicon  was  used  between  the  mounting  plate  and  the 
test  section  wall  to  prevent  water  leakage.  Also,  all  exposed  countersunk  screw  heads  in 
the  mounting  plate  and  in  the  model  were  filled  with  epoxy  and  sanded  flush  with  400- 
grit  sandpaper.  The  model  was  mounted  so  that  the  nose  of  the  submarine  was 
approximately  6  body  diameters  (6  inches)  downstream  of  the  test  section  entrance.  The 
strut  was  mounted  on  the  submarine  model  in  the  position  corresponding  approximately 
to  the  position  of  the  strut  on  the  LCC  submarine  model.  The  leading  edge  of  the  strut 
was  located  at  x/LSUb  =  0. 183. 
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CHAPTER  IV 


EXPERIMENTAL  PROCEDURE  AND  RESULTS 

4.1  LDV  Configuration 

4.1.1.  LDV  Alignment  with  Water  Tunnel 
The  lightweight  traverse  that  was  used  with  the  LDV  had  three  axes  at  right 
angles.  The  x-direction  was  chosen  to  be  the  direction  of  flow;  the  y-direction  was  the 
vertical  axis  transverse  to  the  direction  of  flow,  and  the  z-direction  was  the  horizontal 
axis  transverse  to  the  direction  of  flow.  The  first  alignment  was  to  position  the  traverse 
system  so  that  the  x-axis  of  the  traverse  was  parallel  to  the  centerline  of  the  water  tunnel. 
This  was  achieved  by  hanging  plumb  lines  over  the  side  of  the  test  section  and  taping  a 
length  of  string  on  the  floor  between  the  two  points.  A  second  length  of  string  was  then 
taped  onto  the  floor  so  that  it  was  parallel  to  the  first  string,  with  about  20  inches  of 
distance  between  them.  The  traverse  was  then  positioned  so  that  plumb  lines  that  were 
hung  from  the  traverse  were  lined  up  on  the  second  string.  Next,  the  probe  was  aligned 
on  the  horizontal  plane  of  the  test  section  centerline.  To  do  this,  the  green  beams  were 
first  centered  on  the  outside  wall  of  the  test  section  in  the  vertical  direction  (y-direction) 
by  positioning  the  entering  laser  beams  and  the  exiting  reflected  beams  in  a  straight, 
horizontal  line  at  a  point  near  the  middle  of  the  test  section.  This  adjustment  was  made 
by  moving  the  probe  with  the  traverse.  The  probe  was  then  traversed  in  the  x-direction  to 
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the  end  of  the  test  section.  At  this  new  position,  the  beams  were  adjusted  in  the 
y-direction  by  raising  or  lowering  the  traverse  platform  on  that  end,  if  it  was  necessary. 
This  procedure  accounted  for  any  angular  displacement  the  water  tunnel  may  have  had 
from  the  true  horizontal. 

The  origin  of  the  experiment  coordinates  was  a  point  on  the  centerline  of  the 
outside  test  section  wall  and  7.5  body  diameters  downstream  of  the  submarine  body 
propeller  hub.  For  the  x-direction  reference,  the  blue  beams  were  touched  to  the  aft  end 
of  the  hub;  the  traverse  was  then  moved  downstream  7.5  body  diameters,  which  was  a 
distance  of  190.5  mm.  A  piece  of  masking  tape  with  a  small  black  dot  was  placed  at  this 
point,  so  that  when  the  measurement  volume  was  on  the  origin,  the  laser  beam 
“disappeared”  when  viewed  through  the  laser  safety  goggles. 

To  center  the  measurement  volume  onto  the  outside  wall  in  the  z-direction,  the 
following  procedure  was  used.  First  the  z-position  was  adjusted  until  the  laser  beams 
formed  a  fine  point  on  the  test  section  wall.  The  traverse  coordinates  were  then  set  to 
zero.  A  refinement  was  necessary  to  insure  that  the  middle  of  the  measurement  volume 
was  on  the  outer  wall  surface.  Using  the  green  beams,  the  measurement  volume  was 
stepped  in  increments  of  0.5  mm  in  the  z-direction,  starting  from  a  point  a  short  distance 
off  the  wall.  The  anode  current  reading  in  the  LDV  software  output  was  recorded  at  each 
point.  The  point  where  the  highest  anode  current  was  recorded  was  the  point  where  the 
center  of  the  measurement  volume  was  on  the  outer  surface  of  the  wall.  The  data  rate 
was  a  maximum  at  this  point. 
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It  was  found  that  this  alignment  procedure  had  to  be  performed  for  every  new 
experiment  because  the  water  tunnel  shifted  position  under  operation.  The  masking  tape 
with  the  black  dot  aided  in  locating  the  x  datum  easily.  The  y  and  z  datum  had  to  be 
reconfirmed  using  the  procedures  outlined  above. 

4. 1 .2  Laser  Beam  Refraction  Correction 

The  refraction  of  the  laser  beams  when  changing  media  must  be  accounted  for 
when  positioning  the  measurement  volume  within  the  test  section.  Fortunately,  the  index 
of  refraction  for  water  and  Plexiglas®  are  nearly  the  same,  so  that  only  one  refraction 
index  needed  to  be  considered.  The  overall  effect  of  the  change  in  media  from  air  to 
water  was  that  the  ratio  of  the  travel  distance  of  the  probe  volume  within  the  test  section 
to  the  travel  distance  of  the  probe  traverse  outside  of  the  test  section  was  equal  to  the 
index  of  refraction  of  the  new  medium.  The  index  of  refraction  for  water  is  4/3. 

4.2  Test  Section  Flow  Characteristics 

A  number  of  velocity  measurements  were  made  of  the  water  tunnel  test  section 
with  no  model  to  determine  the  flow  characteristics  of  the  water  tunnel.  Three  positions 
were  selected  for  velocity  profile  measurements:  x/L  =  0. 125,  x/L  =  0.492,  x/L  =  0.738; 
where  L  is  the  length  of  the  test  section.  At  each  of  these  positions,  several  velocity 
profiles  were  recorded  for  three  different  motor  input  speeds.  For  x/L  =  0.492,  the  motor 
input  speeds  were  500,  550,  and  600  rpm.  It  was  noted  during  this  time  that  the  water 
tunnel  was  being  excited  near  its  natural  frequency  at  motor  input  speeds  that  were 
multiples  of  200  rpm.  The  water  tunnel  shook  very  badly  at  600  rpm,  so  that  for  the 
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velocity  profiles  at  x/L  =  0. 125  and  x/L  =  0.738,  the  motor  input  speeds  were  500,  550, 
and  630  rpm. 

The  water  was  seeded  with  5-micron  silicon  carbide  during  these  tests  to  increase 
data  rates.  A  laser  current  of  16  Amps  was  used  throughout  the  tests.  The  PMT  high 
voltage  setting  had  to  be  increased  as  the  probe  measurement  volume  was  moved  further 
into  the  test  section.  Table  4.1  below  lists  the  optimum  high  voltage  setting  for  a  given 
range  of  z  for  a  moderately  seeded  flow. 


Table  4. 1 :  PMT  high  voltage  settings  for  test  section  regions. 


Range  in  z  (mm) 

PMT  High  Voltage 

15-40 

700 

45-75 

800 

80-100 

1000 

105-120 

1200 

The  plots  of  the  water  tunnel  velocity  profiles  are  shown  below  in  Figure  4.1 
through  Figure  4.4  for  each  of  the  motor  speeds.  The  presence  of  the  wall  boundary  layer 
is  noted  on  the  right  hand  side  of  each  plot.  It  is  evident  that  the  tunnel  speeds  are  more 
consistent  at  higher  velocities  and  near  the  entrance  to  the  test  section.  At  x/L  =  0.125,  all 
velocity  profiles  fall  right  on  top  of  each  other,  even  with  the  distinctive  kinks  in  the 
profile.  These  kinks  and  shifts  can  be  seen  in  other  profiles  as  well.  Figure  4.5  shows  the 
representative  velocity  profile  kinks  and  shifts  for  x/L  =  0.492  and  a  motor  input  speed  of 
500  rpm. 
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Figure  4.1 :  Velocity  profiles  in  the  test  section  at  motor  input  speed  of  500  rpm. 
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Figure  4.2:  Velocity  profiles  in  the  test  section  at  motor  input  speed  of  550  rpm. 
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Figure  4.3:  Velocity  profiles  of  the  water  tunnel  at  motor  input  speed  of  600  rpm. 
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Figure  4.4:  Velocity  profiles  of  the  water  tunnel  at  motor  input  speed  of  630  rpm. 
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Figure  4.5:  Velocity  profiles  taken  at  x/L  =  0.492,  500  rpm,  showing  typical  velocity 
shift  and  non-uniform  profile. 

During  the  initial  testing  of  the  LDV  system  on  the  water  tunnel,  numerous 
measurements  were  taken  at  a  fixed  point  within  the  test  section  to  determine  the  optimal 
LDV  system  configuration  settings.  During  this  time  the  motor  was  set  to  a  constant 
speed.  It  was  observed  that  the  mean  velocity  would  occasionally  shift  up  or  down  by 
approximately  0.5  m/s.  After  several  occurrences  of  this  phenomenon,  it  was  noted  that 
the  shifts  seemed  to  occur  suddenly  after  an  approximate  20-minute  period.  No  motor 
speed  fluctuations  were  noted  in  the  motor  controller  feedback.  Several  long  period  tests 
were  performed  to  try  to  capture  the  exact  period  of  the  velocity  shift,  but  these  proved 
inconclusive.  Besides  being  a  physical  phenomenon  of  the  tunnel,  the  velocity  shift 
could  have  been  caused  by  a  faulty  LDV  Bragg  cell  that  was  only  working  intermittently. 
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The  kinks  in  the  velocity  profiles  that  were  repeatable  from  data  set  to  data  set, 
even  when  the  velocity  shift  occurs,  were  most  likely  caused  by  the  local  differences  in 
the  optical  properties  of  the  test  section  wall.  As  the  measurement  volume  was  moved 
across  the  test  section,  the  actual  area  of  the  wall  that  the  beams  passed  through  changed. 
The  repeatability  of  these  kinks  at  the  same  r/R  positions  indicates  that  the  local  optical 
properties  of  the  test  section  were  affecting  the  crossing  angle  or  the  location  of  the  beam 
intersection  that  formed  the  measurement  volume.  The  standard  deviation  from  the  local 
mean  velocities  was  found  to  be  no  more  than  0.06  m/s  at  any  point  in  the  test  section, 
which  amounted  to  a  systematic  uncertainty  estimate  of  ±0.12  m/s  in  velocity 
measurement  due  to  the  optical  properties  of  the  test  section. 

Figures  4.6,  4.7,  and  4.8  show  the  local  turbulence  intensity  factor  urms/u  at 
x/L  =  0.492  and  motor  input  speeds  of  500,  550,  and  600  rpm.  The  typical  value  for  the 
freestream  turbulence  intensity  is  1.5  %,  which  is  slightly  higher  than  is  desirable. 

Figures  4.9,  4.10,  and  4.1 1  illustrate  the  typical  uncertainty  associated  with  the  velocity 
profile  measurements.  Figure  4.12  plots  the  same  data  as  Figure  4.1 1,  but  the  velocity 
axis  has  been  expanded  to  full  scale  to  give  a  better  sense  of  the  relative  magnitudes  of 
the  error.  The  error  in  the  velocity  measurements  includes  the  fossilized  systematic  error 
from  the  calibration,  the  systematic  uncertainty  estimate  due  to  the  test  section  optics,  and 
the  random  uncertainty  of  each  of  the  point  measurements.  The  radial  position  errors  are 
due  to  uncertainties  in  the  traverse  calibration,  the  finite  length  of  the  measurement 
volume,  and  the  uncertainty  in  the  test  section  radius  R. 
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Figure  4.6:  Local  turbulence  intensity  factors  at  x/L  =  0.492,  500  rpm. 
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Figure  4.7:  Local  turbulence  intensity  factors  at  x/L  =  0.492,  550  rpm. 
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Figure  4.8:  Local  turbulence  intensity  factors  at  x/L  =  0.492,  600  rpm. 
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Figure  4.9:  Typical  velocity  profile  with  error  bars,  x/L  =  0.492,  500  rpm. 


6 


(z-zC|)/R 


Figure  4.10:  Typical  velocity  profile  with  error  bars,  x/L  =  0.492,  550  rpm. 
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Figure  4. 1 1 :  Typical  velocity  profile  with  error  bars,  x/L  =  0.492,  600  rpm. 
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Figure  4.12:  Typical  velocity  profile  with  error  bars,  x/L  -  0.492,  600  rpm,  repeated  with 
full  velocity  scale  to  show  relative  magnitude  of  errors. 


Motor  Input  RPM 

Figure  4.13:  Mean  test  section  velocity  at  x/L  =  0.125  versus  motor  input  speed.  Error 
bars  are  shown  on  the  data  points,  and  uncertainty  bands  are  shown  around  the  curve  fit. 
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Figure  4.13  demonstrates  the  linear  relationship  between  motor  input  speed  and 
test  section  velocity  as  measured  atx/L=  0.125.  The  error  bars  are  shown  on  the  three 
data  points  at  500,  550,  and  630  rpm.  The  uncertainty  bands  that  are  shown  apply  to  the 
curve  fit,  and  are  derived  from  the  “classical  approach”  to  regression  uncertainty  as 
described  in  Coleman  and  Steele  1 .  A  least  squares  linear  regression  to  this  data  is  given 
by  the  equation: 

t/.  =0.0161  Q  + 0.5241,  (4.1) 

where  O  is  the  motor  input  speed  in  revolutions  per  minute  and  U„  denotes  the  test 
section  velocity  at  x/L  =  0.125.  The  uncertainty  associated  with  the  regression  is  given 
by  the  equation: 

Uv,.regress  =  0.040878 

The  uncertainty  of  U„  from  the  curve  fit  could  be  reduced  by  a  factor  of  six  by  taking 
velocity  profiles  at  nine  different  motor  speeds  rather  than  just  using  three  motor  speeds. 

4.3  Submarine  Model  Wake  Experiments 

For  the  submarine  wake  measurements,  the  tunnel  motor  speed  was  set  to  500 
rpm.  Velocity  profiles  were  taken  at  several  x/d  stations  behind  the  model,  where  x/d=  0 
is  the  end  of  the  model.  Velocity  profiles  were  measured  at  x/d=  10, 13,  16,  19,  22,  25, 
and  28.  The  time  required  to  complete  one  complete  set  of  profiles  was  approximately  5 
hours.  Most  of  that  time  was  spent  measuring  the  wake  velocity  deficits,  which  had 
accompanying  low  data  rates.  Optimum  LDV  settings  for  the  wake  profile  measurements 


1  (0-560) 

3+  8600 


(4.2) 
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were  found  through  trial  and  error.  The  laser  current  setting  for  the  tests  was  1 6  Amps. 
The  PMT  high  voltage  setting  was  1050,  with  no  seeding  in  the  tunnel.  A  drape  was  used 
over  the  test  section  and  traverse  to  block  light  noise  from  the  overhead  lights.  A  dark 
drape  was  also  hung  vertically  behind  the  test  section  to  reduce  glare  off  of  the  wall. 

The  wake  velocity  profiles  and  the  turbulence  intensity  factors  from  the  five  days 
of  testing  are  shown  in  Figure  4.14  through  Figure  4.23.  The  point  spacing  for  the 
5/16/02  through  5/22/02  data  was  3.33  mm,  which  corresponded  to  a  traverse  movement 
of  2.5  mm.  For  the  remaining  days,  the  point  spacing  was  reduced  to  2  mm 
(corresponding  to  a  traverse  movement  of  1 .5  mm)  and  the  measurement  of  the  velocity 
profile  at  x/d  =  28  was  discontinued.  The  point  spacing  was  reduced  to  increase  the 
number  of  points  in  the  wake.  The  effect  of  measurement  volume  overlap  was  ignored. 
The  measurements  at  x/d  =  28  were  discontinued  for  an  operational  time  savings  of 
nearly  half  an  hour. 
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Figure  4.14:  Mean  velocity  profiles  from  5/16/02. 
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Figure  4.15:  Turbulence  intensity  from  5/16/02. 
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Figure  4.16:  Mean  velocity  profiles  from  5/17/02. 
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Figure  4.17:  Turbulence  intensity  from  5/17/02. 
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Figure  4.18:  Mean  velocity  profiles  from  5/22/02. 
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Figure  4.19:  Turbulence  intensity  from  5/22/02. 
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Figure  4.20:  Mean  velocity  profiles  from  5/30/02. 
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Figure  4.21 :  Turbulence  intensity  from  5/30/02 


The  mean  velocity  profiles  are  asymmetric  and  are  generally  not  centered  along 
the  measured  centerline  of  the  test  section.  The  measurements  of  the  wake  beyond  the 
centerline  were  characterized  by  higher  turbulence  intensity  and  lower  data  rates.  The 
peculiar  behavior  of  the  velocity  profiles  suggests  that  the  model  was  not  aligned 
properly  with  the  flow.  After  the  model  was  mounted,  the  accuracy  of  alignment  was 
checked  and  found  to  be  within  0.5°.  However,  previous  studies  had  discovered  that  the 
test  section  becomes  distorted  when  pressurized.  According  to  measurements  by 
Lindsey  6,  the  test  section  bulged  0.009  in  (0.23  mm)  when  pressurized  to  50  psi.  It  is 
possible  that  the  deflection  of  the  test  section  while  under  pressure  caused  a  misalignment 
of  the  model  with  the  mean  flow.  Also,  it  was  observed  during  operation  that  the  strut 
would  occasionally  deflect  from  side  to  side.  These  strut  deflections  may  have 
contributed  to  the  asymmetry  of  the  wake  profiles. 

The  wake  centerline  drift  seen  in  the  plots  represents  an  offset  of  approximately 
0.5°.  During  later  measurements  of  the  test  section  wall  surface  irregularities,  it  was 
discovered  that  the  traverse  x-axis  was  misaligned  by  0.7°  in  the  same  direction  as  the 
centerline  shift.  This  misalignment  did  not  affect  the  final  drag  calculations  as  the  wake 
centerline  was  corrected  in  the  curve  fitting  process.  The  uncertainty  that  the  angular 
misalignment  imposed  on  the  LDV  velocity  measurements  was  approximately  ±0.007 
m/s. 

Figures  4.24,  4.25,  and  4.26  illustrate  the  typical  uncertainty  associated  with  the 
velocity  measurements.  The  data  is  taken  from  the  series  of  measurements  from  5/22/02, 


56 


and  at  x/d  =  10,  16,  25.  The  velocity  and  position  error  sources  are  taken  to  be  the  same 
as  those  used  in  the  test  section  characteristic  profiles. 


Figure  4.24:  Velocity  profile  with  typical  uncertainties,  5/22/02,  x/d  -  10. 
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Figure  4.25:  Velocity  profile  with  typical  uncertainties,  5/22/02,  x/d  =  16. 
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Figure  4.26:  Velocity  profile  with  typical  uncertainties,  5/22/02,  x/d-  25. 


4.3.2  Curve  Fits  to  Wake  Profile  Data 


The  wake  profile  data  was  fit  using  the  numerical  form  of  equation  (2.59)  so  that 
the  drag  coefficient  of  the  hull  form  could  be  calculated.  The  curve  fits  were  done  in 
MathCAD  using  the  “genfit”  function.  The  genfit  function  finds  the  least  squares  fit  of 
an  arbitrary  function  to  a  given  data  set.  First,  the  wake  velocity  profile  data  z  and  u  were 
read  from  a  file  and  the  Gaussian  curve  fit  function  was  applied  in  the  form: 

— (z-*o)2 

u(z)  =  U0-u0e  1,2  (4.3) 

where  Uo,  uo,  zo,  and  b  were  generated  from  the  curve  fit.  zo  is  the  wake  centerline 
location,  b  is  the  wake  width  parameter,  uo  is  the  velocity  scale  factor,  and  Uo  is  the  mean 
velocity  outside  of  the  wake.  Since  the  curve  fit  process  did  not  include  the  (In  2)  factor, 

the  wake  width  parameter  was  corrected  after  the  curve  fit  by  dividing  by  In  2) .  This 


curve  fitting  process  was  repeated  for  each  x/d  location. 

After  the  wake  curve  fit  was  completed,  the  fits  of  length  scale  and  velocity  scale 
were  performed  to  obtain  the  coefficients  cu  and  Ci  that  were  introduced  in  equations 
(2.61)  and  (2.62).  The  set  of  data  x/d,  uo,  and  b  were  again  used  with  the  genfit  function 
to  fit  the  velocity  defect  function: 
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The  curve  fit  yielded  c„  and  xf/d.  The  x</d  value  was  then  used  in  the  function  fit  for  the 
wake  width  parameter  to  find  c/ : 


b(x )  ( x  —  x o  Y* 

- =  c,  - 

d  l  d  ) 


(4.5) 


In  order  to  use  these  power  law  fits,  the  measurements  had  to  be  made  in  a  region 
of  the  wake  where  the  velocity  deficit  and  width  parameters  follow  this  behavior.  In 
experiments  conducted  by  Chevray  8,  the  wake  of  a  6:1  prolate  spheroid  achieved 
similarity  at  x/d  ~  10 .  All  measurements  in  the  current  experiments  were  made  at  an 
x/  d>  10  ,  and  since  the  submarine  model  had  a  larger  fineness  ratio  than  that  of 
Chevray’ s  spheroid,  the  wake  of  the  submarine  model  should  have  been  self-similar  at  all 
the  measured  positions.  The  self-similarity  of  the  submarine  wake  profiles  is  illustrated 
in  Figure  4.27  and  Figure  4.28.  The  wake  velocity  profiles  are  clearly  self-similar  and  fit 
the  Gaussian  function,  although  the  region  outside  the  wake  is  noisy  in  the  5/22/02  data. 
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Figure  4.27:  Similarity  profiles  for  5/22/02 
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Figure  4.29  illustrates  the  variation  of  the  freestream  velocity  Uo  outside  the  wake 
as  x/d  changes.  The  general  trend  is  that  U0  increases  as  x/d  increases;  however,  there  are 
local  variations  in  this  trend  that  were  not  repeatable.  Figure  4.30  plots  the  changes  in 
nondimensionalized  velocity  deficit  uo/Uo  with  x/d.  The  results  are  similar,  if  not 
consistent.  Figure  4.31  plots  the  changes  in  the  wake  width  parameter  b/d  over  x/d.  The 
fits  are  not  very  repeatable. 


Figure  4.29:  Freestream  velocity  Uo  outside  of  wake. 
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Figure  4.30:  Velocity  deficit  uq/Uq. 


Figure  4.31:  Wake  width  parameter  b/d. 
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4.3.3  Drag  Coefficient  of  Model 

The  drag  coefficient  may  be  calculated  from  each  day’s  data  through  the  use  of 


equation  (2.65).  The  drag  coefficient  can  then  be  corrected  for  the  pressure  gradient  that 
exists  in  the  test  section.  This  pressure  gradient  was  measured  by  Lindsey  6  for  various 
tunnel  speeds.  For  a  tunnel  speed  setting  of  500  rpm,  the  pressure  gradient  is  -0.0058 
psi/in,  with  an  uncertainty  of  ±0.0009  psi/in.  The  drag  coefficient  due  to  the  pressure 
gradient  is  given  by: 


Cd,p  ~ ' 


8  F. 


pTufU^ 


(4.6) 


where  Fp  is  the  reaction  force  on  the  body  due  to  the  pressure  gradient,  and  U„  is  the 
freestream  velocity  as  given  by  equation  (4.1).  The  pressure  reaction  force  Fp  is  found  by 
multiplying  the  pressure  gradient  with  the  volume  of  the  body: 

Fp  (4J) 


This  is  the  standard  form  for  the  so-called  “buoyancy  drag”  as  reported  in 
AGARDograph  109  9.  The  resulting  pressure  force  in  any  case  is  approximately  0.04  lbf, 
and  the  resulting  contribution  to  the  drag  coefficient  is  A CD  p  ~  0.0 1 0 ,  where  the 


reference  velocity  is  U„=  8.575  m/s. 

The  results  of  the  fits  are  reported  in  Table  4.1.  As  can  be  seen,  the  results  were 
varied  and  unrepeatable.  The  range  of  drag  coefficients  were  from  0.122  to  0.351.  The 
uncertainty  in  the  drag  coefficient  due  to  the  uncertainty  in  the  pressure  gradient  is 
±0.002.  The  drag  coefficient  was  expected  to  be  approximately  0.2.  This  expectation 
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was  based  on  previous  experiments  with  the  same  model  at  the  LCC  and  on 
computational  studies  by  Cash  l0.  The  result  that  most  nearly  approaches  the  expected 
value  is  from  the  5/16/02  data  set;  yet  it  is  still  15.5%  higher  than  the  expected  value. 


Table  4.2:  Curve  fit  results  for  each  day’s  experiment. 


Experiment 

Date 

x«/d 

cu 

Cl 

CD 

(neglecting  pressure 
gradient) 

cD 

5/16/02 

2.501 

0.566 

0.192 

0.241 

0.231 

5/17/02 

4.377 

0.49 

0.253 

0.361 

0.351 

5/22/02 

5.55 

0.475 

0.169 

0.157 

0.147 

5/30/02 

3.595 

0.565 

0.208 

0.283 

0.283 

6/4/02 

3.511 

0.594 

0.139 

0.132 

0.122 

The  test  conditions  varied  from  day  to  day,  which  may  account  for  some  of  the 
scatter  in  the  results.  The  data  from  5/16/02  and  5/17/02  were  some  of  the  first 
measurements  to  be  taken  of  the  wake  profiles.  During  these  tests,  the  tunnel  fill  valve 
was  not  left  open.  The  tunnel  was  pressurized  to  50  psi,  and  then  the  fill  valve  was 
turned  off.  Normally,  the  tunnel  would  remain  at  constant  pressure.  However,  it  had 
been  noticed  that  the  water  tunnel  leakage  had  increased  during  the  first  tests.  After  the 
tests  were  completed  on  5/22/02,  it  was  noticed  that  the  static  pressure  in  the  tunnel  had 
dropped  to  25  psi.  This  had  not  been  observed  before.  It  was  decided  to  keep  the  fill  line 
on  for  all  the  tests  thereafter  to  keep  the  tunnel  pressurized  at  50  psi.  The  velocity 
profiles  measured  on  5/30/02  and  6/4/02  were  taken  under  a  constant  pressure  of  50  psi. 
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and  demonstrated  distinct  differences  from  the  profiles  of  the  three  previous  days, 
although  the  input  speed  of  the  motor  was  the  same  for  all  tests.  This  difference  in 
behavior  is  illustrated  in  Figure  4.32. 


Figure  4.32:  Comparison  of  wakes  profiles  at  x/d=  16  for  all  five  data  sets. 

After  the  wake  measurements  on  6/4/02,  a  study  was  performed  to  test  the  effect 

of  decreasing  pressure  on  the  flow  velocity  within  the  test  section.  The  LDV 

measurement  volume  was  placed  in  the  freestream  and  velocity  measurements  were  taken 

continuously  while  the  tunnel  static  pressure  was  allowed  to  bleed  off  from  50  psi  to  25 

psi.  A  plot  of  the  time  histories  for  the  velocity  and  pressure  readings  is  shown  in 

Figure  4.33.  The  results  were  unexpected — the  velocity  fluctuated  wildly  with  no 
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apparent  relation  to  the  pressure  drop.  A  similar  test  was  performed  with  the  tunnel  fill 
valve  left  open,  so  that  the  tunnel  remained  at  a  constant  static  pressure  of  46  psi.  The 
results  are  shown  in  Figure  4.34.  Except  for  a  sudden  drop  in  velocity  observed  near  the 
beginning  of  the  test,  the  velocity  was  fairly  constant  for  nearly  an  hour.  This  test  came 
to  an  abrupt  end  when  the  third  comer  of  the  water  tunnel  catastrophically  failed  under 
pressure.  The  failure  prevented  further  study  into  the  effects  of  decreased  tunnel  pressure 
and  the  velocity  shifts  that  were  noted  earlier. 


Figure  4.33:  Test  section  freestream  velocity  behavior  over  a  period  of  decreasing  static 
pressure. 
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Figure  4.34:  Test  section  freestream  velocity  behavior  over  a  period  of  constant  static 
pressure  of  46  psi. 

The  profiles  of  5/16/02  and  5/17/02  were  taken  under  similar  conditions.  The 
tunnel  water  was  seeded  with  5-micron  silicon  carbide,  which  was  found  during  the  test 
section  studies  to  improve  data  rates.  However,  with  the  presence  of  the  model  and  its 
wake,  the  seeding  was  found  to  be  detrimental.  The  water  was  drained  after  the  tests  on 
5/17/02,  and  clean  water  was  used  on  5/22/02  and  every  day  thereafter,  with  much  better 
readings  within  the  wake.  The  data  taken  on  5/16/02  and  5/17/02  were  not  originally 
intended  to  be  analyzed  and  reported,  but  they  were  analyzed  after  the  failure  of  the 
tunnel  prevented  any  further  testing  to  see  if  any  other  clues  to  the  inconsistent 
measurements  could  be  found.  The  tunnel  pressure  was  not  observed  during  the  tests  on 
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5/16/02  or  5/17/02.  Another  difference  in  operating  conditions  was  the  absence  of 
shading  from  the  overhead  lighting,  which  was  present  from  5/22/02  and  afterwards. 

After  5/22/02,  the  fill  line  was  left  on  during  the  tests  to  maintain  the  static 
pressure  in  the  tunnel.  The  tests  of  5/30/02  and  6/4/02  were  run  with  constant  static 
pressure  of  approximately  50  psi,  consistent  ambient  lighting,  and  clean  water.  It  was 
expected  that  these  later  tests  would  yield  the  most  acceptable  data.  The  drag  coefficients 
from  these  later  tests  were  0.283  and  0.122,  which  are  approximately  40%  higher  and 
lower  than  the  expected  drag  coefficient  of  approximately  0.2. 

To  further  explore  the  effect  of  the  different  operating  conditions,  the  data  was 
divided  into  two  groups.  The  first  three  days’  data  represented  the  low  pressure  operating 
condition,  and  the  last  two  days’  data  represented  the  high  pressure  operating  condition. 
The  velocity  profile  similarity  parameters  of  the  individual  days’  data  were  averaged  to 
form  composite  velocity  profiles  representing  the  low  pressure  and  the  high  pressure 
operating  conditions.  The  power  law  fits  were  then  applied  to  these  two  sets  of  data,  as 
shown  in  Figure  4.35  through  Figure  4.38. 
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Figure  4.35:  Comparison  of  averaged  velocity  deficit  parameter  from  the  first  three  days’ 
experiments  to  the  power  law. 


Figure  4.36:  Comparison  of  averaged  wake  width  parameter  from  the  first  three  days’ 
experiments  to  the  power  law. 
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Figure  4.37:  Comparison  of  averaged  velocity  deficit  parameter  from  the  last  two  days’ 
experiments  to  the  power  law. 


Figure  4.38:  Comparison  of  averaged  wake  width  parameter  from  the  last  two  days’ 
experiments  to  the  power  law. 


The  fits  for  uq/Uq  and  b/d  for  the  composite  decreased  pressure  measurements  are 
very  close  to  the  forms  that  should  be  expected.  For  the  composite  constant  pressure 
measurements,  the  fit  for  uo/Uo  is  fairly  close  to  proper  form,  but  the  fit  for  b/d  is  erratic. 
The  second  fit  parameters  xo/d,  cu,  ci  and  Cp  are  tabulated  in  Table  4.3.  The  drag 
coefficient  results  are  very  close  to  each  other,  and  are  in  good  agreement  with  the 
expected  results  from  the  LCC  experiments  and  the  computational  results  of  Cash  10.  The 
computational  results  and  the  averaged  experimental  results  for  this  hull  shape  are  plotted 
in  Figure  4.39.  The  individual  and  averaged  experimental  results  are  shown  for  two 
different  temperatures,  because  the  viscosity  of  water  is  very  sensitive  to  the  water 
temperature.  The  temperature  of  20°C  represents  room  temperature  conditions.  The 
temperature  of  30°C  is  included  because  the  water  always  seemed  to  be  rather  warm.  No 
temperature  of  the  water  was  taken,  so  the  exact  Reynolds  number  is  unknown. 
Regardless,  the  averaged  results  are  in  good  agreement  with  the  trend  of  the 
computational  results.  This  is  again  shown  in  Figure  4.40,  which  shows  the  averaged 
results  along  with  Cash’s  computational  results  for  three  different  submarine  lengths. 


Table  4.3:  Curve  fit  results  from  composite  data  sets. 


Experiment 

Date 

Xo/d 

Cu 

Cl 

CD 

(neglecting  pressure 
gradient) 

CD 

Composite 
5/16,  5/17, 
5/22 

6.930 

0.445 

0.195 

0.195 

0.185 

Composite 
5/30,  6/4 

3.552 

0.580 

0.174 

0.201 

0.191 
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Figure  4.39:  Comparison  of  experimental  results  with  computational  results  of  Cash. 


Small  Model  (0.94  ft),  CFD  -b-  LCC  Model  (22.964  ft),  CFD 
^  Full  Scale  (383.5  ft),  CFD  x  MSU  Model  (0.94  ft),  20-C  Exp. 

o  MSU  Model  (0.94  ft),  30-C  Exp. 


Figure  4.40:  Averaged  experimental  results  for  MSU  model  and  computational  results  of 
Cash  for  three  different  submarine  lengths. 
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CHAPTER  V 


CONCLUSIONS  AND  RECOMMENDATIONS 

5.1  Conclusions 

The  water  tunnel  was  repaired  after  numerous  problems  developed  after  its  last 
operational  period.  The  comers  of  the  water  tunnel  and  the  drive  section  were  repainted 
with  several  coats  of  protective  paint.  The  drive  system  of  the  tunnel  was  modified  to 
solve  a  longstanding  problem  of  bearing  bum-out.  After  these  renovations  were 
completed,  the  water  tunnel  was  ran  nearly  daily  for  a  period  of  two  months. 

A  calibration  device  was  constructed  for  the  Aerospace  Engineering  Department’s 
new  LDV  system.  An  uncertainty  analysis  of  the  calibration  process  was  completed. 
Once  calibrated,  the  LDV  system  was  found  to  be  accurate  to  within  0.75%  for  the 
velocity  range  of  5-10  m/s. 

The  water  tunnel  test  section  velocity  profiles  were  measured  with  high  accuracy 
using  the  LDV  system.  These  velocity  profiles  are  the  first  detailed  measurements  of  the 
tunnel  flow  quality.  Generally,  the  turbulence  intensity  of  the  freestream  is  1.5%.  A 
highly  accurate  tunnel  speed  curve  was  developed  so  that  for  a  given  motor  input  speed, 
the  freestream  velocity  at  x/L  =  0.125  is  accurately  known.  Also,  it  was  discovered  that 
the  velocity  profiles  of  the  undisturbed  freestream  had  repeatable  kinks  at  a  given  x/L. 
These  kinks  are  believed  to  have  been  caused  by  local  optical  properties  in  the 
Plexiglas®  test  section.  During  repeated  measurements,  it  was  noticed  that  the 
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freestream  velocity  shifted  occasionally  by  almost  0.5  m/s.  The  cause  or  possible  period 
of  this  shift  was  not  determined.  It  is  possible  that  static  pressure  drop  due  to  leakage 
was  a  cause,  or  that  the  LDV  system  had  a  faulty  Bragg  cell.  The  effects  of  tunnel 
pressure  on  velocity  measurements  were  studied  later  after  the  measurements  of  the 
submarine  wake  had  been  completed.  It  was  during  these  studies  that  the  water  tunnel 
ruptured  under  pressure. 

An  1 1.278-inch  long  submarine  model  was  mounted  along  the  tunnel  centerline. 
The  wake  velocity  profiles  were  measured  at  10,  13,  16,  19,  22,  25,  and  28  body 
diameters  downstream  of  the  stem  of  the  model  on  five  different  days.  The  velocity 
profiles  exhibited  asymmetry  and  a  centerline  shift.  The  velocity  profiles  were  fit  to  a 
Gaussian  function,  and  self-similarity  parameters  were  obtained.  The  wake  velocity 
profiles  were  found  to  be  self-similar  beyond  x/d  =  10.  Power  law  fits  were  applied,  and 
the  drag  coefficient  of  the  submarine  model  was  calculated.  The  individual  data  sets 
were  very  inconsistent,  with  Co  ranging  from  0.122  to  0.351.  However,  when  the  data 
sets  were  combined  into  two  sets,  one  representing  measurements  under  decreasing 
tunnel  pressure  and  one  representing  measurements  under  a  high  constant  tunnel 
pressure,  the  drag  coefficients  were  0.185  and  0.191,  respectively.  Although  no  firm 
conclusions  can  be  made  about  the  drag  coefficient  results,  it  is  believed  that  the 
distortion  of  the  tunnel  test  section  under  high  pressure  may  have  caused  the  submarine 
model  to  become  misaligned  with  the  flow.  This  could  have  caused  premature  boundary 
layer  separation  and  a  higher  drag  coefficient  as  compared  to  the  low-pressure  results. 
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5.2  Recommendations 


The  third  comer  of  the  water  tunnel  must  be  repaired  or  replaced  before  the  tunnel 
is  operated  again.  There  is  a  2-ft  long  rupture  in  the  fiberglass,  as  well  as  several  bubbles 
in  the  inner  coating.  It  is  likely  that  water  penetrated  into  the  fiberglass,  causing  the 
layers  to  debond  and  thus  weaken  the  structural  integrity  of  the  section.  Since  the 
bubbles  are  located  throughout  the  damaged  section,  it  is  likely  that  water  penetrated 
everywhere  and  the  existing  comer  must  be  scrapped.  Discussion  with  Raspet  Flight 
Research  Laboratory  personnel  who  performed  the  repainting  work  revealed  that  the  base 
coats  did  not  adhere  well  to  the  tunnel  parts.  In  light  of  this,  all  the  comers  and  the  drive 
section  should  be  inspected  for  water  damage. 

Besides  repair  work,  the  one  improvement  that  would  greatly  expand  the  research 
capability  of  the  water  tunnel  would  be  to  construct  a  rectangular  test  section  so  that  the 
LDV  is  not  restricted  to  velocity  measurements  on  the  horizontal  plane  of  the  centerline. 
This  effect  can  be  achieved  by  building  a  transparent  box  around  the  existing  cylindrical 
test  section,  and  filling  it  with  water.  A  ray-tracing  program  may  be  needed  if  this  is 
done,  to  account  for  four  changes  of  media  instead  of  only  two. 

A  number  of  smaller  modifications  can  be  made  to  the  water  tunnel  to  expand  its 
capability.  A  thermometer  to  measure  the  temperature  of  the  water  inside  the  tunnel 
would  be  useful  for  accurately  determining  Reynolds  number  for  experiments.  It  is  also 
recommended  that  the  tunnel  drive  system  be  improved  for  smoother  operation.  It  was 
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found  that  the  water  tunnel  vibrated  excessively  at  operating  speeds  higher  than  550  rpm. 
This  limits  the  test  section  velocity  to  about  9.7  m/s  (31.8  ft/s). 

A  more  rigid  and  adjustable  platform  for  the  LDV  traverse  would  also  expand  the 
measurement  capability  of  the  system.  The  current  traverse  can  only  cover  half  the 
length  of  the  test  section  at  one  time.  If  the  traverse  were  mounted  on  a  platform  that 
could  be  easily  rolled  while  maintaining  alignment  with  the  tunnel  test  section, 
measurements  could  be  made  anywhere  in  the  test  section  during  a  single  period  of 
testing.  A  rail-mounted  platform  is  envisioned  for  this  purpose,  with  the  angular 
alignment  of  the  rails  being  adjustable  to  remain  parallel  with  the  tunnel  test  section  wall. 

Finally,  the  LDV  system  needs  to  be  examined  to  determine  that  it  is  in  good 
working  order.  The  sudden  velocity  shifts  observed  during  the  experiment  may  be  due  to 
the  intermittent  operation  of  the  LDV  Bragg  cell.  Successive  measurements  of  a  constant 
LDV  calibration  disk  speed  over  a  long  period  of  time  may  show  whether  the  velocity 
shifts  in  the  data  were  caused  by  an  anomaly  in  the  LDV  system  or  in  the  water  tunnel. 
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